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The  coupled-cluster  equation-of-motion  (CC-EOM)  method,  formally  equivalent 
to  the  Fock-space  coupled-cluster  (FSCC)  and  coupled-cluster  linear-response  (CCLRT) 
methods  for  electron  detachment  energies,  offers  distinct  advantages  over  the  single- 
reference  coupled  cluster  (SRCC)  approach  when  several  electron-detached  states  are 
to  be  calculated.  The  method  is  capable  of  providing  transition  energies  for  several 
states  at  a  computational  cost  less  than  that  of  an  SRCC  calculation  for  one  of  them. 

The  equation-of-motion  ionization  potential  (EOMIP)  method  is  implemented 
using  a  Davidson-type  algorithm  for  extracting  selected  roots  from  a  large  symmetric 
matrix,  making  it  possible  to  selectively  search  for  the  diflacult-to-characterize  "shake- 
up"  satellite  peaks  often  found  in  inner  valence  structure,  and  several  approximations 
to  the  ionization  intensities  are  available,  making  it  possible  to  simulate  electron 
detachment  spectra. 
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The  EOMIP  method  is  caUbrated  against  experimental  results  for  several  small 
molecules,  where  the  outer-  and  inner-valence  regions  are  explored,  with  special  at- 
tention to  shake-up  transitions.  Results  compare  favorably  with  the  coupled-cluster 
singles  and  doubles  (CCSD)  approach  in  most  cases.  The  method  is  also  applied 
to  the  interpretation  of  the  electron  detachment  spectra  of  the  linear  carbon  cluster 
anions,  C2  through  Cg  . 
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CHAPTER  1 
INTRODUCTION 

U Correlated  Methods  in  Quantum  Chemistry 

In  the  early  days  of  quantum  chemistry,  the  only  ab  initio  method  available  was 
the  self-consistent  field  (SCF)  approximation.  This  approximation  to  the  electronic 
Hamiltonian  neglects  the  interaction  between  individual  electrons  (electron  correla- 
tion), placing  each  one  in  an  average  field  of  all  the  other  electrons.  The  approach 
yields  a  simple  molecular  orbital  picture  of  atoms  and  molecules  and,  in  many  cases, 
yields  qualitatively  good  results  for  the  energy  differences  between  states  of  a  molecule 
or  reactive  system.  The  molecular  orbital  (MO)  picture  can  provide  simple,  useful 
concepts  of  physical  processes,  such  as  reaction,  excitation,  and  ionization.  But  in 
many  cases  this  simple  picture  is  not  even  qualitatively  correct.  And  there  are  many 
physical  processes  for  which  it  is  not  capable  of  providing  any  explanatory  tools. 

The  case  in  point  is  the  photoelectron  process.  When  electromagnetic  radiation 
in  the  proper  frequency  range  strikes  an  atom  or  molecule,  an  electron  will  be  ex- 
cited into  the  free-electron  continuum.  Photoelectron  spectroscopists  then  observe 
the  absorption  spectrum  or  the  spectrum  of  the  free  electron  kinetic  energies.  Koop- 
mans'  theorem[l]  is  an  SCF-based  approach  which  states  that  the  detached  electron 


is  removed  from  a  specific  MO,  and  that  the  energy  required  is  approximately  the 
energy  of  that  orbitaJ.  This  approximation  gives,  in  many  cases,  a  qualitatively  useful 
prediction  of  the  principal  peaks  of  the  photoelectron  spectrum,  but  it  is  incapable 
of  making  any  prediction  for  the  satellite  peaks  which  are  close  in  energy  to  these 
principal  peaks.  These  so-called  shake-up  peaks  are,  in  the  MO  picture,  the  result 
of  a  mixing  of  the  primary  hole-state  process  (ejection  of  an  electron  from  a  previ- 
ously occupied  orbital)  with  a  two-hole,  one-particle  process  (simultaneous  ejection 
from  an  occupied  orbital  and  excitation  of  an  electron  from  an  occupied  to  a  virtual 
orbital).  In  order  to  predict  these  states,  it  is  necessary  to  improve  upon  the  SCF 
approximation  —  to  introduce  electron  correlation. 

The  configuration  interaction  (CI)  method  has  been  widely  used  to  introduce 
correlation.  In  the  single-reference  formulation,  it  is  simple  in  concept.  Moving  an 
electron  from  an  "occupied"  orbital  (creating  a  "hole" )  to  an  unoccupied  or  "virtual" 
orbital  (creating  a  "particle")  defines  a  single  excitation.  Doing  this  twice  defines  a 
double  excitation,  etc.  The  correlated  wavefunction  is  defined  as  a  linear  combina- 
tion of  the  reference  function  and  all  possible  single,  double,  etc.,  excitations  from  it. 
The  coefficients  are  the  eigenvectors  of  the  matrix  representation  of  the  full  electronic 
Hamiltonian  in  the  basis  of  all  these  functions.  The  lowest-energy  solution  yields 
the  correlated  energy  and  wavefunction  corresponding  to  the  SCF  reference,  and  the 
higher  roots  correspond  to  excited  electronic  states.  The  "full  CI"  solution,  as  it  is 
called,  is  the  best  possible  wavefunction  for  a  given  basis  set  and  satisfies  the  varia- 
tional principle.  In  practice,  however,  the  full  CI  problem  can  rarely  be  solved.  As 


the  number  of  basis  functions  (M)  increases,  the  number  of  possible  excited  determi- 
nants increases  as  a  function  approaching  M\.  Standard  implementations  therefore 
truncate  the  expansion  space,  typically  to  only  a  subset  of  single  and  double  excita- 
tions from  the  reference.  But  this  truncation  introduces  a  particularly  troublesome 
error  —  the  energy  of  the  system  no  longer  scales  properly  with  the  size  of  the  system 
being  studied.  The  percentage  error  in  the  calculated  correlation  energy  increases 
rapidly  as  the  size  of  the  system  increases,  more  than  doubUng  each  time  the  size  of 
the  system  doubles.  The  error  in  the  energy  is  reflective  of  errors  in  the  wavef unction, 
and  calculations  of  observable  properties,  which  depend  on  the  wavefunction,  will  not 
scale  properly  with  size.  Tj'pically,  these  truncated  CI  methods  cannot  be  applied 
reliably  to  systems  larger  than  benzene. 

The  many-body  perturbation  theory  and  coupled-cluster  (CC)[2]  methods  have 
been  developed  specifically  to  avoid  the  deficiencies  of  the  CI  method.  These  methods 
are  "size-extensive";  they  scale  properly  with  the  size  of  the  system.  Implementations 
of  these  methods  have  become  increasingly  available  in  the  past  10  years  and  have 
been  shown  to  be  very  powerful,  capable  of  predicting  experimental  results  for  many 
types  of  molecular  systems  to  a  high  degree  of  accuracy.  Second-order  many-body 
perturbation  theory  (MBPT-2)  calculations,  which  generally  provide  80%  or  more 
of  the  correlation  energy,  require  no  more  computational  resources  than  do  SCF 
calculations.  Coupled-cluster  singles  and  doubles  (CCSD)  calculations  are  routinely 
performed  on  systems  up  to  the  size  of  benzene.  With  the  increasing  availability  of 
supercomputer  time,  it  is  easy  to  consider  such  calculations  on  molecules  up  to  twice 
that  large. 


1.2    The  Electron  Detachment  Problem 

Since  the  mid-60s,  "Molecular  Photoelectron  Spectroscopy  has  been  recognized  as 
an  especially  unambiguous  method  for  the  study  of  the  molecular  electronic  structures 
of  substances  in  the  vapor  state"  [3].  In  these  experiments,  hght  from  an  intense 
source  is  directed  into  a  molecular  gas  or  beam.  The  resulting  detached  electrons  are 
collected  and  their  kinetic  energies  recorded.  At  one  time,  it  was  common  to  record 
the  absorption  spectrum  as  well.  In  some  cases  fluorescence  studies  of  the  gas  are 
performed  as  well. 

Prom  the  beginning,  theoretical  methods  have  been  an  important  tool  in  inter- 
preting the  photoelectron  spectra.  Koopmans'  theorem  predictions  of  ionization  en- 
ergies made  it  possible  to  identify  principal  peaks,  and  later  CI  singles  and  doubles 
(CISD)  and  electron  propagator  approaches  were  used  with  some  success  to  iden- 
tify the  satellites.  In  the  late  seventies,  several  approaches  were  used  to  predict 
the  intensities  of  these  principal  peaks  and  their  satellites[4,  5,  6].  Nakatsuji  and 
coworkers[7,  8,  9]  introduced  the  symmetry-adapted-cluster  (SAC)  and  symmetry- 
adapted-cluster  CI  (SAC-CI)  approaches  to  calculating  the  ionization  energies  and 
intensities  of  molecules.  The  Pock-space  coupled-cluster  approach,  introduced  by 
Mukherjee  and  Lindgren[10,  11]  and  developed  by  others[12,  13,  14,  15]  has  recently 
provided  a  convenient  formalism  for  a  computationally  advantageous  approach  to 
the  direct  calculation  of  energies.  The  CC  EOM  approach,  closely  related  to  CC 
linear-response  theory  (LRT)[16],  is  based  upon  a  coupled-cluster  reference.  Por  the 


calculation  of  principal  ionization  energies,  this  theory  can  be  shown  to  be  formally 
equivalent  to  the  Fock-space  coupled-cluster  (FSCC)  method. 

In  this  dissertation,  derivations  of  the  CC  EOM  and  FS  EOM  methods  for  cal- 
culation of  ionization  energies  will  be  given,  and  the  formal  equivalence  of  these  two 
methods  will  be  demonstrated.  The  CC  linear-response-theory  will  be  derived  for 
the  ionization  problem,  and  expressions  for  calculation  of  transition  intensities  will 
be  developed  as  well. 

Finally,  the  CC  EOMIP  method  will  be  applied  to  several  small  molecules  and  to 
the  line£ir  carbon  cluster  anions,  C^  through  C^. 


CHAPTER  2 
ELECTRON  CORRELATION 

2A The  Schrodinger  Equation 

All  electronic  structure  methods  seek  solutions  or  approximate  solutions  to  the 
Born-Oppenheimer  approximation  to  the  Schrddinger  equation, 

H^    =    E^  (2.1) 

E    =    

where  H  is  expressed  in  Hartree  units  as 

^  =  -^E^?-E^or;/  +  E-i^^  (2-2) 

i  i,a  i>j 

which  is  the  sum  of  a  one-electron  and  a  two-electron  portion, 

i  ),Q 

y2-i/(two)     =     5:r^.\ 

where  i,j  are  indices  over  all  electrons  in  the  system  and  a  over  all  nuclei,  and  Za  is 
the  charge  of  nucleus  a.  Unfortunately,  the  kinetic  energy  term  of  the  Hamiltonian 

-^?'--J?((l:)'*(0*(i)T        - 


taken  with  the  two-electron  term,  gives  a  second-order  differential  equation  in  many 
variables,  for  which  there  is  no  general  analytic  solution.  It  is  possible  to  solve  certain 
approximations  to  the  Schrodinger  Equation. 


2.1.1     The  Independent  Electron  Approximation  and  Fock  space 


AH  the  methods  to  be  discussed  are  based  on  a  finite  set  of  spin  orbitals,  0p(C), 
one-electron  functions  of  the  combined  space  and  spin  coordinate  C  =  {r;<u^}.  For 
molecular  calculations,  the  spin  orbitals  are  referred  to  as  molecular  orbitals  and 
generated  using  an  independent  particle  method  -  in  which  the  V2  portion  of  the 
Hamiltonian  is  replaced  by  an  independent  particle  potential  Vp.  The  ultimate  form 
of  the  (p  orbitals  is  generally  arrived  at  by  a  self-consistent  field  method,  such  as  the 
Hartree-Fock  theory.  In  the  Hartree-Fock  approximation,  the  approximate  Hamilto- 
nian is  referred  to  as  the  Fock  operator. 

F    =    h  +  Vp  (2.4) 

H    =    F+{V2-Vf) 

The  basic  component  of  an  N-electron  wavefunction  is  a  single  Slater  determinant 
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^^^!        \r,=i,N 
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|01i/(Cl)^2i/(C2)  •  •  •  <f>NACN)\  , 


where  N  can  be  any  number  from  0  to  the  total  number  of  molecular  orbitals  available, 
and  the  index  u  distinguishes  among  the  f  ^j  possible  choices  of  A^  orbitals  from  a 
set  of  M.  $°  is  the  trivial  case,  with  no  electrons  chosen,  and  is  referred  to  as  the 
vacuum  state. 

The  Slater  determinant  is  an  antisymmetrized  product  of  A^  one-electron  func- 
tions, ensuring  that  fermion  statistics  are  obeyed  -  the  sign  of  the  wavefunction  will 
change  if  the  C  coordinates  of  any  two  electrons  are  exchanged.  Any  many-electron 
wavefunction  within  this  truncated  space  can  be  represented  as  a  linear  combination 
of  Slater  determinants. 

The  finite  space  spanned  by  the  Slater  determinants  is 

Nu 

which  is  generally  a  subset  of  the  space  spanned  by  the  eigenfunctions  of  the  Hamil- 
tonian.  The  Slater  determinants  span  the  eigenspace  of  the  truncated  Hamiltonian  - 
the  Hamiltonian  projected  into  the  space  of  the  Slater  determinants. 
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2.1.2    Second  Quantization 


The  second  quantized  notation  is  a  convenient  one  which  will  be  used  throughout 
this  work.  The  Hamiltonian  can  be  written  in  second-quantized  form  as 


pq  ^  pgrs 


(2.5) 


where 


V  =  /<^;(1)(-^V?  -J^ZaT^iy^dr^da,  (2.6) 

{pq\\rs)  =  j  j  ct>l{l)(p;{2){\  -  Pn)rl^(t>r{l)W^)dT,dT2da,da^.  (2.7) 

The  creation  (p^)  and  annihilation  (p)  operators  represent  the  creation  and  annihila- 
tion of  orbital  p.  The  fundamental  relations  between  these  operators  are 


ptgtj  ^     pt^t_^^tpt  =  Q 


(2.8) 

(2.9) 

(2.10) 


■'PQ- 


When  the  orbitals  are  orthonormal,  equation  2.10  becomes   p\q]     =  Sp 

Any  Slater  determinant  can  be  constructed  through  the  action  of  creation  opera- 
tors on  the  vacuum  determinant  ($°  =  |-)). 


g^..p^ 


^N\ 


i0p(Ci)-0,(C. 
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The  properties  of  these  operators  ensure  the  Pauli  principle  is  obeyed 


p^p^  =  0  —  p^p^    =    0 


pVI-)    =    0; 


a  given  spin  orbital  cannot  be  occupied  more  than  once  in  the  same  determinant. 

If  one  Slater  determinant  is  chosen  as  a  reference,  the  operators  can  be  separated 
into  occupied  or  hole  {i\  i)  and  virtual  or  particle  (a^  a)  creation  and  annihilation 
operators.  Any  Slater  determinant  can  created  from  a  given  reference  determinant 
|0)  by  first  annihilating  zero  or  several  occupied  orbitals  (creating  zero  or  several 
holes),  then  creating  0  or  several  virtual  orbitals  (creating  zero  or  several  particles), 
for  example  a'^z  |0)  represents  the  excitation  of  a  single  electron  from  occupied  orbital 
(pi  to  virtual  orbital  0a,  while  the  operator  string  a^ji  operating  on  the  same  reference 
excites  one  electron  while  annihilating  another,  reducing  the  number  of  electrons  in 
the  wavefunction  by  one.  A  thorough  discussion  of  the  nature  of  second-quantized 
operators  can  be  found  in  the  first  chapter  of  J0rgensen  and  Simons'  work[17]. 

In  order  to  simplify  further  discussions  and  diagrammatics,  it  is  useful  to  write 
operators  in  normal  ordered  (or  normal  product)  form[2,  18],  represented  by  placing 
"{  }"  around  the  product  of  operators.  The  normal  order  of  a  product  of  operators  is 
that  ordering  of  the  operators  whose  expectation  value  with  respect  to  the  reference, 
|0),  is  zero.  In  practice,  this  means  "moving"  all  particle  annihilation  operators 
(a)  and  hole  annihilation  operators  (z^)  to  the  right,  using  the  fundamental  anti 
commutation  relations  given  in  equations  2.8  through  2.10.  Interchanging  any  two 
single-particle  operators  within  a  string  of  operators  results  in  a  change  in  sign  of  the 
normal-ordered  product, 

{pqsr^}     =     —{pqr^s}  =  {pr^qs}  =  —{r^pqs} 


{aU}    =    aU. 

The  last  of  these,  aH,  composed  only  of  hole-  and  particle-creation  operators,  is 
the  only  type  of  normal-ordered  string  which  may  act  on  the  reference  |0)  without 
destroying  it.  With  normal  order  defined,  the  contraction  between  any  two  operators, 
designated  as  {pq)c  or  PQ  is  defined  as  the  difference  between  the  product  and  the 
normal  order  of  the  product 

P9  =  {pq}+PQ- 

Using  the  fundamental  anti  commutation  relations,  the  contractions  can  be  evaluated 
as 


pq  —  p^q^  =  0 

Jjf  =  ^  =  0 

i^j  =  5ij 

ab^  =  Sab- 


Another  common  normal  ordering  is  that  with  respect  to  the  vacuum  state  |— ). 
In  this  ordering,  all  orbitals  are  "virtual",  and  the  second  quantized  operator  in 
equation  2.5  is  already  in  normal  order.  This  is  not,  however,  the  case  with  the 
normal  ordering  used  here.  The  normal-ordered  Harailtonian,  Hj^  is  that  part  of  the 
Hamiltoniaji  whose  expectation  with  respect  to  the  reference  |0)  is  non-zero 

Hn    =    H-{0\H\0)  (2.11) 


In 


=    //  -  Escf 

=    In  +  Wn 

=    //at  (one)  +  if  AT  (two) 

PQ 
P9rs 


(2.12) 
(2.13) 


Writing  //^a-  in  terms  of  normal  ordered  operators  makes  this  partitioning  into  F  and 
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Figure  2.1:  The  Hf^  operator 

non-F  pieces  possible,  and  is  convenient  when  it  comes  to  developing  diagrammatics. 
A  diagrammatic  representation  of  these  Hj^  elements  is  given  in  figure  2.1. 


2.2     Coupled- Cluster  Theory 
2.2.1     The  Coupled- Cluster  Operators 

The  coupled-cluster  theory  defines  a  wavefunction  \CC)  which  is  an  eigenfunction 
of  the  full  Hamiltonian 

H\CC)    =    E\CC)  (2.14) 

(i/-(0|/f|0))|CC)    =    H^\CC)  =  Ecc\CC),  (2.15) 

where  Ecc  is  the  correlation  energy  obtained  from  the  application  of  coupled-cluster 
theory. 

The  CC  wavefunction  is  defined  in  terms  of  excitation  operators  T  as 

\CC)  =  e^  |0)  (2.16) 

The  operator  T  is  defined  as 

Ti  +  72  +  Ta  +  •  •  ■ 

ia 
ijub 

tjkabc 

etc.  The  inverse  of  this  exponential  is  simply  exp(— T).  Therefore 

e~^H\CC)    =    e~^Ee^\0) 
=    E\0) 
E    =    (0|  e-^i/e^  10) . 


The  correlation  energy,  Ecc  =  E  —  Ec,Qf,  is 


Ecc    =    (0|e"^ife^|0)-(0|//|0) 

=     (0|  e-^//;ve^  |0)  (2.17) 


2.2-2    Thg.^  Operator 


The  quajitity  exp(— r)iirexp(T),  referred  to  as  H,  shows  up  again  and  again  in 
coupled-cluster  based  methods,  so  it  is  worth  looking  at  it  in  detail.  H,  expanded  in 
terms  of  T  in  the   Baker-Campbell- Hausdorff  expansion,  is 

H    =    H  +  HT-TH  +  ^HT^-^T^H  -THT+---  (2.18) 

=    H+[H,T]  +  ^[[H,T],T]  +  ^[[[H,T],T],T]  +  ^[[[[H,T],T],T],T]. 

The  commutators  can  be  written  in  terms  of  their  connected  and  disconnected  parts, 
and  it  can  be  shown  that  the  disconnected  parts  cancel,  ah  owing  H  to  be  expressed 
as  a  sum  of  only  connected  pieces,  in  which  H  is  contracted  to  any  T  operators  in 
the  expression 

H    =    H  +  [HT),+  \{HT')^  +  ^X»T')c  +  l{HT')^  (2.19) 

=    [He^]^^  (2.20) 

The  connectedness  of  the  H  elements  is  what  establishes  the  size-extensivity  of  the 
coupled-cluster  method. 
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Figure  2.2:  .ff/v  elements  required  for  CC  EOMIP-SD  and  FSIP  equations 


2.2.3     The  Coupled-Cluster  Equations 


Similar  to  H  is  if  at  =  e  ^  H^e^  =  {HNe^)c:  which  appears  in  the  coupled-cluster 
energy  expression  in  equation  2.17.  Using  this  formulation,  we  have 


Ecc    -    (0|(if;ve^)c|0) 


(2.21) 


Since  H[^  |0)  =  £■  |0),  it  is  clear  that 


^|^;v|0)  =  0 


(2.22) 


for  any  state  ip  orthogonal  to  the  |0)  state.  This  fact  is  used  to  define  the  T  amplitudes. 
If 


a 

i 

ab 


=    aH  |0) , 
=    aUb^j  |0) , 


(2.23) 
(2.24) 


then  the  singles  and  doubles  equations  are  defined  by 


a 

i 

ab 


Hn\0)    =    0 
Hf,\0)    -    0 


(2.25) 
(2.26) 


Figure  2.2  contains  the  diagrammatic  expressions  for  the  H  terms  which  are  necessary 
for  the  methods  to  be  discussed  in  this  work.  For  example,  the  algebraic  expression 
for  the  term  "  (=  {i\f\a))  is 


(^|  /  \a)  +  Y^t^  {im\  \ae) 
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2.3     Configuration  Interaction 


2.3.1     The  CI  Operators 


In  configuration  interaction,  electron  correlation  is  introduced  by  constructing  a 
wavefunction  from  all  possible  Slater  determinants  with  the  same  number  of  occupied 
orbitals  as  the  reference.  This  means  each  additional  determinant  is  related  to  the 
reference  by  an  excitation  operation  -  aH  represents  a  single  excitation,  while  a^b^ji 
represents  a  double  excitation.  The  CI  wavefunction  is  represented  as 

\CI)  =  {1  +  C)  |0) ,  (2.27) 


where  C  is  defined  as 


Ci  +  C2  +  C3  + 


ia 
ijah 
ijkabc 


etc. 


2.3.2     The  CI  Equations 


The  CI  energy  is  calculated  from 


H\CI)  =  E\CI) 

E  =  (0ii/|C7)=E((0||0)  +  (0|C|0)) 

Eci  +  Escf  =  (0|(if;v  +  (0|i/|0))(l  +  C)|0) 

Eci  =  {0\Hn{1  +  C)\0) 


and  the  C  operator  satisfies 


a 

i 
ah 


H^il  +  C)\0)    =    EciCt 
Hn[1  +  C)\^)    =    EciCtj' 


etc. 
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(2.28) 
(2.29) 


(2.30) 


which  can  be  solved  via  a  matrix  eigenvalue  approach. 

In  the  CI  approach,  there  is  no  Baker- Campbell-Hausdorff  expansion  to  force  the 
connectedness  of  H  to  the  wavefunction  operator  C. 


CHAPTER  3 
THE  COUPLED-CLUSTER  EQUATION-OF-MOTION  METHOD 

3.1     The  EOM  Equations 

Beginning  with  the  exact  coupled-cluster  ground  state  \CC),  an  EOM  wavefunc- 
tion  is  created  using  a  linear  operator  fi,  the  full  set  of  which  spans  the  space  of 
excitations  from  the  SCF  reference,  and  therefore  the  space  of  eigenfunctions  of  the 
second-quantized  R^r  of  the  appropriate  number  of  electrons.  A  particular  EOM  state 
will  be  referred  to  by  the  index  /C. 

Hn\CC)    -    Ecc\CC)^Ecce'^\0)  (3.1) 

HN\fCcc)    =    {E''-Escf)\)Ccc) 

=    {E''-EscF)n''\CC)  =  {E''-Escf)n^e'^\0).  (3.2) 

When  Sl^  is  limited  to  contain  excitation-type  operators  only,  it  commutes  with 
T,  and  we  derive  the  expression 

i/^•  \}Ccc)  -  n'^H^-  \CC)    =    [E'^  -  E)  \Kcc)  =  [i/;v,  n""]  \CC) 
u'^e^Q.^  |0)     =     [Hn,  ^'^]  e^  |0)  =  [if/ve^,  ^^]  |0) 


e-^a;'^e^f]'^  |0)    =     [e'^  H^.e'^  ,^'']\Q) 
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uj'^n''\Q)    =     [Hr,-,n'^]\0)  (3.3) 

where  lj^  =  E*^  -  E  and  Hn  =  e-^H^-e'^  =  (i/yve^)c- 

Projection  of  equation  3.3  on  the  left  by  |h),  the  orthogonal  complement  to  the 
reference  function  |0).  and  by  |0)  itself  yield 

iu'^{h\n'^\0)    =    {h\[Hr,.n>^]\0)  (3.4) 

This  equation  defines  an  eigenvalue  problem  which  is  solved  to  yield  the  energies 
(eigenvectors)  to  (fi). 

The  second  term  of  the  commutator,  Q'^H^,  contains  only  uncontracted  terms. 
Dividing  the  commutator  into  contracted  and  uncontracted  terms  gives 

u'^n'^lO)    =    [Hn^'')^\0).  (3.5) 

The  commutator  in  the  EOM  equations  has  the  effect  of  removing  all  the  uncontracted 
terms  from  the  equations.  The  EOMIP  equations  become 

u;'^{h\n''\0)  =  {h\{Hr,n'')^\0),  (3.6) 

an  eigenvalue  problem  in  terms  of  the  transformed  Hamiltonian  Hs. 

The  expansion  of  the  H^^  terms  is  limited  to  a  maximum  Hi^  excitation  level  of 
(Ln  -  1),  where  Lq  is  the  highest  level  of  excitation  included  in  17'^,  and  again  by 
the  maximum  number  of  up-going  lines  in  Q,'^  available  for  contraction.  This  will  be 
2Lq  —  1  for  the  IP  problem  and  2Ln  for  excitation  energies. 


3.2     Size-Extensivitv  and  the  EQM  problem 

The  EOM  problem  is  similar,  from  a  computational  standpoint,  to  an  ordinary  CI 
problem  except  that  Hf^  is  not  symmetric  and  is  not  limited  to  one-  and  two-body 
terms.  FVom  a  theoretical  standpoint  however,  it  differs  considerably  from  the  CI 
problem,  in  that  all  the  terms  of  the  H^^'^  product  are  contracted. 

Within  this  framework,  a  perturbation  theory  which  is  connected  to  every  order, 
and  at  any  level  of  truncation  of  the  theory,  can  be  constructed  [10]  for  the  principal 
ionization  energy  and  electron  affinity  cases  (see  appendix  B).  Since  the  Hn  operator, 
being  based  on  the  CC  solution,  is  also  connected,  the  EOMIP  and  EOMEA  princi- 
pal energies  are  connected,  and  the  solutions  are  therefore  size-extensive.  Cl-based 
methods,  being  unconnected,  do  not  offer  this  advantage. 

3.3     The  ROM  Wavefimrtmn 

For  the  EOM  problem,  Q.^  is  a  linear  operator, 

^^    =    E^.S  (3.7) 

^' 

where 

^JO)    =    Im)-  (3.8) 

The  operator  r^  consists  of  a  normal-ordered  product  of  creation  and  annihilation 
operators  which,  acting  on  the  Hartree-Fock  reference  |0),  create  new  determinants. 
For  excited  wavefunctions,  the  set  {r^}  contains  the  identity  operator,  as  well  as 


excitation  operators  aU,  aHb^j,  ■  ■  ■.  For  the  electron-detached  case,  {r^}  contains 
i,  aHj  •  •  •,  operator  products  which  will  create  N-1  electron  determinants. 

Because  the  equation-of- motion  problem  is  treated  as  an  eigen  problem  of  the  H^ 
operator,  the  solution  wavefunctions, 

\ic)  =  n'^io) 

are  eigenfunctions  of  H^  within  the  truncated  space  defined  by  the  highest  level  of  ex- 
citation in  {r^}.  But  in  general,  we  are  interested  in  wavefunctions  of  the  Hamiltonian 
H,  and  so  must  go  back  to  the  expression 


\)Ccc)    =    e^f^'^lO) 


(3.9) 


an  approximation  to  an  eigenfunction  of  H . 

If  we  define  a  matrix  C  with  vectors  C^  corresponding  to  the  C^  which  appear  in 
the  definition  of  fi'^,  then  Q  (  the  entire  set  of  Q'^  solution  operators)  can  be  written 


as 


Q    ^    tC 


(  C\    Cl    ■■■    Cf  \ 


C\    Cl 


=     (  ri    T2    •  •  •    r^  ) 

\Cl    Cl    •••    C^  ) 
the  solution  space  for  a  given  EOM  problem  can  be  represented  as 

(0|r^/vrC|0)    =    uC 
C-'{^i\Hr,\^i)C    =    u 


(3.10) 


D^{^i\H!,\^i)C  =  u. 


(3.11) 


The  matrix  D^  =  C"\  defines  a  new  operator  Y  =  tD.  D  and  Y  can  also  be  used 
to  represent  the  solution  space  of  the  EOM  equations: 

DM/i|^;v|/i)    =    Du.  (3.12) 

It  is  clear  that  these  D  vectors,  the  left-hand  eigenvectors  of  the  non  Hermitian  H^, 
yield  left-hand  eigenfunctions  of  H;^, 

D^(0|r^NT|0)    =    D^w  (3.13) 

{Q\Y'^^Hn    =    (Oir'^V,  (3.14) 

and  therefore  the  projection  of  the  full-space  left-hand  eigenfunction  into  the  reduced 
space, 

{0\Y'^^e-'^Hi,    ^    {0\Y''^e-'^E'',  (3.15) 

just  as  the  right-hand  eigenvectors,  C,  yield  the  projection  of  the  full-space  right-hand 
eigenfunction  into  the  reduced  space, 

Hj^Q'^e^lO)    ^    E^Q'^e'^lO).  (3.16) 

It  should  be  noted  that,  since  the  untransformed  H  is  a.  Hermitian  operator,  its 
right-  and  left-  hand  eigenfunctions  are  of  course  merely  the  transpose  of  one  other: 

H\R,)    =    E,\R,) 
{U\H    =    {U\E, 

Up  to  this  point,  there  has  not  been  much  attention  paid  to  the  T  and  ^^/Y'^ 
operators  actually  used  in  practical  calculations.   Although  it  is  easy  to  show  that, 


in  the  full  expansion  of  the  operators,  the  left-hand  side  ground  state  coupled  cluster 
wavefunction  and  the  transpose  of  the  right-hand  side  wavefunction  are  identical  up 
to  a  normalization  constant, 

(0|(1  +  A)e-^    =    (0|e^^e^|0)-^(0|e^\ 

qualitative  differences  between  the  two  wavefunctions  are  apparent  in  the  usual  CCSD 
approximation  (John  Stanton,  private  communication),  because  of  the  truncation  in 
the  operator  space. 

3.4     The  Ionization  Energy  Equations 

In  an  ionization  potential  problem,  the  reference  state  |0)  is  an  N-particle  function, 
whereas  the  solution  Jl'^  |0)  is  an  (A'^  —  1)  particle  function.  The  operator  Q'^  must 
therefore  contain  terms  which  annihilate  on  electron  in  addition  to  any  excitation 
they  might  also  perform: 

Q'^    =    C'^  +  C^  +  ...  (3.17) 

=  E^^  +  ^Ecf>V  +  ---  (3-18) 

1 


2 

I]  a 


;v  (3.19) 


\)c)  =  n''\o)  =  ^cf  J  +  T^E^; 

The  IP  coefficients  c^,  cf'j  are  indexed  by  /C  for  convenience.  The  similarity  of  this 
notation  to  that  for  the  coupled-cluster  amplitudes  tf  and  tf^  implies  that  an  electron 
is  being  excited  out  of  orbital  i  and  into  fC,  but  this  /C  does  not  represent  one  of  the 
reference  orbitals.  It  may  be  taken  to  indicate  that  the  electron  has  been  excited 
into  an  "unbound  orbital,"  but  this  interpretation  is  not  important  to  the  solution  of 


the  EOMIP  equations.  It  is  important  to  remember  that,  when  constructing  EOMIP 


equation  terms,  no  contractions  may  be  made  to  the  "free-electron"  index. 
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Figure  3.1:  The  EOMIP  equation  diagrams. 


The  EOMIP  singles  and  doubles  (EOMIP-SD)  equations,  for  which  the  terms  are 
shown  diagrammatically  in  figure  3.1,  can  be  written  as 


)C„K 


(siNGLES)a;    c- 


-E("l/K)Cn+E(«l/l«)^' 
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(3.20) 
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-W(i:i^^\M^^i]^nL  (3-21) 

nmb  \   e  I 

Where  the  notations  /  and  (•  •  •  ||  •  •  •)  indicated,  respectively,  one-  and  two-body  ti 
elements.  The  final  term  of  the  doubles  equation  involves  the  3-electron  Ht^  term 
lnraa\H\ij^  (figure  2.2).  These  equations  form  an  eigenvalue  problem  which  may 
be  solved  by  any  number  of  diagonalization  procedures. 

3.5    The  EOMIP  Wavefunction 

The  EOMIP-SD  wavefunction  expressed  in  equation  3.19  is  an  eigenfunction  of 
the  Ht^  operator  projected  into  the  reduced  EOMIP-SD  space.  The  full  wavefunction, 
e^fi'^  |0),  is  an  approximation  to  an  eigenfunction  of  if,  and  can  be  written 


i..c)  =  Ecf|J.(ii:c-.^(.)Ecf<")|.;) 


^  \  ijkab  ijkab  / 

+  l-(p{i/j/kl)P{a/bc)    Y:   cftpl^,  +  \p{i3kl)Piabc)    ^   cftplt^ 


ijk^ 
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CHAPTER  4 
THE  FOCK-SPACE  COUPLED-CLUSTER  METHOD 

4.1      Fock-Space  Theory 

The  Fock-space  coupled-cluster  method  (FSCC)  uses  an  exponential  operator  to 
describe  electron  attached,  electron-detached,  and  excited  states  relative  to  one  single- 
determinant  reference.  Like  EOMIP,  it  gives  direct  differences  between  the  correlated 
reference  energy  and  the  energy  of  the  calculated  state.  The  FSCC  approaches  to 
the  calculation  of  single  electron  ionizations  and  attachments  are  formally  equivalent 
to  the  EOMIP  approaches  to  the  same,  as  will  be  demonstrated  for  the  IP  case  in 
section  4.2.  The  Fock  space  coupled-cluster  method  is  based  upon  a  valence  universal 
operator  fi  that,  when  acting  upon  the  appropriate  model  space,  yields  the  desired 
wavefunction(s).  If  we  begin  with  a  single  determinant  function,  such  as  a  Hartree- 
Fock  wavefunction,  a  model  space  can  be  defined  by  annihilating  n  electrons  from 
occupied  orbitals  (also  referred  to  as  creating  holes)  and  creating  m  electrons  in 
virtual  (unoccupied)  orbitals  (also  referred  to  as  creating  particles).  This  is  called  a 
rank  {m,n)  model  space,  one  with  particle-hole  rank  {m,n).  The  rank  {m,n)  part  of 
the  problem  is  also  referred  to  as  the  {m,n)  sector.  For  exEimple,  the  usual  coupled- 
cluster  problem  has  particle-hole  rank  (0,0),  while  the  single  electron  detachment 
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problem  is  of  rank  (0, 1).  An  active  space  of  occupied  and  virtual  orbitals  is  chosen  in 
which  electrons  can  be  created  and  annihilated  to  produce  model  space  determinants. 

The  operator  f2  is  defined  by  (e-^|,  where  the  braces,  {},  denote  normal  ordering. 
Normal  ordering  of  a  product  of  operators  prevents  the  operators  from  contracting 
with,  or  operating  on,  each  other.  For  this  reason,  model  spaces  of  different  hole- 
particle  ranks  are  decoupled,  and  the  solution  of  the  equations  is  hierarchical,  start- 
ing with  the  usual  coupled-cluster  equations,  and  building  solutions  to  sectors  with 
successively  higher  rank  (the  sum  of  m  and  n)  upon  the  solutions  of  lower  sectors. 

The  valence  universal  operator  ft  can  be  written 

ntot={e'''}  (4.1) 

where  X  —  5^°'°'  -I-  5^°-^'  -h  5-^'°'  +  ■  •  ■.  Because  of  the  properties  of  normal  ordered 
products,  QiqI  can  be  written  as  a  product  of  terms 

n^^t  =  {n(°'°U  fi(°'iu . . .}  (4.2) 

where  f^t'"'")    =    {exp(5("^'"))}  =  1  +  5("^'"' +  ^{(5("^'"^)'}  +  •  ■  • . 

g{m,n)  jg  always  chosen  so  that,  when  acting  on  the  (m,  n)  model  space  F('"'"\  it  does 
not  cause  scattering  between  model  determinants.  That  is,  p('".")5^"''")p('"'")  =  0, 
or  s^"''''^P^"''''^  =  Q[m,n) g{m,n) p{m,n)  ^  ^^^^.^  Q{m,n)  jg  ^^g  complement  to  p("''").  The 

(0, 0)  sector,  as  was  noted  above,  corresponds  to  the  usual  single- reference  coupled- 
cluster  problem.  The  operator  exp(5^°'°')  is  identical  to  the  usual  e^  operator  encoun- 
tered in  the  coupled-cluster  equations.  Because  T-to-T  contractions  are  impossible, 
imposing  normal  order  does  not  change  the  operator. 


The  operator  ^2^^^  has  an  inverse, 

f^^^^-i  =  {f^(o.o)-%Q(o.i)-%...}.  (4.3) 

In  the  (0, 0)  case  this  is  Q^^''^'  =  e"^  but  in  the  higher  sectors,  the  form  becomes  more 
comphcated  because  of  the  normal-ordering  constraint.  These  terms  can  easily  be 
derived  using  a  Taylor  series  expansion  of  the  inverse  of  Q^"^'''\  or  by  other  approaches, 
but  are  not  important  to  the  following  discussion. 

To  begin,  an  iV-€lectron  single  determinant  reference  |0)  is  chosen,  which  may  be 
an  SCF  wavefunction.  The  reference  wavefunction  defines  a  complement,  Q,  which 
contains  all  iV-electron  functions  orthogonal  to  the  reference.  The  single-reference 
coupled-cluster  problem  is  solved,  yielding  H^,  which  is  an  effective  Hamiltoni.a'' 
for  the  (0, 0)  sector.  This  is  a  canonical  (eigenvalue-preserving)  transformation  of  t '  e 
Hamiltonian,  chosen  so  that  it  is  the  eigen-operator  of  the  reference  function,  yielc.ng 
the  desired  coupled-cluster  energy  {Ecc)  as  its  eigenvalue. 

The  coupled-cluster  equations  can  now  be  expressed  in  a  compact  form. 

(Q|^^|0)  =  0,  (4.4) 

and  the  equations  are  solved  for  the  T  amplitudes. 

In  further  discussions,  the  following  quantity  will  be  useful: 

[o<;+it  J 

While  ri^'7j"'  contains  all  the  contributions  from  all  the  f2's  of  different  sectors,  the 
old 

projection  of  equations  by  P^""'")  and  Q^""'"'  limits  the  contributions  to  the  energy 
and  wavefunction  expressions  to  only  certain  lower-ranked  sectors    The  (0, 0)  sector 


contributes  to  all  higher- ranked  sectors,  and  is  the  only  one  contributing  to  the  (0, 1) 
and  (1,0)  sectors.  The  (1, 1)  sector  has  contributions  from  all  lower-ranked  sectors, 
while  the  (0,2)  has  contributions  only  from  (0, 1)  and  (0,0). 

For  a  given  sector,  an  effective  Hamiltonian  is  also  defined,  W^\  defined 

which  will  prove  useful  in  the  following  discussion.  Now,  with  the  model  space  and 
the  f2  operator  defined,  we  can  begin  to  define  the  {m,n)  problem. 

While  H^°^'P°'°  =  EccP^^'^^  works  fine  for  the  one-dimensional  problem,  it  is 
clear  that  something  more  is  needed  to  deal  with  the  multi-determinant  model  space. 
It  is  necessary  to  introduce  C^'"'"\  or  C  for  short.  The  matrix  C  causes  a  linear 
transformation  of  the  model-space  functions.  If  F^^"'"'  includes  a  set  of  functions  \pi), 
then  a  linear  transformation  of  these  functions  results  in  another  set 

or 

I 

The  X:*^  wavefunction  of  the  sector  then  is  defined  to  be  fiio^C{;'"'"'P^'"'"'  and  should 
satisfy 

If  the  full  C  matrix  is  used,  the  result  is 

HVLtotP^"''''^C^"'''"^  =  Q^oi^''"'"'C^^'"'"^W^'"'"\  (4-8) 

where  W^"*'"^  is  a  diagonal  matrix  containing  the  eigenvalues. 


This  equation,  like  the  analogous  (0,  0)  equation,  can  be  put  into  an  "-W^g/f "  form, 

where 

Hgff  has  the  property 

Q{m,n)  TT{m,n) p[m,n)^(m,n)      __      Q{m,n) p{rn.,n) Q(m,n)y^{Tn,n) 

=    0  (4.10) 

This  fundamental  result,  unfortunately,  is  not  nearly  as  useful  as  the  similar  equation 
from  the  (0, 0)  sector  is  in  computing  the  solution  to  the  coupled-cluster  equations. 
The  complicated  form  of  57"^  makes  for  a  large  number  of  algebraically  messy  and 
computationally  difficult  terms  in  the  expression  Q^"^'"^H^feCP^"^'"-K  Upon  conver- 
gence of  Q, 

Tj(m,n) p[m,n)  /j{m,n)  _  / p{m,n)    ,    /^(m,n)-i  rr("i,rj)p(m,n)y-i(m,n)  _  p{m,n)  tj        p{m,n) ^(m,n) 

(4.11) 
Since  the  matrix  C  only  mixes  the  P^""'"'  functions  among  themselves,  and  cannot 
promote  pf'"'"'  functions  to  the  Q^"*-"'  space,  it  is  equally  valid  to  say 

^{m,n) p(m,n)  _  p{m,n) ^(m,n) p{m,n)  /"^  ■^2'] 

(^ff  eff  ■  ^  ■     ^ 

Combining  this  with  the  identity  Hj'^Q  —  QH^rr,  we  arrive  at  an  algebraically  and 
computationally  more  tractable  form,  Q('"'"'(.H',vfi  -  f]F^'"'")i/g«-)F('"'"'.  Once  Q 
has  been  determined  from  the  solution  of  this  equation,  p("»-n)//'|^'^")p("»'")  can  be 
diagonalized  to  yield  the  C^  and  the  energies  uj^^    . 


Hff    =    -{Af\j)  +  ESin{m\f\j) 

m 

-  E  S%  {m\  f  \e)  +  Y:  S^mnWje).  (4.13) 

me  mne 

0    =    {m\J)-ESln{m\f\j)  +  ESrm{rri\f\e) 

m  me 

-  ^  E  S'^mn\\je)  +  E-^f/Zf/  (4.14) 

■^  mne  m 

0    =    {ia\\jk)-ESin{'^\\jk)-PwESTm{m\f\k) 

m  m 

+  E  S%  {a\  I  \e)  -  P(,fc,  E  S%{ma\\ej)  +  W  S-^{mr,\\3k) 

e  me  ^  mn 

-\Y.S'L  tT,{mn\\bc)  -  Y:  Sf,^Hf^  (4.15) 

mnbc  "* 

An  important  feature  of  these  equations  is  that  the  ionization  energies  are  invariant 
to  the  choice  of  active  space.  This  is  because  the  S\  amphtudes  and  the  C  coefficients 
are  serving  essentially  the  same  purpose;  they  mix  the  model  space  functions.  If  all 
occupied  orbitals  axe  made  active,  there  will  not  be  any  Si  amplitudes,  and  an  n  x  n 
matrix  of  C  coefficients.  If  the  number  of  active  orbitals  is  decreased,  there  will  be  a 
decrease  in  the  number  of  C  coefficients  and  a  corresponding  increase  in  the  number 
of  5i  amplitudes.  In  fact,  it  is  easy  to  show  that  there  is  a  simple  relationship  between 
C  and  S  for  the  two  cases.  The  proof  of  the  invariance  of  the  energy  to  the  choice  of 
active  space  can  be  found  in  appendix  A 


4.2    Comparison  of  EOM  to  FSCC 

When  FSCC-IP  ajid  EOMIP  calculations  are  performed  using  the  same  active 
space — the  complete  active  space,  consisting  of  ionizations  from  all  occupied  orbitals — 
certain  correspondences  can  be  shown. 


4.2.1     Definitions 


r2^°'°'  is  the  usual  CC  wave  reaction  operator  e^. 

p(o>i)^  written  P  for  convenience,  is  a  projector  over  the  model  space,  which 
contains  all  (A^  —  1)  -  electron  determinants  obtained  by  annihilating  a  single 
occupied  electron. 

Q^^'^\  or  Q,  is  the  complement  to  P  (1  -  P)  in  the  {N  —  1)  -  electron  space. 

^(0,1)  =  {exp(5)}  =  1  +  5  +  \{S'^}  +  •  •  •  is  the  solution  to  the  FSCC  equation 
for  the  (0,1)  (IP)  sector. 

fi(o-i'"'  satisfies  fi(o-i)"'f^Co.i)p  =  Q(o-i)f2(o.i)-'F  =  P. 


4.2.2    FSCC  Method 


The  Fock-space  method  defines  a  model  space  consisting  of  several  (TV  —  1)  elec- 
tron determinants  and  a  projector  P  =  p(°'^)  over  this  space.  An  electron-detached 
wavefunction  is  defined  as 

i 

P  can  be  represented  as  a  sum  over  these  |/C)  states 

p  = -£1)0  mfcrU^- 

K 

The  Jl(°'^'  and  C'^  operators  are  defined  so  that 

^^('^•i)  |X:)  =  E^^Q.^^''^  \1C) ,  (4.16) 

which  leads  to  a  useful  definition  for  H^tr , 

H^ffP    =    f2(°'^)"'^fi(°'i'P;  (4.17) 

=    cu>^\)C),  (4.18) 

which  has  the  property 

QH^ff  \)C)  =  QE'^  \K.)  =  0  (4.19) 

H^ff  =  f](o>i)"'^<^(o-i>F  =  Ff](0'i'~'^J^(°'i)p.  (4.20) 

H^a  acting  on  any  representation  of  the  model  space  can  be  diagonahzed  to  yield 
the  entire  spectrum  of  ionization  energies  u:^ 

C~'H^ffC\0)^cu\0)  (4.21) 


This  operator  eigenvalue  equation  is  easily  replaced  by  the  matrix  eigenvalue  equation, 


C-%^C  =  u, 


4.2.3     EQM  Method 


The  diagonalization  of  H  yields  coefficients  C  and  eigenvalues  uj.  The  space  of 
eigenvalues  and  eigenvectors,  and  hence  the  H  matrix  itself,  can  be  partitioned  into 
principal  and  shake-up  portions.  The  principal  IPs  correspond  to  the  solutions  of  the 
Fock  Space  method  and  can  be  labeled  as  the  P  portion  of  the  solution  space. 

(  Hpp    HpQ  \  I  Cpp    CpQ  \     _     f  Cpp    CpQ  \  (  Up     0  A 

\   Hqp      Hqq    J    \   Cqp      Cqq   j  \   Cqp      Cqq   )    \     0        UJq   J 

or 

(HppCpp  +  HpqCqp    HppCpQ  +  HpqCqq  \     _     I  Cppuip    CpqUq  \ 
HqpCpp  +  HqqCqp      HqpCpQ  -h  HqqCqq    J  \   CgpUp      CqqUq   J 

The  P-P  portion  of  this  equation  can  be  rewritten  as 

Cpp  I  Hpp  +  HpqCqpCpp]  Cpp  ~  Up. 

Since  up  has  been  chosen  to  match  the  solutions  to  the  FS  equations,  there  is 
quite  naturally  an  identification: 

C    HgjrC  =  Up  =  Cpp  (Hpp  -\-  HpqCqpCppJ  Cpp, 

which  will  lead  to 

Hg^  =  (Hpp  -H  HpQCqpCppj 

as  long  as  the  model  functions  chosen  for  both  methods  are  identical. 


36 
4.2.4     Comparison 


The  eigenvalues  and  eigenvectors  obtained  from  the  diagonalization  of  the  effective 
Hamiltonian  i^gff  are  the  {PP)  blocks  of  the  respective  matrices  obtained  from  the 
diagonalization  of  the  full  matrix  H 

"eff    =    "^r+{fS),,,  (4.22) 

=  Hpp-¥{HpqS)^  (4.23) 

=  Hpp  +  HpqCqpCpp  (4.24) 

Sqp    =    CqpCpp  (4.25) 

Cqp    =  SqpCpp  (4-26) 

This  means  that  the  eigenvectors  of  principal  IP  states  obtained  from  an  EOM 
calculation  can  be  directly  compared  to  those  from  the  FSCC  calculation  of  the  same 
states.  The  Cpp  vectors  obtained  in  both  examples  are  identical. 


CHAPTER  5 
THE  ELECTRON  PROPAGATOR 

The  electron  propagator [19],  also  called  the  one  electron  Green's  function,  was 
one  of  the  first  improvements  made  on  Koopmans'  theorem  for  the  calculation  of  IPs 
and  EAs  [20],  and  is  still  in  use.  Propagator  methods  are  direct  methods,  >ielding 
energy  differences  rather  than  state  energies.  Unlike  the  methods  described  so  far, 
which  are  based  upon  correlated  wavefunctions  for  the  reference  and  final  states,  most 
Green's  function  methods  provide  transition  properties  rather  than  wavefunctions  for 
the  finaJ  states.  And  it  is  feasible  to  begin  a  correlated  calculation  of  ionization  en- 
ergies with  an  uncorrelated,  single  determinant  reference,  if  a  suitably  large  manifold 
of  annihilation/excitation  operators,  such  as  those  used  in  the  the  construction  of 
the  EOMIP  wavefunction,  is  chosen.  However,  the  symmetric  form  of  the  electron 
propagator  requires  that  the  manifold  include,  in  addition  to  IP-type  operators,  i 
and  aUj,  the  EA-type  operators  a^  and  a^bH  be  included  as  well.  Energy  differences 
will  correspond  to  the  difference  between  a  correlated  reference  and  a  correlated  final 
state,  neither  of  which  wavefunction  is,  in  general,  directly  constructed. 

Despite  these  apparent  differences,  there  are  similarities  between  the  results  of 
electron  propagator  and  CI  calculations.  More  recent  formal  developments  have 
bridged  the  theoretical  gap  between  Green's  function  and  coupled-cluster  methods. 
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and  propagator  researchers  have  independently  arrived  at  equations  which  are  iden- 
tical to  those  derived  from  coupled-cluster  EOM  theory. 

The  general  form  of  Green's  function  in  the  energy  representation  is 

//.   ou         r      i^±{0\A\k){k\B\0)      ^  {0\B\k){k\A\0)\ 

where  the  states  \k)  are  eigenfunctions  of  the  full  Hamiltonian,  H.    The  electron 
propagator  is  a  matrix  G  of  elements 


'P9 


f  ^   {ON\p^\kN-^i){kN-i\q\0)       ^^   {0\q\kN+i){kN^i\pUO)]      .52) 


where  the  l/c^-i)  and  I^at+i)  functions  are,  respectively,  electron-detached  and  electron- 
attached  functions.  The  poles  of  this  expression  occur  at  £'0  -  -Efcjv-i  +  E  =  0  and 
Ekj^  ^  -  Eq  +  E  =  0;  the  residue  at  a  particular  pole  yields  the  transition  moment 
{On\p^  \kN-i)  {kN-i\q  |0)  or  (0|  q  l/c^+i)  {kN+i\p^  |0).  This  form  of  the  Green's  func- 
tion is  not  directly  useful  for  obtaining  energies  or  wavefunctions,  so  another  approach 
is  used  which  calculates  the  energy  differences  and  transition  moments  directly,  with- 
out first  determining  [0),  \k),  or  any  of  the  state  energies. 

5.1     Superoperator  Expression  of  the  Electron  Propagator 

Propagator  equations  are  most  compactly  written  in  terms  of  superoperator  no- 
tation. For  any  operator  C,  the  superoperator  C  acts  on  other  operators,  and  a  new 
"binary  superoperator  product"  is  defined  between  any  two  operators: 

CD    =    [C,  D]__=CD-DC  (5.3) 


ID    =    D 


{C\D)  =  {o\c^d\o)t{o\dcUo), 


(5.4) 
(5.5) 


where  the  sign  in  equation  5.5  depends  on  the  number  of  creation  and  annihilation 
operators  in  C  and  D.  In  this  notation,  the  Green's  function  can  be  written  as 


Inserting  the  spectral  resolution  of  the  identity, 

1    =    JF')(F'|F')"'(Ft| 

=  EiF')(F'i''*):;(^'i. 

kl 

into  the  expression  for  the  propagator  gives 


(5.6) 


(5.7) 
(5.8) 


{{A;B))^    =    {B^Ei  +  Hy'\A) 

=    (B+|Ft)  (Ft|Ft)''  (Ft|  [Ei  +  h)''  \f')  {F^IF^Y'  {f^\A) 

=    {B^\F^){F^Ei  +  H)\F^y\F^A),  (5.9) 

which  are  the  working  equations  from  which  propagator  approximations  are  derived. 


The  expression 


[ei  +  h)  ' 


is  called  the  superoperator  resolvent.  The  poles  and  residues  of  the  resolvent  give  the 
energies  and  transition  moments. 

For  the  one-electron  propagator,  propagator  matrix  elements  become 

G,,{E)    -     {{p^;q))^  =  ip\{Ei  +  Hy'\q^) 
=    (p|F)(f|(M  +  ^)f)"'(F|9). 


(5.10) 


The  manifold  F  can  be  separated  into  a,  which  contains  all  the  one-electron  operators, 
p  and  q,  and  its  orthogonal  complement  f ,  which  contains  all  higher  odd-electron  op- 
erators. Because  the  f  portion  of  the  manifold  is  orthogonal  to  the  p  and  q  operators, 
only  the  a  -  a  projected  portion  of  (f\  (eI  +  Hj  |Fj  contributes.  This  portion  can 
be  written  as 

G-\E)    -    (a|(Ei  +  ^)|a)- (a|(£i  +  ^)|f)(f|(M-h^)|f)"'(f|(M  +  ^)|a) 
=    £1- (a|^|a)  -  (a|i/|f)(f|(£:i  +  i/)|f)"^(f|if|a)  (5.11) 

5.2    The  Dyson  Equation 

The  Dyson  equation  is  an  expression  of  the  difference  between  the  inverses  of  the 
propagator  G~^{E)  and  the  zeroth  order  propagator  Go^{E),  the  Green's  function  in 
terms  of  the  zeroth  order  Hamiltonian,  Hq. 

G-\E)    =    G-,\E)-nE)  (5.12) 

where 

G-,\E),,    =    {E-e,)5,,  (5.13) 

E(^)    =    (a|Ha)^^^^  +  (al^f)(f|(£i  +  Ar)f)"^(f|^a)        (5.14) 

(p|^l«?)corr    "    (pI^I?)  -  h^Pi- 
In  this  formulation.  S(£'),  the  self-energy,  is  the  correction  to  the  Koopmans'  theorem 
approximation  to  the  propagator  given  by  Gq[E).   The  self  energy  is  itself  divided 
into  two  parts,  the  energy  dependent  part  T!{E)  and  the  constant  part  i;(oo) 

E(E)    =    S(oo)  +  E'(£:)  (5.15) 


E(cx))    =    (al^a) 
^     ^  V  '      /corr 

E'{E)    =    (a\m){{\[Ei  +  H)fy\f\Ha). 

(5.16) 

A  great  many  approaches  to  the  electron  propagator  are  concerned  with  successive 
orders  of  approximation  to  the  energy-dependent  part  of  the  self  energy [21], [22],  in 
combination  with  appropriate  choices  for  the  reference  state  |0)[23],  and  the  most 
important  portions  of  the  manifold  F. 

5.3    Non-Dyson  Equation  Approaches 

The  development,  in  the  world  of  CI  theory,  of  the  Davidson  procedure  for  con- 
vergence of  selected  roots  of  a  large  matrix[24]  spurred  Baker  and  Pickup  to  develop  a 
method  for  extracting  the  roots  of  the  one  electron  propagator  matrix  rather  than  us- 
ing self-energy  expressions[25].  More  recently,  Ortiz  has  experimented  with  the  use  of 
a  highly-correlated  reference  based  on  the  coupled  cluster  wavefunction  [26],  a  process 
referred  to  as  renormalization  of  the  ground  state,  in  a  procedure  which  maintains 
the  symmetry  of  the  electron  propagator  matrix.  Snijders  and  Nooijen  [27],  using  an 
RSPT  analysis  of  the  diagrams  contained  in  the  one  electron  propagator,  showed  that 
they  can  be  represented  using  the  coupled-cluster  wavefunction  as  the  reference,  with 
two  major  consequences.  The  svinmetry  of  the  resulting  matrices  is  broken,  and  the 
IP-type  manifold  is  decoupled  from  the  EA-type  manifold,  decreasing  the  dimension 


of  the  manifold  space  to  less  than  half  for  the  IP  case  when  typical  modern  basis  sets 
are  used.  The  resulting  equations  are  identical  to  the  CCEOM  equations. 

One  other  Green's  function-based  approach  which  gives  results  identical  to  the 
EOM  approach  is  the  coupled-cluster  response  function  approach  developed  by  Koch 
and  J0rgensen[16].  In  this  approach,  a  time  dependent  formulation  is  used  to  derive 
expressions  for  the  linear  and  quadratic  response  functions  of  the  coupled-cluster 
equations.  The  coupled-cluster  Jacobian  (the  first  derivative  of  the  coupled-cluster 
equations  with  respect  to  the  cluster  amplitudes)  used  in  their  work  is  identical  to 
the  coupled-cluster  EOM  matrix.  Transition  matrix  elements  between  the  ground 
and  excited  states  and  between  two  excited  states  are  derived  from  residuals  of  the 
quadratic  response  function. 


CHAPTER  6 
OTHER  METHODS  FOR  IONIZATION  CALCULATIONS 

The  many  approaches  to  calculating  ionization  energies  can  be  broken  down  into 
two  main  categories  -  state-by-state  approaches  which  involve  a  sepeirate  SCF  refer- 
ence and  correlated  calculation  for  each  state  of  interest,  and  direct  methods  which 
use  a  single  reference  (SCF  or  correlated)  as  the  starting  point  for  calculation  of  the 
energy  difference  between  the  reference  state  and  several  states  of  interest.  State-by- 
state  approaches  can  use  any  SCF  or  correlated  (CI,  perturbation,  or  coupled-cluster 
theory)  method  in  the  calculation  of  differences.  Among  the  direct  approaches  in- 
clude two  hole-one  particle  CI[28,  6],  the  symmetry  adapted  cluster-CI  approach  of 
Nakatsuji[29],  the  Fock-Space  Coupled  Cluster  approach[lO,  13,  30,  14],  and  one  elec- 
tron propagator  appro  aches  [31,  32,  23,  33],  as  well  as  the  coupled-cluster  equation-of- 
motion  approach,  which  has  been  derived  independently  several  times,  and  is  known 
variously  as  the  coupled  cluster  linear  response  function  [16]  and  the  coupled-cluster 
Green's  function[27]  methods. 
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6.1     State  Methods 

The  conceptually  simplest  way  to  consider  calculating  energy  differences  is  to 
calculate  the  individual  energies  of  the  states  of  interest  and  take  the  differences. 
The  problems  inherent  in  comparing  SCF  energies  of  systems  with  different  spins 
(as  in  the  common  case  of  the  closed-shell  singlet  neutral  reference  and  the  doublet 
ionized  state)  are  well  known  and  tend  to  make  the  SCF  energy  differences  not  even 
qualitatively  useful.  Configuration  interaction,  perturbation  theory,  and  coupled- 
cluster  methods  may  all  be  applied,  and  for  small  molecules  which  are  well-described 
by  a  single  SCF  determinant,  any  of  these  methods  might  give  a  good  result  for 
valence  ionizations.  There  are  a  few  difficulties  which  can  make  the  problem  less 
straightforward,  however.  Modern  SCF  programs  usually  can  be  induced  to  converge 
on  valence  ionized  states,  and  the  fact  that  core  molecular  orbitals  are  generally  well 
described  as  a  linear  combination  of  core  atomic  orbitals  and  do  not  change  much 
when  the  occupation  changes  means  that  core  ionized  states  can  be  converged  as  well. 
But  when  an  inner-shell  ionization  is  desired,  there  is  no  good  way  to  characterize  the 
excited  SCF  state.  SCF  programs  will  quickly  collapse  back  to  the  valence-ionized 
state.  One  remedy  for  this  problem  is  to  use  a  quasi-restricted  Hartree-Fock  (QRHF) 
reference,  in  which  an  SCF  calculation  is  performed  on  suitable  reference  such  as 
the  closed- shell  neutral,  following  which  the  orbitals  are  re-occupied  and  the  Fock 
matrix  recalculated.  This  non-Hartree-Fock  reference  is  then  used  as  the  basis  for  the 
correlated  (typically  coupled-cluster)  calculation [34]. 
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Because  each  state  is  handled  as  a  separate  case,  there  is  considerable  freedom 
in  the  choice  of  the  reference.  For  approximate  correlation  methods,  poorly  chosen 
orbitals  result  in  a  substantial  shift  in  the  energy  of  the  state  considered.  And  if  the 
orbitals  are  chosen  well  for  one  state  but  poorly  for  another,  the  lack  of  balance  means 
the  accuracy  of  the  diiference  will  suffer.  UHF  references  lack  spin  symmetry.  ROHF 
and  QRHF  references  for  high-spin  states  will  always  be  pure  spin  states,  but  the 
orbitals  themselves  may  be  unphysical.  In  some  difficult  cases,  it  is  a  matter  of  trial 
and  error  to  determine  which  set  of  reference  orbitals  is  best  suited  to  the  problem  - 
the  size  of  the  correlation  corrections  for  each  reference  is  often  the  best  guide. 

One  more  factor  to  consider  is  molecular  geometry.  Experimentahsts  may  report 
vertical  or  adiabatic  ionizations,  or  may  be  unsure.  Proper  vertical  ionizations  can 
only  be  approximated  from  the  optimized  geometry  of  the  N-electron  species,  while 
adiabatic  ionizations  involve  finding  the  potential  energ}'  minimum  for  each  of  the 
A^  -  1-electron  states. 

In  most  cases,  a  single-reference  approach  like  CCSD(T),  an  approximation  to 
the  full  CCSDT  method,  is  significantly  flexible  to  provide  adequate  correlation  to 
all  of  the  states  in  question.  For  calculations  on  certain  open-shell  systems,  however, 
it  becomes  necessary  to  consider  a  multi-determinant  reference.  For  example,  the 
excited  states  of  C2  accessible  from  ionization  of  the  C2  ion[35,  36]  include  open- 
shell  singlet  states.  Attempting  single-reference  calculations  on  these  states  results 
in  significant  spin  contamination,  and  the  splittings  between  these  and  the  related 
triplets  (which  can  be  calculated  from  a  single  high-spin  determinant  calculation) 
£ire  too  large  to  ignore.  For  such  states,  the  multireference  coupled-cluster  (MRCC) 
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method[37,  38]  can  be  used.  In  this  approach,  two  (or  more)  determinants  serve  as 
the  reference  function,  and  an  exponential  operator  (anedogous  to  the  single-reference 
coupled-cluster  theory  described  in  chapter  2)  is  applied  to  each  reference  in  turn 
to  create  the  correlated  wavefunction.  The  number  of  amplitudes  required  grows  as 
the  number  of  unique  determinants  used  in  the  reference;  a  typical  apphcation  of  the 
theory  might  require  two  or  four  reference  determinants.  Since  higher  and  higher 
levels  of  SRCC  theory  are  required  to  obtain  better  spin-projection,  approaching  a 
pure-spin  result,  the  MRCC  method  gives  substantial  savings  in  cases  where  a  full 
triples  or  even  CCSDTQ  approach  might  be  required  to  gain  the  desired  degree  of 
projection. 

6.2     Direct  Approaches 

6.2.1     Two-Hole  One-Particle  Configuration  Interaction 


Configuration  interaction  approaches  can  also  be  applied  to  the  direct  calcula- 
tion of  a  spectrum  of  ionized  states  based  on  a  single  SCF  reference.  Martin  and 
Davidson  [6,  28]  used  a  2hlp  CI  approach  in  which  the  CI  expansion  for  each  of  the 
cationic  states  is  based  on  a  closed-shell  SCF  reference.  The  excitation  operators  used 
in  creating  the  ionized  wavefunction  are  identical  to  those  used  in  the  EOM  method, 
but  the  reference  is  the  SCF  wavefunction.  The  ionized  wavefunctions  thus  may  be 
very  different  from  the  corresponding  EOMIP  wavefunction.  Also,  the  ^ohoCI  ^^' 
ergy  will  be  the  difference,  between  the  correlated  energy  of  the  ionized  state  and 
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the  SCF  energy  of  the  reference  state.  This  means  that  only  the  differences  between 
the  ionization  energies  are  considered  "good" ,  and  it  is  common  practice  to  offset  the 
entire  calculated  spectrum  by  a  constant  energy  in  order  to  match  one  theoretical 
peak  (typically  the  lowest  ionization  energy)  to  its  experimental  value. 


6.2.2     Spin-Adapted  Cluster  Methods 


The  Symmetry  Adapted  Clusters  (SAC)  method  of  Nakatsuji[7,  39]  is  similar  in 
most  respects  to  the  CC  method  presented  in  this  work.  The  important  features 
of  this  approach  are  the  explicit  spin  adaptation  of  the  excitation  operators  and  the 
projection  of  the  final  wavefunction  by  a  spin  projection  operator.  The  result,  though 
requiring  fewer  independent  variables,  has  a  much  more  complicated  algorithmic  ex- 
pression, a  necessity  to  customizing  the  program  for  each  type  of  spin  function,  and  a 
proliferation  of  terms  to  be  coded.  As  a  practical  matter,  the  developer  found  it  nec- 
essary to  truncate  the  exponential  expansion  to  no  more  than  quadratic  terms.  For 
open  shell  references,  this  is  not  quite  equivalent  to  truncation  of  the  CC  equations 
to  T^  terms  as  it  includes  some  higher  spin-flip  terms  as  a  consequence  of  the  spin 
adaptation. 

The  SAC-  configuration  interaction  (SAC-CI)  method[7,  8]  is  very  much  akin  to 
the  EOM  approach,  but  using  explicitly  symmetry-adapted  EOM-type  operators,  and 
based  on  a  SAC  rather  than  a  CC  reference,  resulting  in  a  slightly  different  H,  which 
is  further  altered  by  the  neglect  of  certain  contributions  to  H  elements.  In  addition. 
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implementations  of  SAC-CI  have  also  employed  configuration  selection  procedures 
similar  to  those  used  by  CI  practitioners  to  limit  the  size  of  the  matrices  that  must 
be  diagonalized,  unlike  the  EOMCC  implementation  presented  here. 

6.3    Eclectic  Approaches 

When  the  most  straightforward  methods  fail,  it  is  possible  to  approach  a  difficult 
state  from  an  unusual  direction.  For  example,  inner  valence  ionization  states  can 
be  accessible  as  excitations  from  an  outer  valence  ionized  reference,  or  as  electron 
attached  states  from  a  di-ionized  reference.  Bernholdt  [40],  finding  that  the  structures 
of  the  low-lying  E  states  of  the  C5"  were  difficult  to  impossible  to  converge  using  QRHF 
at  the  CCSD  level  of  theory  as  well  as  evidence  of  discontinuities  in  the  ^E^  state, 
was  able  to  use  the  FSCC(0,1)  approach  with  the  well-behaved,  closed-shell  C5  as 
the  starting  point,  in  order  to  examine  portions  of  the  potential  energy  surfaces  of  all 
three  ionized  states. 


CHAPTER  7 
COUPLED-CLUSTER  LINEAR-RESPONSE  THEORY 

The  formalism  developed  by  Koch  and  J0rgensen[16]  for  coupled-cluster  linesir 
response  functions  may  be  as  easily  applied  to  ionization  potentials  as  to  excitation 
energies.  All  that  is  required  is  to  construct  a  perturbation  operator  capable  of  causing 
an  "excitation"  from  the  occupied  to  "free-electron"  space,  and  a  new  wavefunction  of 
the  appropriate  form.  In  the  following  sections,  a  linear  response  theory  formahsm, 
based  on  the  time-dependent  perurbation  theory  (TDPT)  approach  of  Koch  and 
J0rgensen,  will  be  developed. 

7.1     Perturbation  Applied  to  the  CC  Hamiltonian 

An  RSPT  approach  is  used,  with  the  usual  coupled-cluster  quantities  defined  as 
the  zeroth-order  solutions,  and  the  full  wavefunction  expanded  in  terms  of  corrections 
to  the  coupled-cluster  T  operator.  Although  there  is  a  single  zeroth  order  expression, 
a  spectrum  of  solutions  to  the  full  Hamiltonian  will  be  sought,  the  maximum  number 
of  which  can  be  determined  from  the  form  of  the  corrections  to  T.  The  definitions  to 
be  used  are: 
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1.  i^sCF)  ^^6  zeroth  order  approximation  to  the  Hamiltonian,  is  the  usual  Fock 
operator.  The  perturbation  is  a  first-order  correction  to  the  Hamiltonian,  V^^. 
The  full  Hamiltonian  becomes 


IP 


Hscf  +  pvr 


(7.1) 


2.  The  zeroth  order  wavefunctions  and  energies  are  those  from  the  coupled-cluster 


and  A  solutions; 


|X;(°))    =    |0)  =  |cc°)  =  e^i^' |0) 


-C'O  =      -^1 


cc 


=    E^'^ 


K^'^)  ^  Ecc\fC^'^) 


(0|(l  +  A^)e-^<: 


C    E co- 


rn 


(7.2) 
(7.3) 
(7.4) 

(7.5) 
(7.6) 


3.  The  T  operator  for  the  state  labeled  K  is  written  as  the  zeroth  order  T  plus  a 
sum  of  perturbative  corrections: 


Tk.  ^  ri">  +  tH^  '  ^^'^ 


=  r(o)  +  ri^' +  rf  + 


K     + 

(1)      ,     rp{2) 


(7.7) 
(7.8) 


The  full  wavefunction  can  be  written  as 


IK) 


o'^IC 


HF) 


n(0)  -^(1) 


-(2) 


e  >c    *  e  *=    *  e  'c    * 


|HF) 


7.(0) 
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((i  +  4^'  +  ^4^^^---)*(i  +  :4^'  +  ^rf%...)*---)|HF), 


which  works  out  the  zeroth  order  wavefunction  plus  a  sum  of  perturbative  cor- 


rections 


\)C)  =  \k:^'^)  +  \fc^'^)  +  \k'^)  + 


(7.10) 
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^(2)^  =   (Tf +  ^ri^^')|x:(°)) 


(7.11) 

(7.12) 


etc. 


The  A  operator  can  also  be  expanded  in  a  perturbation  series 


A^  =  A(°)  +  A'^'  + 


so  that  the  full  wavefunction  is 


(£d  =  (0|(l  +  Ar  +  A?  +  ---)e-<'-°'*'-i'+-'. 


a  sum  of  the  zeroth  order  wavefunction  and  perturbative  corrections: 
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(7.14) 
(7.15) 


4.  The  perturbation  operator,  V'^,  must  allow  first-order  corrections  to  T  and  A 
to  be  defined  as  follows: 
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rf    ^    ^t^aU  +  -^Y.i^aUb^J  + 


(7.16) 
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.(1)  ^I 


,{0)     ^     v-x(o).   nu 


Af  -  E4";«)^  (7-18) 

At^'    ^    E^^i(-^)'-  (7-19) 

Here,  as  expected,  is  where  the  treatment  of  the  IP  problem  is  different  from 
that  of  the  EE  problem.  Although  the  TJ^^  amplitudes  have  nonzero  contri- 
butions only  to  the  usual  excitations  from  occupied  to  virtual  space,  a  single 
excitation  to  a  "free-electron"  function  xk  must  be  introduced  in  the  first  order 
to  describe  £in  ionization  process  while  preserving  the  number  of  electrons.  In 
the  limit  where  the  "escaping"  electron  has  zero  kinetic  energy,  this  form  of  Tj^ 
allows  the  the  first-order  corrections  to  the  wavef unctions,  IC^^n  and  (£(^'  ,  to 
be  represented  as  (iV  —  l)-electron  eigenfunctions  of  the  zeroth-order  Hamilto- 
nian.  Keeping  in  mind  that  in  the  zeroth-order  wavefunction  coefficients  to  the 
TJ^Xir^y  operators  are  zero,  while  in  the  first-order  correction  the  coefficients 
to  the  r°\(tau°)^  operators  are  zero,  it  is  possible  to  treat  the  perturbation 
problem  in  a  manner  analogous  to  that  of  Koch  and  J0rgensen. 

So  far,  the  only  assumption  made  in  this  development  of  a  perturbation  theory 
is  that  the  only  processes  initiated  by  the  perturbing  operator,  V^^,  are  electron 
detachment  type  processes.  This  is  justified  in  describing,  for  example,  a  photode- 
tachment  process  where  there  is  a  sizable  gap  between  the  frequencies  of  hght  which 


will  cause  electron  detachment  and  those  which  cause  excitations.    One  does  not 
expect  excitations  and  ionizations  to  occur  at  the  same  frequency. 

7.2     coupled-cluster  Response  Theory 

7.2.1     Time  Dependent  Perturbation  on  a  Coupled-Cluster  Reference 

With  these  definitions,  the  coupled-cluster  response  theory  for  electron  detach- 
ment processes  can  be  developed  in  a  manner  analogous  to  that  of  Koch  and  j0rgensen. 
Although  the  determination  of  the  energies  requires  only  the  Tj^  amphtudes,  the 
formation  of  the  hnear  response  function  requires  the  calculation  of  the  first-order 
corrections  to  the  wavefunction. 

Once  the  CC  solution  is  found,  an  operator  Hj^  =  e"^»^  H<^Q^e^i<^  is  defined, 
which  has  right  and  left  eigenfunctions  |0)  and  (0]  (1  -I-  A^) 

Hn\^)    =  Ecc\0)  (7.20) 

(0|(1  +  A^)^;v    =  {0\{l  +  A'Jc)Ecc,  (7.21) 

or 

i/sCPe^r|0)     =  ^cce^"  |0)  (7.22) 

(0|(l+A°:)e-^ri/scF    =  {0\{1  +  ADe-"^'^' Ecc-  (7.23) 

The  next  step  is  to  find  the  response  of  both  the  left  and  right  eigenfunctions. 
Projecting  the  Schrodinger  equation 

if^\CC{t))  =  H\CC{t))  (7.24) 


by  (//]  e  ^'^  and  collecting  terms  by  order  in  /?  leads  to  the  following  expressions  for 
the  time  evolution  of  T/c'. 
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where  (/i|  =  (/i|  e  -^k  ,  if^^  =  e  '^<:  H^Qpe'^ic  .  The  response  of  the  A  amphtudes, 
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z  (HF|  (1  +  Ax:)  e-^'=(i/scF  +  O^^'', 


can  be  derived  with  the  help  of  equation  7.25 


dhK 
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7.2.2    The  Linear  Response  Function 

Following  the  development  of  Olsen  and  J0rgensen[41],  as  apphed  by  Koch  and 
J0rgensen  to  the  coupled-cluster  reference [16],  we  define  the  operator  perturbation 
operator,  V^^ .  Assuming  the  perturbation  vanishes  at  t  =  — oo,  it  can  be  written  as 

/oo 
^y^g(-«^+a)t^  (7.32) 

-oo 

where  a  real,  positive,  and  small  a  ensures  that  V^^  is  zero  at  t  =  — oo,  and  V'^e°'^ 
is  assumed  to  reach  a  finite  value  at  non-infinite  times.  If  V^  contains  one  constant 


firequency  component  Uf,,  this  becomes 

-ex:' 


(7.33) 
(7.34) 


which  could  describe  the  amplitude  of  a  periodic  electric  or  magnetic  field  of  frequency 


CUb- 


The  expectation  value  of  an  operator  0{t)  is  (£(i)|  0{t)  \CC{t))  and  can  be  written 
as  a  perturbation  expansion  of  C  and  \CC) 
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Or,  in  response  function  notation. 
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So  that 
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The  t  and  A  amplitudes  can  also  be  written  in  terms  of  their  Fourier  transforms: 
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Substituting  these  expressions  into  equations  7.27  and  7.31  gives 
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where 
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Now  the  linear  response  function  can  be  written  in  terms  of  the  Fourier  transform  of 
the  first-order  perturbed  wavefunction: 

m  Vr)U    =    E  >S('^i)  (/i|  O  \CC')  -f  E  (>Cl  [O,  T,]  \CC')  X^^lM  (7.43) 
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+  E  (^1  [O,  T,]  |CC«)  E(-^  +  u;,I);J  e  (7.44) 


7.2.3     Energies  and  Transition  Intensities 


The  transition  energies  are  found  at  the  poles  of  the  Hnear  response  function, 
which  are  found  at  +/—  the  eigenvalues  of  the  matrix  A.  This  matrix, 

A^,^  =  {^l\  [Hn,t^]  |0), 

is  exactly  the  matrix  derived  in  the  coupled-cluster  Equation  of  Motion  method,  by  a 
completely  different  route.  The  r^  operators  are  the  EOMIP  wavefunctions.  There  is 
no  quantity  analogous  to  the  t'^''  amplitudes.  What  the  LRT  provides,  from  the  EOM 
practitioner's  point  of  view,  is  an  unambiguous  definition  of  the  transition  intensities. 


/=    lim  (a;i-a;,)((0;F'))  (7.45) 

where  Uk  is  a  simple  pole  of  ((O;  V*)),  an  eigenvalue  of  the  matrix  A.  If 

{S-'AS)nm  =  Vy„^  =  ^nm^n,  (7-46) 

the  bra  and  ket  functions  can  be  transformed: 


The  X^^'^  and  1'^'^  vectors  can  also  be  transformed: 
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The  hnear  response  function  can  now  be  written 
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and  the  transition  intensity  as 
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are  referred  to  as  the  transition  matrix  elements  and 


(0)  (0)  „(0) 

0    =    6^'=  Oe^'c    =  (Oe^*^  )<- 


In  the  usual  propagator  approach  to  transition  expectation  values,  an  expression 
of  the  form  (0|  O  \k){k\  V^*  |0)  would  be  expected.  The  expression  given  above  does 
not  have  the  same  S3Tnmetric  appearance  because  of  the  difference  between  the  right- 
and  left-hand  solutions.  The  calculation  of  the  transition  expectation  value  is  more 


complicated  because  of  this  asymmetry. 


CHAPTER  8 
METHODS  FOR  TRANSITION  MOMENTS  IN  THE  EOMIP  APPROACH 

For  an  experiment  in  which  it  is  proposed  to  observe  an  absorption  spectrum,  the 
dipole  operator  r  is  the  appropriate  choice  for  O.  However,  in  the  case  of  photoelec- 
tron  experiments,  it  is  the  electron  momentum  which  is  observed.  The  proper  choice 
of  O  must  be  p,  the  momentum  operator  for  the  outgoing  electron.  In  either  of  these 
cases,  (9  is  a  one-body  operator.  The  second  term  of  T^^f.  contains  the  integral 

(0|(l  +  A°)(Or„Tfc),|0), 

which  includes  a  product  of  two  one-body  excitation  operators  (r„  and  r^)  contracted 
with  a  one-body  operator  (O).  If  the  outgoing  electron  functions  are  orthogonal  to 
the  occupied  and  virtual  orbitals,  there  is  no  contraction  of  these  three  operators 
which  can  survive  the  left  projection  by  (0|  (1  -I-  A°),  so  that  term  becomes  zero  and 
the  transition  matrix  element  becomes 

ro^,  =  (0|(l  +  A°)(OT,),|0).  (8.1) 

No  matter  what  is  used  for  O,  the  nature  of  the  outgoing  one-electron  function 
must  be  known  in  order  to  calculate  r^_^^.  and  Fj^^g-  Free-electron  functions  are 
dependent  on  the  experimental  method  and  not  included  in  the  basis  sets  for  standard 
electronic  structure  methods;  integrals  of  f  and  p  involving  these  functions  are  not 
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readily  available.  Therefore,  it  is  usual  to  make  some  kind  of  approximation  for 
both  Tq_^i^  and  V^II^q.  The  most  common  approximation  is  the  so-called  monopole 
approximation  as  developed  by,  for  example,  Martin  and  Shirley  in  1976[4]. 

8.1     The  Monopole  Approximation 

A  symmetric  Hamiltonian  is  assumed,  and  the  dipole  approximation  to  the  cross 
section  for  photoemission  by  an  incident  photon  field  is  employed: 

GDA  oc  (hw,)-'  [f  ■  {i;f{N)\p\i;i{N))]'p{Ef),  (8.2) 

where  p{Ef)  is  the  density  of  states  in  the  continuum  at  the  final  energy.  The  authors 
proposed  to  determine  the  intensities  of  the  shaJteup  satellite  peaks  relative  to  the  re- 
lated principal  core  ionization  potentials,  introducing  two  important  approximations. 

1.  The  principal  ionizations  are  characterized  by  annihilation  of  an  electron  from 
a  single  orbital,  and  the  satellites  also  have  a  non-negligible  contribution  from 
this  annihilation.  Therefore  the  momentum  operator,  acting  on  an  initial  state 
dominated  by  an  ionization  from  orbital  i,  is  approximated  as: 

E  ("^i  P  \X)  m^X  \MN))  ~  ^^X  mN))  {i\  P  |X>  •  (8.3) 

m 

2.  Both  p  and  the  momentum  integrals  (i|p|x)  are  approximately  constant  over 
the  range  of  energies  spanned  by  the  principal  and  its  satellites. 


This  means  that  the  ratio  of  the  intensity  of  the  n*'*  sateUite  peak,  I{n),  to  its  prin- 
cipal, 7(0),  is 
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And  the  calculation  of  the  intensities  is  reduced  to  the  calculation  of  transition  density 
matrix  elements  for  the  principal  and  its  satellites. 

Further  simplifications  are  possible.  In  1977,  Martin  and  Davidson[6]  employed  a 
simplification  which  included  only  final  state  correlation  and  only  one  term  from  the 
resulting  transition  density  matrix.  For  the  ionized  states,  they  used  a  (Ih,  2hlp)  CI 
calculation  and  argued  that,  for  the  principal  and  its  satellites,  only  one  contribution 
was  important  —  that  hole  state  which  dominates  both.  So  for  a  CI  wavefunction  of 
the  form 

^„  =  j:Cf$.-f^C%$r,.  (8.5) 

i  ija 

they  approximated  the  ratio  of  intensities  as 

/(n)//(0).<^"l''l°>Nig4.  (8.6) 

In  a  1992  paper  by  Murray  and  Davidson [28],  the  same  features  were  calculated  with 
a  multi-reference  CI  approach,  and  the  same  expression  was  used  for  the  relative 
intensities.  The  authors  report  good  agreement  with  experiment. 


8.2    Transition  Moment  Expression  for  Non-Svmmetric  Eigenstates 


For  the  non  symmetric  eigenvalue  problem,  there  are  two  different  approximations 
to  the  wavefunction  in  the  full  space.  Equation  8.2  must  be  modified  to  recognize 
that  the  left-  and  right-  handed  wavefunctions  are  not  identical.  This  results  in 

GDA  oc  {hu,)-'  [r  •  (^^f  (iV)|p  |^;(iV))]  [f  •  (tAj'(iV)|p|V;f  (iV))]  p{Ef),        (8.7) 


which  yields,  under  the  same  set  of  approximations  as  was  made  above, 
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In  this  case,  the  Murray  and  Davidson  approximation  gives 
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In  cases  where  Hf^  is  close  to  being  symmetric,  there  is  some  justification,  consid- 
ering the  number  of  approximations  already  made,  in  reducing  this  back  to  equation 
8.6.  This  approximation  can,  of  course,  only  be  applied  to  those  cases  where  one  of 
the  states  of  interest  is  a  principal  ionization  state. 

8.3     Transition  Moment  Expression  for  the  EOMIP  problem 

For  the  EOMIP  case,  the  expressions  for  the  transition  matrix  elements  are: 
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Figure  8.1:  Diagrammatic  expression  of  Fo'.^^t  terms 
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The  Other  matrix  element  is: 
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Figure  8.2:  Diagrammatic  expression  of  FJ^Zo  terms 
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In  cases  where  the  outgoing  electron  function  is  orthogonal  to  the  virtual  orbitals, 
all  F  terms,  which  involve  (x|  p)  or  (p|  x)  will  disappear,  leaving  the  much  simpler  D 
expressions  in  terms  of  (p|  O  \x)  and  (xi  O  \p).  Using  the  approximations  of  Martin 
and  Shirley,  these  expressions  reduce  to 
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(8.12) 
(8.13) 


where  i  is  the  occupied  orbital  which  dominates  the  ionization  event. 


That  means  that  the  ratio  of  intensities  of  satellite  state  A/"  and  its  principal  state 
K,  is  approximated  as. 

^    ) ^^ ^.         (8.14) 
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Because  Y'^  is,  to  first  order,  identical  to  C'^,  one  can  make  the  further  approxi- 
mation 

:W      ^  V  ^  rnenf  J 

i{}C)    ~       )  ^^     ^  1        ^,;  ,/  ■       ^^   ^ 

\  fne  ■^  men/ 

The  simplest  possible  approximation,  then,  to  EOMIP  relative  intensities  is 


C; 


and  the  next  simplest  approximation  is: 
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One  more  approach  is  simply  to  use  the  approximations  to  the  transition  matrix 
elements  directly,  as  electron  propagator  practitioners  often  use  pole  strengths  to 
approximate  intensities.  This  procedure  is  most  defensible  when  the  pole  strengths 
are  compared  to  low-incident  energy  photoelectron  spectra. 

lW^C^U-^Y.{^r^-^mti)K]-  (8.18) 


CHAPTER  9 
IMPLEMENTATION 

The  EOMIP  equations  have  been  implemented  in  the  ACESII  framework,  us- 
ing the  direct  product  decomposition  approach  of  Stanton,  Gauss,  and  Bartlett. 
Rather  than  expHcitly  forming  the  matrix  defined  by  equations  3.20  and  3.21,  the 
non-symmetric  generalization  of  the  Davidson[24]  method  developed  by  Hirao  and 
Nakatsuji[42]  for  eigenvector  extraction  has  been  used.  In  this  approach,  matrix- 
vector  products  are  formed  and  stored  on  disk,  and  only  a  smaller,  projected  "mini- 
matrix"  is  kept  in  memory. 

It  is  computationally  convenient  to  form  a  small  matrix  in  memory  and  diagonalize 
it  to  recover  all  the  eigenvectors  and  eigenvalues.  But  in  many  cases,  only  a  few  of 
these  solutions  are  of  interest,  and  the  diagonalization  becomes  increasingly  expensive. 
Therefore,  it  is  advantageous  to  use  methods  in  which  the  matrix  to  be  diagonalized 
is  never  formed  explicitly,  or  at  least  never  held  in  memory,  and  the  matrix  is  never 
diagonahzed.  All  of  these  large  matrix  methods  depend  on  the  fact  that  the  Raleigh 
quotient  of  a  non-symmetric  matrix  A 

P(y,x)  =  ^  (9.1) 

is  stationary  with  respect  to  variation  of  the  elements  of  x  and  y  if  and  only  if  these 
vectors  are  corresponding  left  and  right  eigenvectors  of  A. 
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Davidson's  method,  like  many  other  large  matrix  methods,  uses  a  variation  on  the 
Lanczos  method.  In  both  these  methods,  the  nonsymmetric  eigenvalue  problem 

AXi    =    XiX,  (9.2) 

Y^A    =    A,r/  (9.3) 

is  solved  by  projecting  the  nonsymmetric  matrix  A  into  a  smaller  space  spanned  by 
a  set  of  trial  vectors  B.  Since  the  matrix  is  nonsymmetric,  it  is  formally  correct  to 
use  two  different  bi-orthogonal  trial  spaces,  B  and  L.  The  projected  "minimatrix" 

M  =  VAB  (9.4) 

L^B  =  I  (9.5) 

is  diagonalized  to  yield  approximate  eigenvalues  and  eigenvectors  which  are  expansion 
coefficients  to  approximate  eigenvectors  of  A 

D^MC    =    d  (9.6) 

X,    =    ^Q,B,  (9.7) 

k 

y;  =  ^Afci^-.  (9.8) 

k 

Residual  vectors  X  and  y  are  a  measure  of  the  error  in  the  approximate  eigenvector. 
These  error  vectors  are  defined  to  satisfy 

{A-di){X,-Xi)  =  0  (9.9) 

{A-d,){Y,-y,)  =  0.  (9.10) 

It  is  the  method  of  approximating  the  residual  vector  which  characterizes  these  meth- 
ods. The  Davidson  method  approximates  these  as 

X,,    =    {d,-A,j)-'[{AX),,-diX,,]  (9.11) 
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y,,  =  (d.-.4,,)-i[(,4^r),.-rf.r,,].  (9.12) 

These  new  correction  vectors  are  then  bi-orthonormahzed  and  added  to  the  B  and  L 
spaces,  after  which  the  entire  procedure  is  repeated. 

If  only  the  right-hand  eigenvectors,  A',  are  required,  there  is  another  approach 
which  can  be  used  to  save  computational  time.  The  right-hand  projection  space  B 
can  be  used  in  place  of  the  L  space,  and  the  method  becomes  eilmost  identical  to 
the  symmetric  Davidson  method  except  that  the  algorithm  for  diagonalization  of  the 
minimatrix  M  must  be  suitable  for  nonsymmetric  matrices.  By  the  time  convergence 
is  reached,  B  generally  spans  a  space  sufficient  to  represent,  to  a  good  approximation, 
the  left-hand  eigenvectors  of  A.  In  test  calculations,  these  two  approaches  have  been 
shown  to  have  very  similar  convergence  proper  ties  [42]. 

At  each  iteration,  it  is  necessary  to  choose  some  subset  of  roots  of  the  minimatrix 
to  be  improved.  There  are  several  approaches.  If  one  is  interested  in  the  A^  lowest 
roots  of  the  matrix  A,  one  can  simply  always  choose  the  lowest  roots  of  the  minimatrix. 
If,  however,  it  is  desired  to  converge  directly  to  higher  roots,  some  other  approach 
must  be  used. 

Because  the  principal  IP  eigenvectors  are,  in  almost  all  cases,  dominated  by  the 
annihilation  of  an  electron  from  a  single  occupied  molecular  orbital,  it  is  useful  to 
"pre-load"  the  B  space  with  guess  vectors  having  C^  —  1  for  each  root  K,  desired 
corresponding  to  the  ionization  of  the  electron  in  orbital  k.  Then,  at  each  iteration, 
those  minimatrix  roots  which  are  most  similar  to  the  original  guess  vector  are  chosen 
for  improvement  until  convergence  is  reached.    Alternatively,  those  roots  which  are 
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most  similar  to  roots  saved  from  the  previous  iteration  can  be  chosen,  but  the  au- 
thor has  found  that  the  use  of  this  approach  can  lead  to  choosing  the  wrong  root,  as 
the  desired  root  may  be  poorly  represented  in  the  first  few  iterations,  and  the  algo- 
rithm may  follow  several  different  roots  successively,  finally  settling  on  an  unexpected 
solution. 

Shakeup  roots  which  are  characterizable  as  satellites  of  a  particular  principal  ion- 
ization are  still  dominated  by  the  annihilation  of  a  single  electron,  but  aJso  contain 
large  contributions  from  one  or  several  excitations  of  the  ionized  species.  These  so- 
lutions are  more  difficult  to  characterize  for  two  reasons:  1)  it  is  not  always  known 
ahead  of  time  which  excitations  are  most  important  to  the  desired  shakeups,  and 
2)  even  in  cases  where  the  important  excitations  are  known,  the  unoccupied  orbitals 
available  from  a  typical  Hartree-Fock  calculation  in  a  larger  than  minimaJ  basis  set  are 
not  useful  in  a  "molecular  orbital"  analysis  of  the  system  being  studied  -  they  have 
very  little  correspondence  to  the  orbitals  which  would  be  produced  in  an  excited-state 
or  N-l-1  (or  more)  electron  calculation  which  occupied  portions  of  the  desired  orbital 
space.  In  this  case,  it  is  useful,  once  the  principal  root  is  found,  to  specify  an  energy 
range  which  removes  the  principal  root  from  consideration,  and  then  converge  all 
roots  found  within  that  range,  or  impose  the  additional  constraint  that  chosen  roots 
should  have  a  significant  overlap  with  the  principal  root  of  which  they  are  satellites. 
In  the  author's  experience  that  first  converging  the  principal  roots  in  question  results 
in  a  B  space  containing  significant  contributions  to  the  space  of  shakeup  roots.  Sev- 
eral of  these  satellites  near  in  energy  to  the  principal  may  be  selected  within  the  first 
few  shakeup  iterations,  especially  in  cases  where  these  shakeups  zire  large  compared 


to  the  "principal" ,  and  the  two  wavefunctions  are  similar.  In  other  cases,  the  author 
uses  other  procedures  to  load  the  B  space  with  relevant  vectors,  such  as  relaxing 
energ}'  range  restrictions  for  a  few  iterations,  even  choosing  random  vectors. 


CHAPTER  10 
APPLICATION  OF  THE  EOMIP  METHOD  TO  SMALL  MOLECULES 

10.1     A  Survey  of  FSCCSD  and  EOMCCSD  Calcnlations 

Table  10.1  is  a  survey  of  FSCCSD  for  the  ionization  of  a  range  of  small  molecules. 
Bernholdt  [40]  reported  results  differing  by  almost  zero  to  approximately  .6  eV  from 
experiment  for  valence  ionizations  of  a  range  of  molecules,  mostly  in  DZP  basis  sets, 
mostly  at  experimental  geometries.  There  are  a  few  general  trends. 

•  Energies  in  a  larger  basis  at  the  same  geometry  are  larger. 

•  DZP  ionization  energies  are  generally  smaller  than  experimental  energies. 

•  For  well-behaved  closed-shell  molecules,  the  first  ionization  energy,  even  at  the 
DZP  level  is  generally  within  .2  eV  or  so  of  experiment.  A  notable  exception  to 
this  rule  is  the  s-tetrazine  example. 

•  In  the  larger  basis  sets,  there  is  a  trend  for  higher-energy  EOMIP  ionizations 
to  be  increasingly  high  relative  to  their  corresponding  experimental  energies. 

•  Increasing  the  basis  set  size,  as  in  the  cases  of  the  Stanton,  Bartlett,  and  Rittby 
O2  study  and  the  Pal  et  al.  study  of  formaldehyde,  improves  agreement  with 
experiment. 
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Table  10.1:  A  survey  of  FSCCSD  results 


Final 

Energy  (eV) 

Molecule 

Basis 

State 

Theory 

Expt. 

CH2NH 

DZP' 

''A' 

10.34 

10.52, 

10.70,  10.56 

M" 

12.29 

12.43, 

12.48,  12.44 

M' 

15.08 

15.13, 

15.11,  15.00 

'A' 

17.19 

17.04, 

17.07,  17.00 

s-tetrazine 

DZP'^ 

'B^g 

^B2g 

'B,g 

9.20 
11.87 
12.12 
12.69 
13.01 
13.52 

9.7 
11.9 
12.1 
12.8 
13.3 
13.5 

H2CO 

[5s3pld/2slp]<' 

'B2 

'B, 
Ml 

10.51 
14.36 
15.75 

10.88 
14.38 
16.00 

H2CO 

[5s3j)/2s]V 

'B2 

10.07^ 

10.43'' 

10.88 

[5s3pld/2slpY 

'B, 

14.14^ 

14.13'^ 

14.38 

Ml 

15.20^ 

15.67'' 

16.00 

CH2 

DZP^ 

Ml 

10.20 

10.2e 

;  (DZP-FCI) 

'B, 

14.91 

14.85  (DZP-FCI) 

Ml 

22.55 

22.14 

1  (DZP-FCI) 

CH2PH 

DZP' 

M" 
M' 
M' 

10.08 
10.28 
13.17 

10.30 
10.70 
13.20 

^Sp-  O2 

DZPVPVQZ++9 

^n. 

11.76°, 

12.39» 

12.35 

'Uu 

16.40", 

16.76» 

16.85 

^E; 

17.75^ 

18.265 

18.33 

4^_ 

24.58°, 

24.82» 

24.66 

N2 

[5s4pl(fl'' 

^g 

15.44 
17.14 
18.64 

15.6^ 

16.88'' 

18.6'' 

ketene 

DZP" 

'B, 
'B2 

9.40 
14.26 

9.64(a) 
13.84(a) 

H2O 

DZP= 

'B, 
Ml 
'B, 

11.97 
14.28 
18.75 

12.61 
14.73 
18.55 

diazomethane 

DZP" 

'Br 
'B, 

8.60 
13.79 

9.00(a) 
14.13(v) 

HF 

[5s3p/3s]/ 

2n 

2v; 

15.08 
19.43 
39.75 

16.1 
19.9 
39.7 

"  Bernholdt  [40];  *"  Experimental  results  [34]  '^  Vaval,  Ghose,  Pal,  and  Mukherjee  [43] 

^  Pal,  Rittby,  Bartlett,  Sinha,  Mukherjee  [44],  [5s3p/2s]  basis;  ^  ibid,  [5s3pld/2slp]  basis 

(a)  and  (v)  indicate  adiabatic/vertical  ionizations  of  ketene  and  diazomethane 

f  Smha  et  al.  [45] 

*  Stanton,  Bartlett,  and  Rittby  [46] 

''  Rittby  and  Bartlett  [13] 

'  Watts,  Rittby,  and  Bartlett  [47] 


These  results  are  encouraging;  it  is  easy  to  surmise  that  the  EOMIP  approach  can 
be  appHed  to  a  variety  of  molecules  with  very  good  results.  In  addition,  the  EOMIP 
method  can  be  used  to  examine  shake-up  states,  and  is  not  plagued  (as  is  FSCC)  by 
intruder-state  problems  in  the  inner- valence  region. 

The  EOMIP  approach  is  conceptually  well  suited  to  the  study  of  ionic  spectra, 
including  first,  principal,  and  shakeup  ionization  energies  and  intensities.  For  N2,  for 
example,  there  are  six  molecular  orbitals  (five  a  and  one  tt)  to  be  explored.  Using 
a  well-correlated  single-state  method  such  as  CCSD  or  CCSD(T)  would  require  one 
calculation  for  the  neutral  N2,  and  a  separate  calculation  for  each  of  the  N2  states. 
All  the  Nj"  calculations  would  have  to  be  open-shell,  doubling  the  effective  size  of 
the  orbital  space  involved  in  the  SCF  and  correlated  calculations.  In  addition,  it 
is  not  easy  to  obtain  a  suitable  set  of  molecular  orbitals  for  every  open-shell  state. 
While  it  is  a  simple  matter  to  remove  an  electron  from  the  highest-energy  orbital 
of  any  symmetry,  it  requires  a  bit  of  trickery  and  a  few  more  SCF  iterations  to 
converge  on  a  core-hole  state.  And  it  is  generally  impossible  to  converge  an  SCF 
state  with  an  electron  removed  from  an  inner- valence  orbital.  The  QRHF  reference, 
a  non-SCF  state  can  be  used  in  this  case.  There  is  no  reliable  way  to  describe  a 
state  corresponding  to  a  shakeup  peak,  and  no  way  to  specify  which  of  the  two  2<Tj 
nitrogen  peaks  (table  10.4  and  figure  10.1)  will  be  found  in  this  approach.  When  the 
EOMIP  approach  is  used,  only  one  CCSD  calculation,  on  the  closed-shell  state,  is 
required.  The  EOMIP  calculations  which  follow  are  considerably  less  expensive  in 
terms  of  computer  time.  And  shakeup  peaks  can  be  specified  in  terms  of  energy  ranges 
and  the  molecular  orbitals  with  with  they  are  associated.    The  equation-of-motion 


excitation  energy  (EOMEE)  approach  can  also  be  used  as  a  route  to  finding  ionization 
states  in  the  excitation  spectrum  of  N2'  species.  There  are  two  disadvantages  to  this 
approach.  First,  not  all  of  the  desired  ionization  states  can  be  described  primarily  as 
a  single  excitation  from  one  reference.  Multiple  references,  eind  thus  multiple  CCSD 
calculations  must  be  employed.  Secondly,  although  the  linear  form  of  the  EOMEE 
equations  makes  iteration  less  computationally  demanding  than  CCSD  iterations, 
EOMEE  is  still  more  demanding  than  EOMIP  -  there  are  more  amphtudes  and  more 
and  bigger  terms  to  calculate  in  EOMEE. 

What  follows  in  this  chapter  is  an  EOMIP  study  of  the  ionization  spectra  of 
several  well-understood  molecules,  using  large  basis  sets,  with  attention  to  prominent 
shakeup  features.  All  calculations  were  performed  with  the  ACES  II  program  system^ 
and  the  author's  ACES  Il-based  CC  EOMIP  program^. 

The  EOMIP  intensities  which  appear  in  this  and  the  following  chapter  were  calcu- 
lated using  the  approximation  of  equation  8.18,  except  where  principal  peak  intensities 
have  been  scaled  to  experimental  intensities,  as  noted.  For  the  figures,  Lorentzian 
functions  were  applied  to  the  theoretical  transition  energy/intensity  information  in 
order  to  simulated  the  degree  of  experimentally-observed  broadening.  In  some  cases, 
the  experimental  curves  were  also  created  from  available  energy /intensity  data  by 


^ACES  II  is  a  quantum  chemical  program  package  especially  designed  for  CC  and  MBPT  en- 
ergy and  gradient  calculations.  Elements  of  this  package  are:  the  SCF,  integral  transformation, 
correlation  energy,  and  gradient  codes  written  by  J.  F.  Stanton,  J.  Gauss,  J.  D.  Watts,  W.  J.  Laud- 
erdale, and  R.  J.  Bartlett;  the  VMOL  integral  and  VPROPS  property  integral  programs  written  by 
P.  R.  Taylor  and  J.  Almlof;  a  modified  version  of  the  integral  derivative  program  ABACUS  written 
by  T.  Helgaker,  H.  J.  Aa.  Jensen,  P.  J0rgensen,  J.    Olsen,  and  P.  R.  Taylor[48] 

•^CCEOMIP  is  a  quantum  chemical  program  designed  to  calculate  vertical  electron  detachment 
energies  ajid  intensities,  written  by  Renee  Peloquin  Mattie.  It  requires  ACES-II  CCSD  T  and  H 
lists  [49] 


application  of  Lorentzian  lineshapes.   In  the  carbon  cluster  plots,  the  experimental 
curves  were  created  by  digitization  of  experimental  plots  provided  by  Yang  et  al[35]. 

10.2     Unsaturated  Molecules 

10.2.1     Nitrogen 


Nitrogen  presented  one  of  the  earliest  challenges  to  theoreticians.  The  ionization 
energies  obtained  from  Koopmans'  theorem  are  not  predictive  the  observed  spectrum, 
and  are  in  fact  incorrectly  ordered  (table  10.3).  Self-consistent  field  calculations  do 
not  order  these  ionizations  properly.  Correlated  calculations  are  necessary,  and  highly 
correlated  calculations  are  needed  to  describe  the  breakdown  of  the  molecular  orbital 
picture  in  the  inner- valence  region  correctly. 

Two  basis  sets  were  used  in  the  nitrogen  calculations.  The  smaller  basis  was  based 
on  the  triple-zeta  plus  polarization  basis  set  of  Dunning [50],  with  the  innermost  p  or- 
bital uncontracted,  as  was  used  by  Rittby  and  Bartlett[34].  The  larger  is  the  Dunning 
J9VTZ+  basis  set  [51],  a  polarized- valence  triple-zeta  basis  set  which  includes  diffuse 
s,p,d  and  f  functions  .  Here,  however,  the  CCSD  optimized  internuclear  separation, 
R  =  2.0778  was  used  rather  than  the  experimentally-obtained  geometry  R  =  2.0693 
used  by  Rittby  and  Bartlett. 

In  the  smaller  basis  set  (table  10.2),  results  show  agreement  with  experimental 
results  to  within  .2  eV  for  the  outer- valence  region.  The  2og  inner  valence  hole  state, 
however,  shows  an  error  of  1.3  eV  and  the  lUg  core  hole  state  is  in  error  by  .7  eV. 


Table  10.2:  N2  in  5s4pld  basis  at  /?  =  1.09898  A 


N2 

N- 

Reference 

EOMIP 

0-theo) 

Exp" 

xii:+ 

RHF  -109.358913  au 

X^e;^ 

3ag 

15.376  eV 

15.5  eV 

^ri„ 

ItTu 

17.022  eV 

16.8  eV 

'^: 

2a„ 

18.607  eV 

18.6  eV 
25.0  eV 
28.2  eV 
31.5  eV 

(shakeup) 
(shakeup) 
(shakeup) 

'K 

2ag,  I7rjl7ry2cru 

32.153  eV 

(.5") 

'^: 

2ag 

38.616  eV 

(1.0") 

37.3  eV 

'4 

20g,      bal'^GgZOg 

42.543  eV 

(.2") 

'^i 

lag 

410.623  eV 

409.9  eV 

"  Banna  and  Shirley  [52] 

''  Simple  approximation  relative  to  2ag  principal 

Moving  the  the  larger  pVTZ-l-  basis  seems  to  afford  no  improvement  of  the  valence 

and  core  peaks,  in  fact  shifting  each  peak  toward  higher  energies  by  an  amount  which 

increases  with  the  transition  energy.  As  shown  in  table  10.3,  the  valence  Og  and  a^ 

holes  gire  stiU  within  .2  eV  of  experiment,  but  the  energy  of  the  Itt^  hole  seems  to 

have  gotten  worse  -  it  is  not  quite  within  .5  eV  of  experiment.   However,  the  inner 

valence  region  is  now  much  more  similar  to  the  experimental  spectrum.   The  inner 

valence  2ag  hole  ionization  is  characterized  experimentally  by  a  fairly  broad  peak, 

relative  to  the  outer  valence  or  core  ionization.  In  this  region,  the  simple  MO  picture 

breaks  down,  and  it  is  difficult  to  speak  of  a  single  line  corresponding  to  the  ejection 

of  the  2ag  electron.  The  calculations  show  two  lines  of  similarly  large  intensity,  one 

near  33  and  one  near  39  eV.  Vibrational  and  vibronic  effects  can  be  expected  to  be 

quite  strong  in  this  region  [3l]  as  well. 


Table  10.3:  N2,  pVTZ+  basis 


Koop. 
approx 


Ionization  energy  eV  (Intensity) 
EOMIP  Other"       Exp 


CLOSED 


CORE 


Orbital 
Assignment 


17.277 

15.683 

;.931) 

15.499 

(.880) 

15.666" 

15.5 

2>ag 

16.793 

17.332 

,.958) 

17.159 

(.918) 

17.227 

16.8 

ItTui 

16.793 

17.332 

;.961) 

17.159 

(.964) 

17.227 

16.8 

iTTuy 

21.147 

18.888 

,.895) 

18.674 

(.858) 

18.908 

18.6 
25.0 

2a„ 

— 

29.887 

,.02) 

28.917 

(.026) 

28.2 

2cru,  2iT^Zagl'Ku 
2ag,  27r|2cT„l7r„ 
l7r„,  37r]  l«7u  ScTg 

— 

32.679 

,.320) 

31.983 

(.284) 

29.0734 

31.5 

— 

33.319 

0.004) 

40.189 

38.695     ( 

,.518) 

38.106 

(.470) 

37.3 

2ag 

— 

39.307     ( 

.002) 

2a g,  Aa^Zag'iag 
2a g,  4ap7fu37ru 
2a g,  baHag'iag 
2a g,  27fj3ag27r„ 

— 

42.400     ( 

.034) 

— 

42.519     ( 

.001) 

— 

43.030     ( 

.045) 

— 

43.461     ( 

.067) 

2ag,  4a^37fu37ru 

44.566     ( 

.067) 

2cr<,  +  3(T<,,  27f]^l7f„3a3 



45.452     ( 
46.919     ( 

.001) 
.005) 

2cru,  6cr|^l7r„l7ru 
2(T„,  ha'^^2>ag2>ag 

426.723 

410.587     ( 

.806) 

409.609 

(.753) 

409.9 

ICTu 

426.623 

410.768     ( 

.807) 

409.799 

(.755) 

lag 

"  EOMEE, 
"  ACCSD 


except  where  otherwise  noted 
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In  the  plots  in  figure  10.1,  Lorentzian  lineshapes  have  been  apphed  to  the  EOMIP 
peaks  to  simulate  the  Une  broadening  found  in  low-resolution  experimental  spectra. 
Each  EOMIP  plot  is  superimposed  on  a  plot  of  an  experimentally-derived  photoe- 
mission  spectrum  based  on  MgKa  [52]  and  YMi^  [53]  experiments.  Peak  positions 
were  taken  from  the  YM^  experiment  (at  1253.6  eV),  which  provided  much  better 
resolution  in  the  region  between  20  and  45  eV,  while  relative  peak  intensities  were 
taken  from  the  MgKa  experiment  (at  132.3  eV)  ,  which  is  more  suitable  for  compar- 
ison to  EOMIP  intensities,  as  they  should  most  resemble  threshold  photodetachment 
intensities. 

The  simulated  spectra  in  figure  10.1  indicate  a  good  match  between  the  experi- 
mental peaks  at  28.2  and  31.5  eV  and  the  2cr„  peaks  at  29.9  and  32.7  eV  with  pVTZ-h 
closed  shell  orbitaJs  and  28.9  and  31.9  eV  with  pVTZ-l-  core  hole  orbitals.  The  core 
hole  orbitaJs  give  somewhat  better  agreement  over  the  entire  spectrum,  especially 
at  higher  energies.  In  the  inner-valence  region  and  for  the  core  peaks,  improvement 
becomes  notable,  as  seen  in  table  10.5  No  peaks  were  found  which  can  reasonably 
be  assigned  to  the  25.0  eV  experimental  peak.  Schirmer,  Cederbaum,  Domke,  and 
von  Niessen,  in  their  1977  paper,  saw  a  peak  right  at  25  eV,  but  they  were  using  a 
much  smaller  basis  set  (lis  7p/5s  3p),  which  is  lacking  in  polarization  functions  com- 
pared to  the  pVTZ+  (lis  6p  3d  2f/5s  4p  3d  2f).  The  author's  experience  and  that 
of  others  indicates  that  this  means  some  larger  shakeup  peaks  in  the  smaller  basis 
set,  and  Schirmer  et  al.  show  more  shakeup  peaks  than  are  perhaps  indicated  by  the 
experimental  spectra,  making  a  doublet-with-a-shoulder  of  the  37  eV  2ag  ionization, 
and  giving  them  a  significant  2ag  peak  at  about  40  eV,  similar  to  the  appearance  of 


Table  10.4:  Summary  of  N2  EOMIP  ionization  calculations 


Energies  eV 

5s4p1d 

pVTZ+ 

C.S. 

C.H. 

Exp 

15.376 

15.683 

15.499 

15.5 

17.002 

17.332 

17.159 

16.8 

18.607 

18.888 

18.674 

18.6 
25.0 

29.887 

28.917 

28.2 

32.153 

32.679 
33.319 

31.983 

31.5 

38.616 

38.695 
39.307 
42.400 
42.519 
43.030 
43.461 
44.566 
45.452 

38.106 

37.3 

410.587 

409.609 

409.9 

410.623 

410.768 

409.799 

Table  10.5:  N2  EOMIP  error,  5s4pld  and  pVTZ+  bases. 
Energies  and  errors  are  in  eV. 


Expt. 

5s4pld 

pVTZ+ 
closed     core 

15.5 

-0.1 

0.2       0.0 

16.8 

0.2 

0.5       0.4 

18.6 

0.0 

0.3       0.1 

28.2 

— 

1.0       0.7 

31.5 

0.6 

1.2       0.4 

37.3 

1.3 

1.4       0.8 

409.9 

0.7 

0.8       0.2 

Ng  Ionizations.  EOMIP 


1^     i^     ^     S     S     ze     28     3E     S     3*     ^     S     X     J     ^ 


Ng  Ionizations,  EOMIP 


>i  OMIP 


0.9- 

0.8- 

0.7- 

0.6- 

O.! 

O.' 

0.3- 

0.2- 

O. 


2o: 


1 3    it    75    ?6    i5    S    S    S    S    S    3E    S    S    *    S    45    35    4^ 


Ionization  Energy  (eV) 
No  ionizations,  EOMIP 


3o;     1«u    ^aZ 


2o* 


ze      ^      32      34      3e 


40      ^      *« 


Figure  10.1:  N2  ionizations  in  three  MO  sets 


significant  EOMIP  peak  at  38.6  eV  when  the  5s4pld  basis  is  used.  This  region  shows 
much  less  intensity  in  the  pVTZ+  closed-shell  or  core-hole  orbital  calculations,  in 
agreement  with  experiment.  In  the  core-hole  orbital  set,  the  largest  of  the  2og  peaks 
has  moved  .8  eV  closer  to  the  experimental  peak,  while  its  companion  peak  has  moved 
.3  eV  closer  to  experiment  and  the  ratio  of  the  two  peaks  is  more  similar  to  that  in  the 
experimental  spectrum.  In  the  core  region,  the  loy,  and  lOg  peaks  have  both  moved 
slightly  closer  to  the  experimental  value.  In  general,  the  energies  of  all  the  major 
peaks  are  reduced  in  the  core-hole  basis,  moving  closer  to  experiment.  Interestingly, 
core-hole  orbitals  do  not  improve  the  performance  of  EOMIP  by  decreasing  the  size  of 
correlation  corrections  to  the  wavefunction.  The  sizes  of  the  EOMIP  doubles  contri- 
butions in  the  two  cases  are  similar,  although  there  is  some  difference  in  the  ordering 
of  these  contributions  because  the  largest  Tl  amplitudes  are  approximately  twice  as 
large  (.075  versus  .015)  in  the  core-hole  as  in  the  closed-shell  basis. 


10.2.2    Ethvlene 


Ethylene  has  been  studied  in  two  different  basis  sets  -  Nakatsuji's  Basis  I  [54] 
and  the  Dunning  pVTZ-i-  basis[51].  In  the  smaller  basis,  the  Herzberg[55]  ground 
state  geometry  was  used  and  relative  intensities  were  calculated  using  the  simple 
approximation  of  equation  8.16.  For  the  larger  basis  set  study,  the  geometry  was 
optimized  at  the  CCSD  level  in  the  pVTZ  basis  set,  and  the  improved  approximation 
to  relative  intensities  from  equation  8.17  was  used. 


The  inner-  and  outer- valence  region  of  the  ethylene  ionization  spectrum  has  been 
studied  before,  with  an  eye  toward  the  assignment  of  the  strong  satelhte  peak  at  27.4 
eV[6,  28,  56,  54].  The  valence  ionizations,  as  well  as  the  strongest  of  the  related  shake- 
ups,  are  given  in  table  10.6  along  with  the  experimental  numbers  and  the  theoretical 
results  of  Nakatsuji  in  the  same  basis  set.  The  outer-valence  results  show  excellent 
agreement  with  experiment,  while  the  inner- vaJence  principle  peak  at  w  24  eV,  like 
the  SAC-CI,  appears  at  a  slightly  higher  energy.  The  shakeup  at  27.4  eV  appears 
in  the  experimental  spectrum  as  a  strong,  somewhat  broad  satellite.  Cederbaum[56] 
implies  that  this  peak  should  be  assigned  to  a  2ag  shakeup,  and  in  this  Davidson[28] 
concurs.  Nakatsuji  attributes  this  peak  to  a  sum  of  the  nearby  lb2u  and  2biu  shake- 
ups.  The  difference  in  the  ordering  of  the  lb2u  and  2ag  peaks  just  below  the  the 
experimental  peak  is  puzzling  for  two  methods  which  are  apparently  so  similar,  but 
there  are  several  other  theoretical  shakeup  peaks  with  very  small  intensity  in  the  same 
region.  It  would  seem  that  the  nearby  2aj  shakeup  is  another  hkely  candidate. 

It  is  interesting  to  note  that  both  Cederbaum  and  Davidson  report  significantly 
greater  relative  intensities  for  the  shakeup  peaks  than  have  been  obtained  from  either 
the  EOMIP  or  SAC-CI  calculations.  Davidson  has  noted  that  the  satellite  intensities 
decreased  as  the  virtual  space  size  was  increased,  so  that  the  addition  of  Rydberg 
orbitals  to  relatively  small  basis  sets  such  as  the  Nakatsuji  basis  I  could  be  a  significant 
effect. 

To  clear  up  some  of  the  questions  raised,  the  study  was  repeated  using  the 
larger.  Dunning  ;?VTZ-|-  basis[5l],  which  includes  more  diffuse  s,p,d  and  f  func- 
tions.   The  pVTZ.  CCSD  optimized  geometry  is  Rcc  =  1.327A,  Rch  =  1.077A, 


Ta,ble  10.6:  Ionization  spectrum  of  ethylene,  10-30  eV 


Ionization  energy  eV  (Intensity) 

Experiment'' 


EOMIP 

pVTZ-l-  BASIS  I 


SAC-CP 


Orbital  Assignment 


10.72  (.95) 
13.14  (.93) 
14.97  (.92) 
16.37  (.88) 
19.65  (.85) 


10.39  (.95) 
12.85  (.93) 
14.59  (.92) 
16.03  (.88) 
19.36  (.85) 


24.52  (.007)  23.55  (.008) 
24.51  (.74)     24.25  (.74) 


26.13  (.015) 
26.16  (.052) 


26.73  (.052) 
26.81  (.006) 


10.25  (.95)' 
12.78  (.93)' 
14.50  (.92)' 
15.93  (.88)' 
19.32  (.85)' 
23.45  (.01) 
24.37  (.74)' 

25.86  (.04) 


26.97(.01) 


10.51  (.95) 
12.85  (.93) 
14.66  (.92) 
15.87  (.88) 
19.1   (.85) 


lb3u 
lb39 

lb2„ 
2biu 
2biu, 


23.7   (.74)         2ag 


2b^glb3ulag 


lb2u,  2b^jlb3glb3„ 
2aj,     7a^lb3„lb3u 
lb2u,  Sbjjlbsglbsu 

lb2u, 


2a, 


2h2g  Ib^g  lb3u 


■9' 


5at  lb 


lb. 


3u  i  U3u 


28.18  (.001) 
28.6  (.017) 
29.26  (.044) 

30.03  (.012) 


29.76  (.010)      29.28  (.02) 


30.06  (.094) 


30.77  (.075) 


30.15  (.09) 


27.4   (.29) 


31.2   (.1)*^ 


2biu,  3b|„  Ibsu  lb3u 
2a3,    7a^  Ib3ulb3u 
2biu,  2b2g  3a.g  Ibsu 
Ibsff,  2b2„  lb2u  lb3u 

J. 

Ihsg,  3h^g  lb3u  lb3u 

2a,g,    Ib^g  Ibsu  2biu 

2ag,    6aJ  Ibs^lbsu,  2b^g  2biu  Ibs^ 

2aig,    2h2g  2biu  Ibsu 


"  Nakatsuji  [54] 

'"  Banna  and  Shirley  [52],  Gelius  [57] 

'^  Theoretical  intensities  of  principcd  peaks  scaled  to  experimental  data 

**  Energy  and  intensity  from  Weigold  [58] 


EthylGne  Ionizations.  EOMIP,  Nakatsuji  basts 


Ethylene  lon(zation&.  EOMIP.  pVTZ+  basis 


J"i       i       i      it     lb     iV    lfe     it     *     *     21*     *     A 

Ethytone  Ionizations.  SAC-CI,  Nakatsuji  basis 


Figure  10.2:  EOMIP  and  SAC-CI  ionization  spectra  of  ethylene 


'^HCH  —  121.459°.  The  virtual  space  is  significantly  larger,  but  the  diffuse  functions 
less  so  than  those  in  the  Rydberg-supplemented  basis  used  by  Nakatsuji.  It  was  an- 
ticipated that  this  basis  set  would  represent  an  improvement  over  the  smaller  basis 
and,  since  ethylene  is  well  represented  at  the  CCSD  level  of  theory,  the  end  result 
would  be  improved  ionization  energies.  Instead,  what  is  observed  is  a  spectrum  of 
ionized  functions  which  are  generally  similar  to  those  obtained  in  the  smaller  basis 
(allowing  for  the  presence  of  diffuse  orbitals  in  the  low-ljdng  virtual  space),  but  with 
ionization  energies  which  grow  progressively  larger  than  the  experimental  energies 
and  those  calculated  in  the  smaller  basis. 

A  small  study  of  the  effects  of  basis  set  choice  and  orbital  relaxation  on  the  Sa^ 
ionization  (table  10.7)  reveals  that  removing  the  diffuse  functions  or  replacing  them 
with  Nakatsuji's  molecule-centered  Rydberg  set  has  little  discernible  effect  on  the 
ionization  energies  calculated  by  the  CCSD,  CCSD(T),  or  EOMIP  methods.  It  is 
apparent,  however,  that  orbital  relaxation  is  an  important  effect.  The  ACC  (QRHF 
calculations  for  the  cation)  and  EOMIP  results  in  the  C2H4  orbitals  over-estimate 
the  energy  by  about  .19  to  .34  eV,  while  in  the  C2H^  ROHF  orbitals,  ACC  (QRHF 
calculations  for  the  neutral)  methods  actually  underestimate  the  energy  by  over  1  eV, 
and  the  EOMIP  overestimates  it  by  a  similar  amount.  In  the  more  natural  calculation 
-  both  neutral  and  cation  correlated  calculations  are  based  on  the  appropriate  SCF 
calculation  -  both  of  the  CO  methods  perform  well,  within  .2  eV  of  experiment  in 
either  basis  set.  The  problem  for  the  CC  calculations  is  apparently  in  large  part  one 
of  orbital  relaxation.  Neither  set  of  orbitals  provides  a  reasonable  starting  point  for 
both  the  neutraJ  and  the  cation.  The  coupled-cluster  triples  correction  can  hardly  be 


said  to  be  much  help  in  improving  either  of  the  semi-QRHF  calculations.  In  the  ;)VTZ 
closed-shell  orbitals,  triples  lower  the  energy  of  the  neutral  by  .42  and  the  cation  by 
.38  eV.  In  the  pVTZ  cation  orbitals,  triples  lower  the  neutral  by  .47  and  the  cation 
by  .44  eV. 

Table  10.7:  Effect  of  orbital  choice  on  ethylene  Sa^  ionization 


basis 

MOs 

Energy 

eV 

ACCSD 

ACCSD(T) 

EOMIP 

pVTZ 

C2H4 

14.846 

14.877 

14.935 

C2H4 

13.458 

13.486 

15.902 

MIXED 

14.879 

14.855 

pYTZ 

C2H4 

14.847 

14.880 

14.936 

-l-Rydberg 

C2H4 

13.669 

13.752 

15.886 

MIXED 

14.881 

14.858 

Basis  I 

C2H4 

14.610 

14.627 

14.686 

C2HI 

13.491 

13.533 

15.660 

MIXED 

14.640 

14.604 

°ROHF  performed  with  3a.g  orbital  half  filled. 
''Neutral  in  Neutral  orbitcds,  cation  in  ROHF  orbitals 


The  difference  between  the  ACC  and  EOMIP  calculations  in  either  set  of  orbitals 
are  largely  attributable  to  spectator  triples  -  excitations  which  are  CC  doubles  but 
EOMIP  triple  excitations.  In  the  CCSD  calculation  on  the  neutral  in  these  orbitals, 
there  is  a  large  T2  ampUtude  (-.14272)  describing  the  excitation  of  both  Ibiu  (t^cc) 
HOMO  electrons  to  a  I62U  (ttcc)  orbital,  followed  by  another  double  excitation  n  to 
TT* (-.05468).  The  EOMIP  wavefunction  is  dominated  by  the  annihilation  of  the  Sag 
electron,  with  other  contributions  from  the  tt  ^  tt*  double  excitation,  followed  by  a 
single  excitation  from  the  I63U  HOMO  to  the  2b3u  virtual  orbital  and  one  from  the 


Ib2u  to  2b2u,  giving  a  wavefunction  dominated  by  the  action  of 

-.143  2blg2blglb3^1bl^3ag  +  .0953  2bljb3u3ag  +  .0839  2bljb2u3ag. 

For  the  cation  in  these  orbitals,  the  largest  amphtude  is  again  the  T2  providing  the 
double  excitation  from  the  163^  ttcc  bonding  HOMO  to  the  263„  ttcc  anti  bonding 
orbital,  followed  by  a  double  excitation  from  the  bonding  to  the  anti  bonding  orbital, 
combined  with  an  excitation  from  the  half-full  Sag  orbital  to  the  empty  4ag 

(-.127  2blg2b!,glbs,lbl  -  .0701  3al2biu2blglh,  -  .0547  2b\g•^\g\b^^\V^^)^ag 

It  is  the  second  term  here  that  is  missing  from  the  EOMIP  wavefunction  expression 
above,  where  it  contributes  as  a  T2  *  Ci  term,  and  has  an  amplitude  of  less  than  .001 
relative  to  the  principal  ionization.  EOMIP-SD  has  no  way  to  recover  this  determi- 
nant. Strangely  enough,  the  same  phenomenon  is  observed  in  the  C2H4  orbitals,  but 
the  effect  on  the  energy  is  not  nearly  as  large.  The  ACC  results  are  much  closer  to 
the  "mixed"  results  (around  .2  eV)  and  the  EOMIP  is  closer  as  well  (within  .1  eV  of 
the  ACCSD,  and  within  .3  eV  of  either  "mixed'  result). 

In  Nakatsuji's  Basis  I,  the  very  same  kind  of  behavior  is  again  observed.  In  the 
closed-shell  orbitaJs,  the  leirgest  contributors  to  EOMIP  wavefunction  are 

3a5+.098  3feJ„l63u3ap+.0531  462^  163^36^5 163^309 --0525  36^gl63u3a<,-.0489  \b\^\b2u'io-g, 

while  the  largest  contributors  to  02114"  CCSD  wavefunction  are 

30^  +  .143  '<^2g^^u'^b\g\b^^'Sag  -  .103  36^1^3^30^  -  .078  2b\jLb^^-Sag  -  .076  \b\gZag2h^^. 

The  first  three  determinants  are  the  same  in  either  case,  but  with  very  different 
relative  weights.  Still,  the  EOMIP  energy  is  very  slightly  closer  to  the  ACCSD  energy 
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than  is  the  case  in  the  pVTZ  basis,  with  or  without  Rydberg  orbitals.  Going  from  the 
small  BASIS  I  to  the  larger  pVTZ  basis  lowers  the  energy  of  the  closed-shell  ethylene 
more  than  it  does  the  cation,  which  means  that  all  the  calculated  ionization  energies 
are  lower  in  the  Nakatsuji  basis.  The  ACC  energies  in  the  closed-shell  orbitals  are 
actually  below  the  experimental  value  of  14.66  eV,  and  the  EOMIP  value  is  a  very 
good  match  to  experiment. 

Because  Neikatsuji's  basis  I  lowers  C2H4  ionization  energies  relative  to  more  com- 
plete basis  sets  it  manages,  without  improving  the  agreement  between  EOMIP  and 
ACC  methods,  to  provide  a  very  good  match  between  SAC-CI  or  CC  EOMIP  calcu- 
lations and  the  experiment. 

1Q,2,3    Butadiene 

As  the  experimental  spectrum  in  figure  10.3  indicates,  there  is  noticeably  more 
complexity  in  the  butadiene  than  in  the  ethylene  valence  spectrum.  The  broadening  of 
the  peaks  in  the  region  above  20  eV  indicates  the  possibility  of  shakeup  peaks  as  well 
as  the  onset  of  shake-off  effects  -  double  ionizations.  Because  the  PS  (photoemission 
spectroscopy)  experiment  was  performed  in  the  x-ray  region,  at  1487  eV,  the  experi- 
mental intensities  are  not  directly  comparable  to  those  obtained  using  the  formula  in 
equation  8.18.  However,  there  appears  to  be  good  correlation  between  the  positions 
of  experimental  and  theoretical  peaks.  The  tendency  of  higher  energy  peaks  to  be 
increasingly  higher  in  energy  than  experiment  is  consistent  with  the  trends  observed 
at  this  level  of  theory  for  large  basis  sets  with  several  polarization  functions. 


Butadiene  Ionizations,  EOMIP 


^  20  ^  ^4  3b" 


28 


Figure  10.3:  Ionization  spectrum  of  1,3  ^ra7i5- butadiene 
Simulated  experimental  spectrum  from  1487  eV  XPS,  reference  [59].    Theoretical 
peaks  are  convoluted  with  Lorentzian  functions  with  FWHM  of  leV. 


Table  10.8:  Ionization  of  1,3  ^rans-butadiene,  DZP+sp  basis. 
eV 


Energies  reported  in 


EOMIP 


Expt 


Orb 


EOMIP 


Expt 


Orb 


8.968(0.946) 

9.09(0.06) 

Ibg 

11.552(0.917) 

11.49(0.07) 

lau 

12.320(0.925) 

12.19(0.09) 

7ag 

13.621  (0.925) 
13.675  (0.920) 

13.29(0.11) 

13.76(0.13 

14.29(0.07) 

}bu 
Gag 

15.573(0.897) 
15.786(0.904) 

15.59(0.17) 

Sag 
5au 

16.651  (0.298) 

16.29(0.08) 

lau 

18.396(0.873) 

18.00(0.33) 

4bu 

18.854(0.007) 

Ibg 

19.724(0.839) 

19.19(0.64) 

4ag 

19.816(0.005) 

4bu 

20.903(0.007) 

2ag 

21.479(0.003) 

4bu 

21.796(0.003) 

Sag 

22.172(0.003) 

4bu 

22.483(0.047) 

3bu 

22.892(0.004) 

Ibg 

22.974(0.009) 

3ag 

23.052(0.161) 

3bu 

23.191(0.004) 

3bu 

23.191(0.002) 

3ag 

23.335(0.567) 

22.44(0.93) 

3bu 

23.374(0.016) 

Sag 

23.537(0.011) 

3bu 

23.678(0.006) 

3ag 

23.807(0.007) 

3bu 

22.815(0.012) 

3bu 

23.808(0.007) 

3bu 

23.888(0.005) 

4ag 

24.366(0.020) 

Sbu 

24.367(0.002) 

lau 

24.640(0.012) 

2bu 

24.687(0.009) 

3ag 

24.696(0.001) 

lau 

24.944(0.003) 

3bu 

24.991(0.419)     24.53 

(1.00)     3ag 

25.061(0.006) 

3bu 

25.194(0.026) 

3bu 

25.447(0.019) 

3ag 

25.501(0.007) 

lau 

25.527(0.002) 

3ag 

25.660(0.001) 

Sag 

25.728(0.008) 

Sag 

28.29(0.13) 

293.332(0.809) 

2ag 

293.332(0.809) 

2bu 

293.707(0.816) 

Ibu 

293.731(0.816) 

lag 

One  interesting  feature  of  this  conjugated  molecule  is  the  relatively  large  shakeup 
peak  which  appears  among  strong  principal  ionization  at  16.65  eV.  This  does  not 
occur  in  the  case  of  ethylene  (section  10.2.2),  where  the  first  strong  shakeup  appears 
slightly  higher  in  energy  than  all  the  valence  principal  peaks.  In  addition,  the  ethylene 
peak  shakeup  is  a  lb2u  satellite  of  a  peak  more  than  10  eV  lower  in  energy,  and  has 
an  intensity  less  than  6%  that  of  the  principal  peak.  The  butadiene  lau  sheikeup, 
however,  appears  only  5  eV  above,  and  has  an  intensity  nearly  a  third  that  of  the 
principal.  Lisini,  Pronzoni,  and  DeCleva  [60]  do  not  find  this  strong  shaikeup  using 
the  2hlp  CI  method. 

There  are  a  couple  of  noticeable  discrepancies  between  the  experiment  and  theory. 
Keane  et  al.  find  3  peaks  in  the  structure  centered  around  14  eV,  and  two  in  the 
structure  near  16  eV.  The  EOMIP  calculation  gave  only  two  peaks  near  14  eV,  and 
3  near  16  eV.  Of  the  three  near  16  eV,  the  EOMIP  peak  at  16.65  appears  to  be 
assignable  to  the  experimental  16.29  eV  peak,  leaving  the  pair  at  15.57  and  15.79  eV 
to  be  assigned  to  the  15.59  eV  experimental  peak.  The  nearly-degenerate  pEiir  at  about 
13.6  eV  can  be  assigned  to  the  experimental  structure  centered  around  the  13.76  eV 
peak.  These  are  all  strong  principal  ionizations,  as  evidenced  by  the  pole  strengths 
given  in  table  10.8,  and  it  is  difficult  to  conceive  of  the  energies  being  improperly 
ordered  at  this  level  of  theory.  The  energy  is  too  low  for  double  ionization  processes 
to  play  an  important  role.  The  only  remaining  possibility  is  that  of  relatively  strong 
shakeup  peak  in  near  14  eV.  Still,  this  would  not  account  for  the  "extra"  peak  found 
in  the  EOMIP  calculation  near  16  eV.  A  lower-energy,  UPS  (ultraviolet  photoemission 
spectroscopy)  experiment  or  an  electron  momentum  experiment  would  give  different 


relative  intensities  between  states  of  different  symmetries,  yielding  more  information 
about  the  composition  of  these  two  structures. 

10.3    Saturated  Molecules 
10.3.1     Ammonia 


ACESII  calculations  on  NH3  are  carried  out  in  the  Cj  symmetry  point  group,  the 
largest  Abehan  subgroup  of  the  full  Cat,  point  group  to  which  NH3  belongs.  The 
C3,;  electronic  configuration  is  written  as  laj2ajle*3ai,  while  the  computational  Cj 
configuration  is  la'  2a'  la"  3a'  4a'  ,  where  the  degenerate  pair  (2a',  la")  is  equivalent 
to  the  le  orbital  in  the  C31,  representation. 

Table  10.9:  NH3  EOMIP  error  vs.  experiment 


Expt.     EOMIP        Expt.    EOMIP 
10.85  0.05  27.6  0.30 

16.50  0.16 


The  calculated  outer-valence  EOMIP  energies  (10.94  and  16.62  eV)  are  in  good 
agreement  with  experimental  results  (10.85  and  16.5  eV,  references  [61,  62,  63]).  The 
experimental  peak  at  27.6  eV  is  matched  to  .3  eV  by  the  EOMIP  result  of  27.91 
eV,  and  the  shakeup  roots  on  the  high-energy  end  of  this  peak  are  in  qualitative 
agreement  with  the  "tail"  on  the  experimental  peak.  The  researchers  in  reference  [61] 
fit  the  shape  of  this  undulating  tail  to  4  wide  Lorentzian  hneshapes  (see  table  10.10). 


NH3  Ionizations.  EOMIP  and  approx  Expt. 
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Figure  10.4:   NH3  Ionization  Spectrum.    Lorentzian  lineshapes  applied  to  reported 
experimental  energies  and  intensities. 

In  the  current  study,  the  structure  found  in  this  region  was  qualitatively  similar  in 
appearance  but  smaller  in  intensity. 

The  trend  of  higher-energy  ionizations  to  be  progressively  further  from  the  exper- 
imental values  is  consistent  with  results  for  ethylene  and  N2  (section  10.2). 


10.3.2    Water 


Compared  to  ammonia,  water  shows  a  much  richer  inner-valence  structure,  with 
several  shakeup  peaks  surrounding  the  inner  valence  peak  structure  between  32  and 
35  eV.  The  EOMIP  calculation  shows  good  agreement  with  experiment  in  the  valence 
region,  including  the  2ai  shakeup  peaks,  but  does  not  reproduce  the  amount  of  in- 
tensity in  the  weak  satellites  near  27  and  38  eV.  Bawagan  and  Davidson,  in  reference 
[64],  found  more  shakeup  intensity,  using  a  109  Gaussian  function  basis  and  a  2hlp  CI 


Table  10.10:  NH3,  pVTZ+  basis 


Energ\'  eV 

(Intensit)^) 

EOMIP 

Experiment" 

Orbital  Assignment 

10.902  (.941) 

10.85  (.956) 

4a' 

16.616(1.56) 

16.5  (1.58) 

le 

27.907  (.81.3) 

27.6  (.813) 

2a' 

29.709  (.009) 

30.3  (.13) 

2a' 

8a'^3a'3a' 

30.628  (.0005) 

2a' 

10a'+3a'3a' 

32.361  (.003) 

33.2(.ll) 

la' 

7a'+3a'3a' 

33.453(.016) 

2a' 

7a'+3a'3a', 

33.639  (.0005) 

la" 

8a'Ua"4a' 

35.252  (.0002) 

la" 

3a"Ua'4a' 

34.631  (.0002) 

la" 

3a"+4a'4a' 

35.074  (.0004) 

35.5(.ll) 

la" 

(10a'Ua"4a'  -  3a"t3a'4a') 

37.595  (.0008) 

la" 

8a/Ua"3a' 

37.939  (.0003) 

41.8  (.08) 

la" 

(3a"^4a'3a'  +  9a't4a'la") 

"  Davidson  et  al.  [61] 


Table  10.11:  H2O  EOMIP  error  vs.  experiment 
Expt.     EOMIP         Expt.     EOMIP 


12.61 
14.73 
18.55 


0.00 
0.13 
0.38 


32.30 

34.90 

539.70 


0.46 
0.19 
1.36 


Table  10.12:  Ionization  of  H2O 


Ionization  energy  eV 

(Intensity) 

Koopmans' 
theorem 

EOMIP 

Exp 

Orbital 
AssiKnment 

13.884 

12.609 

f.940) 

12.61 

(1.0) 

I61 

15.935 

14.858 

.940 

14.73 

0.8 

3ai 

19.432 

18.925 

(.950) 

18.55 
27.0 

(0.6) 
(0.02) 

I62 

36.786 

32.760 

[•707) 

32.30 

(0.44) 

2ai,  5ah62l62 

— 

35.092 

[.160) 

34.90 

(0.15) 

2ai,  5a|l62l62 

— 

37.828 

;.oo6) 

37.8 

(0.078) 

2ai,  Sajloilai 

— 

40.421 

[.007) 

2ai,  2623ail62 

— 

40.609 

;.ooi) 

8ajl62l62 
269162301 
76,301 1 61,  2a  1 
26h62l62,  2oi 
4a|l62l62,  3oi 
8ah62l62,  2ai  +  3ai 

— 

41.524 

,.001) 

— 

41.759 

;.oo2) 

— 

44.758 

,.001) 

— 

45.117 

,.0007) 

— 

45.496 

,.001) 

— 

46.348 

,.003) 

7a|3ai3oi,  2ai  +  3oi 

— 

46.749 

,.001) 

7aJ3ai3ai,  3ai 

— 

47.077 

.012) 

5aj.3ai3oi,  3ai 
36J3ail6i,  3ai 
36|l6i3ai,  2ai 
9a|l62l62,  2ai 

— 

47.256 

.002) 

— 

48.713 

.003) 

— 

48.728 

.002) 

— 

49.243 

.0004) 

8aj3ailoi,  2oi  -  3ai 
26j3ail6i,  3oi 

— 

50.069 

.001) 

— 

50.119 

.003) 

10aj362362,  2ai 

— 

50.936     ( 

.0004) 

5a]36i36i 
46^301 162 

— 

51.077     ( 

.0004) 

559.638 

541.062     ( 

.830) 

539.7 

lai,  42ol2ai2ai 

approach.  However,  they  report  energies  of  33.1  and  33.6  eV  for  the  2ai  ionizations, 
an  average  of  about  .6  eV  in  error  compared  to  the  average  of  .3  eV  in  error  for  the 
EOMIP  calculation  on  the  same  two  peaks. 

HjO  Ionizations,  EOMIP 


Expt 

EOMIP    


~X         it       ife       1^4       i¥       i¥       ^       S       ^4       S      S       ^       3?       3^       ^       35       ^b 

Ionization  Energy  (eV) 

Figure  10.5:    H2O  EOMIP  ionization  spectrum.    Lorentzian  lineshapes  apphed  to 
reported  experimental  energies  and  intensities. 


10.4     Nitrogen  and  C?:  A  study  in  contrasts 


10.4.1     Single-Determinant  and  Multideterminant  References 


The  C2X^E^  —^  C2"a^nu  transition  is  a  difficult  test  for  any  method  because  both 
species  have  significant  multireference  character.  The  C2  X'  S^  ground  state  is  char- 
acterized by  two  electronic  configurations,  Ia'^lal2a^2all7i^  and  laglal2a'^3a'^li:^. 
The  CJ  state  is  dominated  by  the  lcr^lcr^2cr|2(7„l7r^  configuration,  with  a  major  con- 
tribution from  I(7jlc7^2(7g3c7gl7r;^,  and  smaller  contributions  from  Iajla^2a^2aul7ril7rj 
and  other  electronic  configurations.  The  N2X^EJ"  -^  N^X^Eg  problem,  by  contrast,  is 
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computationally  much  less  demanding.  The  N2  neutral  state  is  well  characterized  by 
only  one  electronic  configuration,  and  the  N^X^E^  ionized  state  has  a  much  smaller 
multi-reference  character  than  Cja^n„,  with  l(Tglcr^2cr^2ai3cr^l7r^l7r]  weighted  only 
one  tenth  of  the  major  configuration  in  the  wavefunction.  A  comparison  of  the  CCSD, 
CCSD(T),  and  EOMIP  energies  and  wavefunctions  of  these  two  molecules  near  their 
minima  should  provide  some  insight  into  the  nature  of  the  EOMIP  solution  at  these 
two  extremes.  A  pVTZ+sp  basis  was  chosen  for  all  calculations,  and  the  optimum 
geometries  of  the  molecules  and  the  studied  cations  were  determined  in  that  basis  at 
the  CCSD  level  of  theory. 

The  EOMIP  energy  calculated  for  NJ  differed  from  the  CCSD  energy  by  an  aver- 
age of  .011  eV  over  the  range  shown  in  table  10.13,  with  a  maximum  of  .03  eV  at  the 
longest  bond  length,  1.150  A.  The  difference  between  EOMIP  and  CCSD  energies  of 
C2"  over  the  range  1.180  to  1.360  A,  shown  in  table  10.14,  differ  by  an  average  of  .025 
eV,  over  twice  as  large.  The  EOMIP  and  ACCSD(T)  electron  detachment  energies 
for  N2  differ  by  an  average  of  .12  eV,  while  for  C2  these  differ  by  an  average  of  .55 
eV. 

For  both  molecules,  EOMIP  state  energies  (CCSD  energy  of  the  neutral  -h  EOMIP 
electron  detachment  energy)  were  closer  to  the  CCSD  than  to  CCSD(T).  In  figure 
10.6a,  the  Nj  EOMIP  and  CCSD  curves  appear  nearly  superimposed,  about  .02  au 
(.5  eV)  removed  from  the  CCSD(T)  curve.  In  figure  10.6b,  the  EOMIP  and  ACCSD 
energies  again  track  closely,  differing  from  the  ACCSD(T)  energies  by  about  .01  eV. 
Although  the  EOMIP  electron  detachment  energies  for  C2  are  closer  to  the  ACCSD 
(figure  lO.Gd),  than  to  the  ACCSD(T)  energies,  they  cannot  be  said  to  be  particularly 


close  to  either,  differing  by  about  .5  eV  from  the  ACCSD  energies,  which  in  turn  were 
separated  from  the  ACCSD(T)  energies  by  about  .15  eV. 

As  a  general  rule  of  thumb,  CCSD  underestimates  the  stabihty  of  multiconfigura- 
tion  references,  while  CCSD(T)  might  overestimate  it,  relative  to  Full  CI  (FCI)[65]. 
For  the  simple  case  of  N2,  where  the  reference  state  is  very  much  a  single  configu- 

Table  10.13:  N2  X^Eg  and  Nj  X^Sg  in  pVTZ+sp  basis 


N2  XI E+ 

N+  X2S+ 

Ionization 

R-NN 

CCSD 

CCSD(T) 

EOMIP 

CCSD 

CCSD(T) 

EOMIP 

CCSD 

CCSD(T) 

A 

-109. 

(au) 

-108.  (au) 

(eV) 

1.040 

.398  874 

.416  579 

.818  616 

.817  765 

.839  868 

15.790 

15.813 

15.693 

1.060 

.404  342 

.422  754 

.825  668 

.825  069 

.847  923 

15.747 

15.763 

15.642 

1.080 

.406  948 

.426  097 

.829  808 

.829  505 

.853  123 

15.705 

15.713 

15.592 

1.090 

.407  296 

.426  825 

.830  905 

.830  766 

.854  771 

15.685 

15.723 

15.566 

1.091" 

.407  299 

.426  862 

.830  974 

.830  851 

.854  891 

15.683 

15.686 

15.564 

1.100 

.407  066 

.426  982 

.831  411 

.831  449 

.855  844 

15.665 

15.663 

15.542 

1.107* 

.406  511 

.426  736 

.831  422 

.831  605 

.856  307 

15.649 

15.644 

15.522 

1.110 

.406  297 

.426  610 

.831  367 

.831  593 

.856  382 

15.645 

15.639 

15.517 

1.130 

.403  295 

.424  423 

.829  781 

.830  422 

.856  006 

15.606 

15.590 

15.468 

1.150 

.398  570 

.420  545 

.826  428 

.827  536 

.853  928 

15.569 

15.539 

15.419 

AEA 

15.67r 

15.666 

15.526 

"VEDE 

'■VEA 

^Calculated  as  CCSD(N2,  1.107)  -  CCSD(N2,  1.091)  +  EOMIP(1.107) 


ration  function  and  the  ionized  state  is  very  slightly  multiconfigurational,  it  is  easy 
to  conclude  that  the  ACCSD  (T)  energy  will  be  lowered  relative  to  the  ACCSD  en- 
ergy. And  this  is  the  case.  The  energy  of  the  cation  is  generally  expected  to  be  more 
stabilized  than  the  energy  of  the  neutral  by  higher  excitations  (such  as  the  triples 
approximation),  lowering  the  ionization  energ>-.  The  EOMIP  energies  closely  follow 
the  ACCSD  energies,  which  seems  reasonable. 

The  C2  case  is  different,  however.  For  C2  the  triples  are  more  important  in  the 
neutral  than  in  the  ionized  state.  The  CCSD  and  CCSD(T)  energies  of  the  CsX^EJ- 


Table  10.14:  C2  X^EJ"  and  C^  a^n„  in  pVTZ+sp  basis 


C2  XI E+ 

c+  a2n„ 

Ionization 

Rcc 

A 

CCSD 

-75. 

CCSD(T) 
(au) 

EOMIP 

CCSD 

-75.  (au) 

CCSD(T) 

EOMIP 

CCSD    CCSD(T) 

(eV) 

1.180 

.767  465 

.801  316 

.279  333 

.304  090 

.333  627 

13.283 

12.609    12.727 

1.200 

.771  083 

.805  085 

.286  498 

.311  277 

.340  748 

13.186 

12.512    12.635 

1.220 

.773  176 

.807  351 

.291  977 

.316  807 

.346  242 

13.094 

12.419    12.548 

1.240 

.773  918 

.808  287 

.295  954 

.320  863 

.350  286 

13.006 

12.328    12.463 

1.242 

°       .773  924 

.808  313 

.296  268 

.321  185 

.350  608 

12.998 

12.320    12.455 

1.260 

.773  467 

.808  048 

.298  598 

.323  610 

.353  041 

12.922 

12.200    12.347 

1.280 

.771  965 

.806  774 

.300  061 

.325  199 

.354  657 

12.841 

12.157    12.303 

1.300 

.769  544 

.804  595 

.300  481 

.325  766 

.355  265 

12.764 

12.076    12.227 

1.302 

''       .769  220 

.804  298 

.300  469 

.325  772 

.355  276 

12.755 

12.067    12.219 

1.320 

.766  319 

.801  624 

.299  984 

.325  437 

.354  988 

12.690 

11.997    12.154 

1.340 

.762  399 

.797  969 

.298  684 

.324  322 

.353  936 

12.618 

11.921    12.083 

1.360 

.757  878 

.793  723 

.296  683 

.322  524 

.352  209 

12.550 

11.847    12.014 

AEA 

12.883^ 

12.187    12.328 

°VEDE 

*VEA 

^Calculated  as  CCSD(C2,  1-302) 
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Figure  10.6:  N2  and  C2  potenriaJs  and  ionization 


neutral,  are  further  apart  than  those  for  N2,  but  those  for  the  C2"a^nu  state,  which  has 
even  more  multiconfiguration  character,  are  closer  together.  Thus,  the  ACCSD(T)  is 
greater  than  ACCSD.  And  the  EOMIP  is  now  greater  than  ACCSD(T).  It  is  apparent 
that  this  example  defeats  simple  rules  of  thumb.  The  EOMIP  energy  is  in  fact  worse 
than  the  Koopmans'  approximation  of  12.45  eV. 

For  the  N2  ionization,  the  experimental  vertical  electron-detachment  energy  (VEDE) 
is  15.5  eV.  The  CCSD(T)  VEDE  energy  is  in  error  by  less  than  a  tenth  of  an  eV, 
while  the  CCSD(T)  AEA  result  is  in  error  by  about  a  quarter  of  that.  Any  of  the 
EOMIP  results  differ  by  between  .1  and  .2  eV,  still  within  "experimental"  accuracy, 
and  accurate  enough  to  interpret  with  confidence  the  doublet  ionization  spectrum  of 
N2XIE+  (table  10.3). 

For  the  C2  ionization,  ionization  peaks  have  not  been  experimentally  assigned, 
but  according  to  Watts'  best  calculations[66],  the  EOMIP  calculation  is  not  accurate 
enough  to  distinguish  between  the  C2X'S+  -)•  Cj^a^IIu  (at  12.51  eV)  and  C2X^E+  -^ 
C2'X'*i;^  (at  12.60  eV)  transitions.  The  EOMIP  calculation  of  the  former  weighs  in  at 
12.998  eV,  hardly  within  "experimental  error"  of  either  peak.  The  latter  transition, 
of  course,  is  fundamentally  2hlp,  which  means  there  is  no  relaxation  at  all  to  be 
provided  by  a  EOMIP-SD  calculation.  The  test  calculation  provided  a  VEDE  of  over 
15  eV  (table  11.1)-  not  useful  to  this  task  (there  is,  however,  another  approach  to 
the  EOMIP  calculation  of  the  energy,  which  will  be  presented  in  section  11.1). 


10.4.2    The  Spectator  Triples  Problem 


At  the  equihbrium  distance  for  the  neutral  species,  the  most  important  config- 
urations in  the  N2"  CCSD  and  EOMIP  wavefunction  are  lcrgla„2cr^2cr„3a^l7r^  and 
lallal2al2al3all7TllTrl,  with  coefficients  in  the  ratio  1:-0.1  for  CCSD,  and  l:-.077 
for  EOMIP: 


There  is  qualitative  agreement  in  the  general  form  of  the  wavefunction.  In  the  CCSD 
calculation,  the  second  determinant  is  created  by  the  T2  element  associated  with 
^ixu27'  operating  on  the  cation  reference  determinant  3a g  N2X^E^y  while  in  the 
EOMIP  calculation  it  is  created  by  the  action  of  the  Citt^^^I  operator  -  double  exci- 
tations for  both  methods. 

For  the  C2^  CCSD  wavefunction,  the  major  components  are 

C^^a^n^)    =    {l  +  tlZ%[3al2;al2a,2a, 

+  t;::Sl7rJl7fJl7f„l7r„  +  4;l:n7fi3al2a,mg 
+  tl'/Japau  +  g:3aj2a„)  l7f„  |C2X^E+) 

with  the  determinants  in  the  ratios  1:-. 351:-. 106:. 098:-. 083:-. 081,  resulting  from  the 
action  of  the  wave  reaction  operator  {1  +  Ti  +T2  +  •  •  ■)  on  the  N-1  electron  reference 
l7f„  02X^2^  y  In  the  EOMIP  wavefunction,  these  determinants  appear  in  the  ratio 
l:-.278: 


(<  .06):— .078:— .119:— .046.  In  addition,  there  are  three  more  determinants: 

(.081  l7ril7r„l¥„  +  .075  3^il7r„l7f„  +  .072  37f|,l7f„l7f„)  |C2X^E+) 

which  are  at  least  as  important  to  the  EOMIP  wavefunction  as  some  of  the  deter- 
minants above,  but  which  are  not  as  important  to  the  €2"  CCSD  wavefunction.  The 
most  striking  difference  between  the  CCSD  and  EOMIP  wavefunctions  is  the  coeffi- 
cient to  the  iTrhTfUTfulTTu  determinant.  It  appears  in  the  EOMIP  wavefunction  as 
a  product  of  the  C2  CCSD  tj^^lj  coefficient  (|ampntude|  <.06)  and  the  EOMIP  t^^ 
coefficient.  The  difference  between  the  sizes  of  this  t  amplitude  in  the  two  CCSD 
calculations  is  quite  large,  and  could  only  be  made  up  by  the  EOMIP  calculation  if 
it  included  the  ii¥„i^'i^*  operator,  a  triple  excitation  in  EOMIP  parlance.  In  FSIP 
parlance,  triples  of  the  form  ^j^^j^^j^^  are  commonly  referred  to  as  spectator  triples 
because  an  active  hole  is  both  annihilated  and  created,  and  therefore  appears  as  a 
sort  of  "spectator"  line  to  a  double  excitation  process.  As  discussed  in  chapter  4.2, 
there  is  a  simple  correspondence  between  the  FS  S  and  the  EOMIP  t  amphtudes 
which,  in  cases  like  this  where  only  one  of  the  EOMIP  Tj  amplitudes  is  important,  is 
nearly  identity.  The  desired  iij„i^*ij'  is  therefore  the  EOMIP  equivalent  of  the  FSIP 
spectator  triple  excitation.  A  perturbation  theory  analysis  of  the  EOMIP  equations 
indicates  that  EOMIP  triples  contributions,  which  appear  at  third  order  in  W  (ap- 
pendix B)  are  possibly  more  important  in  cases  like  this  than  triples  corrections  to 
the  CCSD  energ}'.  This  is  borne  out  by  the  importance  of  the  spectator  triples  to  the 
description  of  the  C  J  wavefunction.  It  is  clear  that  the  absence  of  the  spectator  triples 
is  principally  responsible  in  this  case  for  the  inabihty  of  the  EOMIP  wavefunction  to 
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mimic  qualitatively  the  €2^  CCSD  wavefunction,  and  the  author  asserts  that  the  at- 
tempt to  reproduce  the  proper  wavefunction  in  the  absence  of  this  important  part  of 
the  manifold  space  also  results  in  the  relatively  large  weights  of  the  three  odd  EOMIP 
determinants  which  are  not  important  to  the  02"  CCSD  wavefunction.  Attempts  to 
improve  FSIP  energies  by  adding  triples  approximations  have  had  mixed  results[40], 
but  it  would  be  an  interesting  course  of  extended  further  study  to  examine  the  effects 
of  various  levels  of  approximation  to  the  addition  of  spectator  triples  effects  to  the 
wavefunction. 

For  the  single-reference  N2  case,  of  course,  the  spectator  triples  do  not  make  an 
important  contribution,  and  the  EOMIP  calculation  does  a  good  job  of  reproducing 
the  CCSD  result. 


CHAPTER  11 
LINEAR  CARBON  CLUSTER  ANION  PHOTOELECTRON  SPECTRA 

Carbon  clusters  have  been  found  in  interstellar  space,  where  they  are  suspected  of 
being  important  in  the  formation  of  man}^  carbon-containing  compounds  [67,  68],  and 
in  hydrocarbon  flames,  where  they  are  important  in  the  formation  of  soot  particles 
[69].  They  can  be  created  in  the  laboratory  by  super  heating  graphite,  separated  by 
mass  spectrometry,  and  the  smaller  clusters  have  been  studied  with  more  and  more 
exacting  techniques  in  recent  years  [70,  71,  35,  72,  36].  Moving  from  larger  to  larger 
clusters,  researchers  hope  to  observe  the  transition  from  cluster  to  extended  material. 
Yet  there  are  still  many  questions  about  the  structure  and  properties  of  small  clusters. 
State-of-the  art  quantum  theory,  employing  high-level  ab  imtio  methods,  permits  ac- 
curate predictions  of  the  properties  observed  in  spectroscopic  experiments.  Matching 
theoretical  calculations  of  vibrational,  excitation,  and  ionization  energies  to  experi- 
mentally obtained  spectra  can  and  has  provided  information  about  the  structure  of 
the  species  in  question.  Early  experimental  analyses  of  carbon  clusters  produced  in 
hot  carbon  vapors  [73]  and  among  molecular  dissociation  products  [74]  were  invari- 
ably comphcated  by  the  mixture  of  electronic  and  vibrational  states  produced  by 
such  methods.  Early  theoretical  approaches  based  on  MO  theory  or  low-level  correla- 
tion approaches  [75]  yielded  energies  roughly  in  agreement  with  the  experiments,  but 
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electronic  configuration  assignments  which  have  been  discounted  by  more  recent  ex- 
perimental and  theoretical  results.  Thirty  years  of  improved  results  from  increasingly 
sophisticated  experimental  and  theoretical  designs,  laid  out  in  Weltner  and  Van  Zee's 
1989  review  [70],  led  to  the  present  situation.  During  an  extended  debate  between 
experimentalists  and  theoreticians  [76,  77,  78,  79,  80],  a  picture  began  to  form  of 
the  structures  of  the  small  carbon  clusters.  Recent  photoionization  experiments,  in- 
volving the  separation  by  time-of-flight  mass  spectrometry  of  a  beam  of  mixed  cairbon 
cluster  anions,  have  produced  some  very  high  quality  spectra  of  the  small  carbon  clus- 
ters [71,  35,  72,  36].  These  have  been  followed  up  by  some  very  high-level  theoretical 
calculations  on  the  geometries  and  electronic  structures  of  the  neutral  and  anionic 
clusters  [66,  65,  81,  82,  83,  84],  in  which  the  geometries  and  ground  electronic  states  of 
the  neutral  carbon  clusters  C2  -  Cio  and  their  anions  have  been  calculated  very  accu- 
rately [65].  The  vibrationally  resolved  spectra  from  Neumark's  group  [72,  36,  85]  are 
highly  resolved,  reveaUng  vibrational  and  electronic  structure  of  the  transition  within 
a  window  of  2  to  3  eV  above  the  first  ionization  energy  of  each  cluster  anion,  C^  - 
C^,  with  which  the  results  of  Watts  [66]  are  in  very  good  agreement.  The  spectra  of 
C^  -  Cjo  produced  in  Small ey's  group,  while  much  less  resolved,  reveal  a  wealth  of 
structure  in  a  range  of  3  to  5  eV  above  the  first  ionization.  C^  in  particular  exhibits  a 
very  rich  structure,  and  the  other  small  linear  species  display  several  broad  features. 
Most  of  these  peaks  are  as  yet  unassigned;  many  must  be  due  to  transitions  to  excited 
states  of  the  neutral  clusters.  Considering  the  main  strength  of  the  EOMIP  method  - 
an  economical  and  conceptually  simple  approach  to  calculating  several  ionizations  at 
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once  -  the  author  has  applied  the  method  to  the  interpretation  of  the  photoelectron 
spectra  of  these  novel  molecules. 

Anions  of  linear  carbon  clusters  present  a  challenge  to  any  theoretician.  For  chains 
of  three  carbons  and  longer,  the  valence  orbitals  of  the  neutral  are  of  tt  symmetry. 
The  anions  involve  either  one  or  three  electrons  in  a  tt  orbital.  With  a  code  that  is 
limited  to  Abehan  symmetries,  there  are  plenty  of  opportunities  for  both  spin  con- 
tamination and  symmetry  breaking  at  the  Hartree-Fock  level.  C2  exhibits  significant 
multireference  character.  While  the  CCSDT  approach  can  be  expected  to  overcome 
these  defects  of  the  reference  determinant,  and  has  been  shown,  in  large  basis  sets,  to 
yield  good  agreement  with  experiment  [65,  66],  the  CCSD  method  cannot  be  expected 
to  do  so  well.  This  has  been  demonstrated  with  the  odd-numbered  neutral  chains, 
where  there  is  the  indication  that  non-iterative  triples  corrections  are  not  sufficient 
to  describe  the  wavefunction  [65];  some  type  of  iterative  triples  approach  is  indicated. 

Despite  these  difficulties  Watts  and  Bartlett  reported  ACCSD  calculations  within 
.2  eV  of  the  experiments  for  the  first  ionizations  of  the  linear  carbon  clusters  C^ 
through  Cg  ,  using  DZP+sp  or  DZP+p  bases  and  UHF  orbitals.  Since  EOMIP  re- 
sults can  come  very  close  to  reproducing  CCSD  energies  in  many  cases,  it  is  not 
unreasonable  to  apply  the  EOMIP  method  to  the  study  of  the  outer  valence  spectra 
of  these  carbon  clusters. 
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11.1     Co,  Cj-.  and  C^ 

Recent  coupled-cluster  studies  of  C2  and  its  anion  have  demonstrated  good  agree- 
ment with  experimental  results  when  pVTZ  and  pVQZ  basis  sets  were  used  at  the 
CCSD(T)  level  of  theory  and  higher [66,  65].  Here,  the  pVTZ  basis  set  is  used,  and 
geometry  is  optimized  at  the  CCSD  level  for  both  the  neutral  and  the  anion.  These 
CCSD  geometries  are  quite  close  to  that  obtained  at  the  CCSD(T)  level  although  the 
reference  energies  differ  somewhat. 

ILIA Ionization  of  Co 

A  detailed  study  of  the  C2  X'E^  and  C^  a^IIu  transition  in  the  pVTZ+sp  basis 
was  presented  in  table  10.14.  Here,  that  ionization  is  presented  again,  along  with 
some  2hlp  transitions,  in  table  11.1. 

These  vertical  detachment  energies  from  the  ground  state  of  C2  are  not  experimen- 
tally verified,  but  the  comparison  to  calculations  of  Watts  and  Bartlett[66]  indicate 
that  EOMIP  has  not  done  especially  well  in  this  very  difficult  case.  Especially  diffi- 
cult are  the  ^E^  and  ^E"  calculations.  These  fundamentally  2hlp  transitions  cannot 
be  very  well  represented  at  the  CC  EOMIP-SD  level  of  theory  except  in  the  most 
fortuitous  of  circumstances. 

The  C2  a^Hu  — >  C^  X'^^p  transition  should  be  much  more  amenable  to  EOMIP 
treatment.  The  neutral  is  well  characterized  by  the  single  ROHF  determinant 
Icr^lc7^2cr^2a„l7r^3(7gn  and  the  cation  by  Ia^l(j22cr^2(7^l7r„3cr^,  the  result  is  an  energy 


12.991 

12.998 

12.45 

15.366 

12.51 

16.961 

12.231 

11.67 

12.506 

12.23 

Table  11.1:  Ionization  of  02-  Energies  reported  in  eV 

EOMIP  CCSD(T) 

C2         C^     Reference pVTZ     pVTZ+^p     pVTZ° 

a^Ilu     ItTu 
'^T,-     3aU7ruil7r„j, 
^E-     3crjl7ruxl7ru3/ 

a^n^ 

X^S;     1¥„ 
a^riu     Sag,  'Sa\7a^2a^ 

"  Watts  and  Bartlett  [66] 

.56  eV  greater  than  the  CCSD(T)  result,  which  is  itself  perhaps  an  overestimation 
of  the  true  value  [66].  The  case  of  the  transition  to  a^IIu  state  of  the  cation  should 
be  a  more  difficult  one.  With  significant  contribution  from  the  a  doubly-excited 
determinant  (.27  in  lcr^lcT^2o-^l7r23cr^),  it  would  be  expected  that  the  EOMIP  energy, 
like  the  ACCSD,  would  overestimate  the  ionization  energy,  while  the  ACCSD(T) 
should  underestimate  it.  This  is  consistent  with  the  results,  but  here  the  EOMIP 
energy  is  only  about  .3  eV  from  the  ACCSD  energy,  apparently  better  performance 
than  that  obtained  from  the  seemingly  more  tractable  C2  a^IIu  — >■  C^  ^^^3  transition. 


11.1.2     Ionization  of  Ct 


The  most  accurate  experimental  [36,  71]  results  for  the  photodetachment  of  C^ 
indicate  that  the  EOMIP  value  of  the  X^E+  -^  X^E^  vertical  electron  detachment 
energy  falls  within  .3  eV  of  the  experimental  result  (see  table  11.2).  Previous,  highly 
accurate  theoretical  results  [66,  65]  indicate  that  better  agreement  with  experiment  is 


Table  11.2:  Ionization  of  C2  in  p\TZ->rsp  basis,  UHF  orbitals.  Energies  reported  in 
eV. 

C2     Reference  EOMIP         ACCSD(T)"     Experiment 


X^E^ 

UHF/CCSD 

-75.900280  au 

X1E+ 

^^t 

3.529     3.115 

3.273'' 

^n^" 

ItTu 

4.019 

'K' 

20g 

5.500 

15^+5 

H 

15.752 

lV+6 
^9 

H 

285.167 

a^n. 

!¥„ 

3.627 

3.35r 
?s4.3?'= 

'K 

^t 

4.761 

"  Watts  and  Bartlett  [66],  in  pVQZ  basis 
"  Arnold  et  al.  [36];  '  Yang  et  al.  [35] 

possible,  but  that  the  results  are  affected  both  by  increasing  the  basis  set  from  pVTZ 

to  pVQZ,  and  by  the  inclusion  of  triples  in  the  coupled-cluster  calculation  (both  cause 

changes  to  the  energy  on  the  order  of  .leV)  .  As  a  rule,  CCSD  underestimates  the 

stability  of  a  multi-configuration  reference  relative  to  the  FCI  result,  while  CCSD(T) 

overestimates  it.  Therefore,  it  is  not  surprising  that  the  EOMIP  result  (which  should, 

except  in  pathological  cases,  be  quite  close  to  the  ACCSD  result)  is  greater  than  the 

experimental  energy,  while  the  ACCSD (T)  result  is  lower.  Simulated  C^  spectra  are 

found  in  section  11.4 

11.2     Ionization  of  C^ 

Electron  detachment  energy  calculations  on  C^  were  performed  using  3  different 
sets  of  orbitals  -  orbitals  from  a  C3  UHF  calculation,  from  a  C3  RHF  calculation,  and 
from  a  Cl'  RHF  calculation.   While  the  UHF  calculation  is  expected  to  reproduce 


the  energy  of  the  anion  fairly  well,  the  calculation  on  a  non  symmetric  wavefunction 
produces  both  a  poor  approximation  to  the  true  spin  eigenfunction  -  the  average 
multiplicity  at  the  UHF  level  is  2.117  -  and  a  set  of  n  orbitals  in  which  the  x  and  y 
(Bsu  and  B2u  in  the  D2/1  Abelian  subgroup  in  which  the  computations  were  performed) 
components  differed  significantly.  Although  the  approximate  spin  multiplicity  is  2.092 
at  the  MBPT(2)  level  and  2.013  when  the  CCSD  calculation  converges,  it  might  be 
expected  that  the  UHF  orbitals  are  unsuited  to  represent  the  closed-shell  C3.  While 

Table  11.3:  Cg-X^n^  -^  C3XIE+  in  pVTZ-\-sp  basis. 

Method E(C3  )  Hartrees       VEDE  eV Source 

ACCSD  -113.915414534294  2.038 

EOMIP,  C3"  orbitals  -113.915414534294  2.157 

EOMIP,  C3"  orbitals  -113.915428019063  -1.133 

EOMIP,  C3  orbitals  -113.913530561817  1.989 

ACCSD(T)  1.61  Watts  and  Bartlett  [65] 

Experiment  1.98  ±  0.02  Arnold  et  al.  [36] 


a  restricted  open-shell  Hartree-Fock  (ROHF)  calculation  in  the  Abelian  subgroup 
D2/1  provides  a  spin  eigenfunction,  it  cannot  provide  tt^  and  iTy  orbitals  of  the  proper 
symmetry.  Therefore,  quasi-restricted  Hartree-Fock  (QRHF)  calculations  were  at- 
tempted, one  using  the  orbitals  from  the  C3  RHF  calculation,  and  one  from  Cl" ,  the 
next  ^E^  species,  which  has  a  filled  iig.  While  the  Cl~  orbitals  reproduce  the  C^ 
CCSD  energy,  there  are  several  large  T  amphtudes  in  the  CCSD  wavefunction,  indi- 
cating these  orbitals  are  not  suitable  for  the  EOMIP  CCSD  calculation.  The  result 
is  a  negative  electron  detachment  energy.  The  C3  orbitals  gave,  to  within  .002%,  the 
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same  Cj  CCSD  energy,  an  electron  detachment  energy  within  .05  eV  of  the  ACCSD 
value,  and  fortuitously  closer  to  the  experimental  energy. 

11.3  First  Ionization  Energies 

The  carbon  clusters  were  studied  in  several  different  basis  sets  as  a  preUminary  to 
the  study  of  their  full  spectra.  These  preliminary  results  indicate  considerable  effects 
of  orbital  choice.  The  orbitals  from  the  restricted  Hartree-Fock  calculation  on  the 
neutral  (RHF  in  the  cases  of  C2,  C3,  and  C5,  ROHF  in  the  cases  of  the  open-shell 
neutral  caxbon  clusters)  seemed  to  give  results  closest  to  the  ACCSD  energies,  and 
so  were  initially  selected  for  the  extended  studies  of  the  ionization  spectra.  But  these 
were  not,  in  all  cases,  the  choices  which  best  reproduced  the  experimental  spectra,  as 
can  be  seen  from  the  simulated  ionization  spectra  presented  in  section  11.4. 

11.4  Interpretation  of  Spectra 

Lorentzian  line  shapes  were  used  to  simulate  ionization  spectra,  with  a  FWHM 
of  .25eV  used  in  all  cases  except  for  the  C^  spectra,  for  which  a  FWHM  of  .2eV  was 
used. 


Table  11.4:  Cn/C". 


Method 

Orbitals 

C-  CCSD  (au) 

VEDE  (eV) 

Source 

C3-  2n,/C3 

■'^t 

pVTZ+sp 

EOMIP 

C3-  UHF 

-113.915415 

2.157 

pVTZ+sp 

EOMIP 

C3RHF 

-113.913530 

1.989 

pVTZ+sp 

CCSD 

C3"  UHF 

-113.915415 

2.05 

Watts  and  Bartlett  [82] 

pVTZ+sp 

CCSD(T) 

C3"  UHF 

1.88 

Watts  and  Bartlett  [82] 

pVTZ+sp 

CCSD 

C^  ROHF 

1.90 

Watts  and  Bartlett  [82] 

Experiment 

1.98+  0.02 

Arnold  et  al.  [36] 

C4-  2n,/C4 

3v- 
^9 

DZP+sp 

EOMIP 

Cf-  RHF 

-151.758857 

4.55 

DZP+sp 

EOMIP 

Cl+  RHF 

-151.756564 

2.08 

DZP+sp 

EOMIP 

C4  ROHF 

-151.756460 

2.78 

DZP+ sp 

EOMIP 

C4  ROHF 

-151.757755 

3.47 

DZP+sp 

EOMIP 

C4"  UHF 

-151.758929 

3.43 

DZP+sp 

CCSD 

3.73 

Watts  and  Bartlett  [65] 

DZP+sp 

CCSD{T) 

3.76 

Watts  and  Bartlett  [65] 

pVTZ+sp 

EOMIP 

C4"  ROHF 

-151.875081 

3.66 

pVTZ+sp 

EOMIP 

C4"  UHF 

-151.876304 

3.63 

pYTZ+sp 

EOMIP 

C4  ROHF 

-151.873807 

3.51 

pVTZ+ 

CCSD 

ROHF 

3.91 

pVTZ+ 

CCSD(T) 

ROHF 

3.95 

Experimeni 

3.882+0.01 

Arnold  et  al.  [36] 

C5-  2n„/C5 

'%" 

DZP+sp 

EOMIP 

C5  UHF 

-189.732655 

2.91 

DZP+sp 

EOMIP 

C5  ROHF 

-189.731051 

2.83 

DZP+sp 

EOMIP 

C5RHF 

-189.730180 

2.71 

DZP+sp 

CCSD 

2.79 

Watts  and  Bartlett  [65] 

DZP+sp 

CCSD(T) 

2.50 

Watts  and  Bartlett  [65] 

pVTZ+sp 

EOMIP 

C5  UHF 

-189.885918 

3.14 

pVTZ+sp 

EOMIP 

C5  ROHF 

-189.884250 

3.06 

pVTZ+sp 

EOMIP 

C5RHF 

-189.883307 

2.93 

pVTZ+sp 

CCSD 

UHF 

3.02 

Watts  and  Bartlett  [82] 

pVTZ+sp 

CCSD(T) 

UHF 

2.74 

Watts  and  Bartlett  [82] 

pVTZ+sp 

CCSD 

ROHF 

2.97 

Watts  and  Bartlett  [82] 

pVTZ+sp 

CCSD(T) 

ROHF 

2.78 

Watts  and  Bartlett  [82] 

Experiment 

2.839+0.008 

Arnold  et  al.  [36] 

Table  11.4:  C„/C;  {Continued) 


Method 


Orbitals         C"  CCSD  (au)      VEDE  (eV)     Source 


Cg-  2n„/c6  'S5 

DZP+sp     EOMIP 

Ce  ROHF 

-227.706504 

5.99 

DZP+5P     EOMIP 

0^+  RHF 

-227.706974 

5.63 

DZP+sp     EOMIP 

Cg  ROHF 

-227.707420 

4.35 

DZP+sp     EOMIP 

Cg  UHF 

-227.708828 

4.31 

DZP+sp     EOMIP 

Ce  UHF 

-227.708080 

4.14 

DZP+sp     CCSD 

4.17 

DZP+sp     CCSD(T) 

4.15 

Experiment 

4.185+0.006 

Cf  ^n^/Cy  ^s+ 

DZP+sp     EOMIP 

C7RHF 

-265.661618 

3.31 

DZP+sp     EOMIP 

Cf  UHF 

-265.693376 

3.56 

DZP+sp     EOMIP 

Cf  ROHF 

-265.692714 

3.43 

DZP+sp     CCSD 

3.4 

DZP+sp     CCSD(T) 

3.01 

Experiment 

3.538+0.014 

Cg-  2n,/C8  ^s- 

DZP+sp     EOMIP 

Cg  ROHF 

-303.633302 

4.50 

DZP+sp     EOMIP 

Cg  ROHF 

-303.635412 

4.63 

DZP+sp     EOMIP 

Cg  UHF 

-303.637055 

4.59 

DZP+p      CCSD 

4.45 

DZP+p      CCSD(T) 

4.39 

Experiment 

4.379+0.006 

Cg-  2n„/C9  ^s+ 

DZP+sp     EOMIP 

C^  UHF 

-341.589514 

4.04 

DZP+sp     EOMIP 

Cg  ROHF 

-341.587993 

3.85 

DZP+sp     EOMIP 

Cg  RHF 

-341.586665 

3.76 

DZP+sp     CCSD 

UHF 

-341.589514 

3.83 

DZP+p      CCSD 

3.80 

DZP+p      CCSD(T) 

3.31 

Experiment 

3.684+0.01 

Cj-Q  2n„/Cio  'S5 

DZP+sp     EOMIP 

Cfo  ROHF 

-379.558974 

4.88 

DZP+sp     EOMIP 

C^o  UHF 

-379.560459 

4.80 

DZP+sp     EOMIP 

Cfo  ROHF 

-379.557551 

4.73 

DZP+sp     MBPT(2) 

UHF 

-379.520429 

4.72 

DZP+sp     CCSD 

UHF 

-379.560459 

4.68 

DZP+p      CCSD 

4.71 

DZP+p      CCSD(T) 

4.59 

Watts  and  Bartlett  [65] 
Watts  and  Bartlett  [65] 
Arnold  et  al.  [36] 


Watts  and  Bartlett  [65] 
Watts  and  Bartlett  [65] 
Arnold  et  al.  [36] 


Watts  and  Bartlett  [65] 
Watts  and  Bartlett  [65] 
Arnold  et  al.  [36] 


Watts  and  Bartlett  [65] 
Watts  and  Bartlett  [65] 
Arnold  et  al.  [36] 


Watts  and  Bartlett  [65] 
Watts  and  Bartlett  [65] 
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11.4.1     Odd-Length  Chains 


For  the  odd  clusters  C^  and  Cf ,  and  C^,  where  there  is  only  one  electron  in 
an  open-shell  tt  orbital,  the  QRHF  orbitals  based  on  the  closed  shell  give  very  good 
results  for  the  first  ionization,  as  might  be  expected.  Where  QRHF  or  ROHF  orbitals 
were  used  for  the  C5  and  Cf,  spectra,  the  EOMIP  results  are  a  good  match  to 
the  experimental  spectrum.  Anion  UHF  orbitals,  used  in  C^  gave  decidedly  poorer 
results.  So  selection  of  appropriate  orbitals  is  important.  Large  T  amplitudes  and 
large  correlation  energies  indicate  that  the  orbitals  are  poorly  chosen,  and  that  the 
EOMIP  results  will  be  poor  as  well. 


11,4,2    C3" 


The  experimental  ionization  spectrum  of  C^  is  very  rich.  The  two  main  features 
(the  peaks  at  1.95  and  4.3  in  the  Smalley  C3  spectrum,  1.995  and  4.113  in  the  Neu- 
majk  spectrum)  are  easily  assigned  to  the  X^E^  and  a^n^  states  of  C3,  in  agreement 
with  the  assignment  in  Neumark.  The  feature  at  about  5.0  eV  seems  to  be  well 
matched  to  the  ^Ylg  state  at  about  the  same  energy,  with  the  succeeding  gap  and 
calculated  ^11+  peak  at  5.286  eV  relatively  well  matched  to  the  gap  and  peak  in  the 
SmaJley  spectrum  at  about  that  frequency.  The  complicated  feature  at  about  6.0 
eV  in  the  Smedley  spectrum  is  difficult  to  assign.  There  are  a  number  of  large  peaks 
in  the  region  which  may  all  contribute.  Chief  among  these  is  the  ^S"  peak  at  5.6  eV. 


Table  11.5:  Ionization  of  C3  X^IIg  in  pVTZ+sp  basis,  QRHF  reference 


Energy 

eV  (Intensity) 

C3  State 

Reference 

EOMIP 

Other 

Experiment 

XiE+ 

1-Kg 

1.989 

(.883) 

1.90" 

1.98  ±  0.02* 
1.95^ 

^n„ 

iw: 

4.172 

(.879) 

W4.113'' 
w4.3= 

in„ 

^i 

4.616 

(.842) 

^n. 

4aJ 

5.056 

(.848) 

'^t 

ItTux 

5.286 

(.903) 

^n. 

H 

5.560 

(.805) 

^s- 

iTTuy 

5.595 

(.758) 

«5.1  -  e.o'^ 

^s- 

iTTuy 

5.683 

(.907) 

'K 

ItTux 

5.821 

(.835) 

^Ug  2a^,  6a+h(7+la^  15.163  (.224) 
3n„  3a+,  4a+t3CT+4CT+  17.498  (.325) 
^Uu     3a+   44^4(7+ l;r„x     17.541     (.495) 


°  Watts  and  Bartlett  [82]  ROHF-CCSD(T)/PVTZ+sp; 
*■  Arnold  et  al.  [36];  "  Yang  et  al.  [35] 

Because  the  EOMIP  code  used  is  not  spin-adapted,  open-shell  singlets  such  as  this 

one  are  not  properly  represented,  and  it  is  expected  that  the  singlet- triplet  sphtting 

is  under- represented  by  a  factor  on  the  order  of  22  (cf.  work  by  Nooijen  et  al.).  It 

is  reasonable  to  expect  that  the  ^'S,'  peak  as  well  as  the  nearby  ^S^  at  5.8  eV  will 

contribute  to  this  feature. 

It  is  important  to  note  that  all  of  the  calculated  transitions  in  this  range  are  very 
strongly  dominated  by  a  single  1-hole  process.  Significant  calculated  shakeup  peaks 
first  appear  in  the  15-18  eV  region. 

Because  the  EOMIP  program  lacks  explicit  spin  adaptation,  the  only  open-shell 
wavefunctions  which  are  properly  spin  adapted  are  those  with  aJl  parallel  spins  in 
the  half-filled  orbitals.   This  means  that  the  open-shell  singlet  states  of  C3  (all  but 


Table  11.6:  Ionization  of  X^n^  C3"  in  pVTZ+5p  basis,  UHF  basis 


Energy  eV 

(Intensity) 

C3  State 

Reference 

EOMIP 

Experiment 

XiE+ 

iTTg 

1.989 

1.98  ±  0.02* 
1.95^^ 

^n. 

3a+ 

4.172 

5=4.113'' 

^n„ 

3<7+ 

4.616 

^n. 

4a; 

5.056 

(l.OOO'^) 

'K 

iTTux 

5.286 

^n. 

4a; 

5.559 

3t- 

iTTuy 

5.595 

^Su 

iTTuy 

5.638 

'K 

ItTux 

5.821 

«5.1  -  6.0^^ 

'^t 

i-Kg,  ml  i¥q  3cr; 

10.351 

(0.002'^) 

^n„ 

2< 

15.163 

^n. 

H 

17.541 

Watts  and  Bartlett  [82]  UHF-CCSD{T)/PVTZ+sj); 
*■  Arnold  et  al.  [36];  '  Yang  et  al.  [35];  ''  Relative  to  ^IIp  state  at  5.056  eV 
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Figure  11.1:  Ionization  spectrum  of  C3  ,  pVTZ+sp  basis 


the  lowest-energy,  closed-shell  X^E^  state)  are  spin-contaminated.  For  example,  if 
the  X^IIu  state  has  one  unpaired  electron  in  the  Itt^  orbital  of  beta  spin,  the  EOMIP 
open-shell  singlet  wavefunction  describing  the  ^U.^  wavefunction  is  dominated  by  the 
annihilation  of  the  Sa^  beta  electron,  leaving  an  open-shell  I'Wg  Sa^^  configuration. 
But  the  proper  representation  of  the  spin  state  demands  that  the  spin-flipped  con- 
figuration iTTg  Sa-u  be  represented  with  equal  weight  in  the  wavefunction.  Because 
this  state  can  be  reached  via  the  EOMIP  doubles  operator  Cj-^^^  iTr^Scr^lTfj,  it  does 
appear  in  the  final  wavefunction.  But  correlation  of  this  2-configuration  reference 
requires  excitations  from  both  configurations,  and  the  required  triple  excitation  op- 
erators are  not  present  in  the  EOMIP  operator.  This  lack  of  balance  -  essentially 
treating  two  equally  important  configurations  at  different  levels  of  theory  -  makes  it 
impossible  to  correlate  the  state  while  at  the  same  time  maintaining  a  spin  eigen- 
function.  Inclusion  of  triple  excitations  in  the  method  would  improve  the  situation 
by  pushing  the  problem  back  one  level  of  theory,  where  correlation  contributions  will 
be  smaller,  in  all  but  very  unusual  cases.  For  larger  molecules,  the  error  induced  in 
the  energy  by  the  incomplete  treatment  of  spin  adaptation  may  be  barely  significant 
compared  to  experimental  errors,  but  molecules  as  small  as  C3  are  known  to  exhibit 
larger  singlet-triplet  splittings. 

Two  sets  of  two-determinant  MRCCSD  calculations  were  performed  -  one  us- 
ing QRHF  orbitals  from  the  closed-shell  C3  SCF  calculation,  and  the  other  using 
orbitals  chosen  to  be  more  suited  to  each  ionization.  An  important  feature  of  the 
two-determinant  calculations  (figure  11.1)  is  the  very  close  grouping  of  ionization  en- 
ergies in  the  5.0  to  5.6  eV  region.  This  results  in  one  tall  peak  at  around  5.5  eV  rather 


120 

than  several  maxima  in  the  4.5  to  6eV  region.  The  EOMIP  spectrum,  although  it  has 
similar  difficulties,  seems  to  represent  the  spectrum  slightly  more  realistically,  with 
a  better  matching  of  triplet  peaks  to  the  experimental  features.  Too-small  singlet- 
triplet  splittings  might  be  held  responsible  for  the  very  close  spacing  of  lEs  in  the 
5.5  to  6  eV  region.  It  is  not  easy  to  account  for  the  lack  of  any  calculated  peak 
corresponding  to  the  feature  above  7  eV. 


11-4.3    c; 


The  calculated  peaks  in  figure  11.2  are  offset  slightly  from  the  experimental  in  a 
manner  consistent  with  results  for  small  closed-shell  molecules. 

The  rising  section  at  the  high-energy  end  of  the  spectrum  is  explained  by  energies 
which  fall  outside  the  0-7.5  eV  window  observed  in  the  Smalley  experiment. 

In  table  11.7,  the  first  ionization  energy  is  within  .2  eV  of  the  experimental  results, 
and  the  theoretical  states  between  5.3  and  6.0  eV  result,  after  the  application  of  Lorenz 
broadening,  in  a  peak  at  5.46  eV,  within  1  eV  of  the  experimental  peak  at  5.39  eV. 


11-4.4    Cj 


Like  C5,  C7  is  a  closed-shell  species,  and  again  the  calculated  spectrum  bears  a 
distinct  resemblance  to  experiment.  Peak  shapes  suggest  that  lEs  greater  than  5.5 
eV  are  too  large. 
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Figure  11.2:  Ionization  of  C^,  ROHF  basis 


Table  11.7: 

C5  Ionization.  ROHF  orbitals 

Energy  eV 

(Intensity) 

C5  State 

EOMIP 

Expt. 

2.839±0.008" 

IS." 

3.058 

(0.945) 

2.8* 

^n,. 

5.359 

(0.873) 

'K 

5.362 

(0.928) 

'n„. 

5.455 

(0.872) 

5.29^ 

^2„- 

5.473 

(0.934) 

^s- 

5.476 

(0.928) 

'n^x 

5.708 

0.854 

'n„^ 

5.816 

0.849 

'K 

5.932 

0.889 

^s; 

7.860 

0.898 

As  reported  by  Arnold  et  al.  [36] 
^  As  reported  by  Yang  et  aJ.[35];  '  Maxima  observed  in  plots  by  Yang  et  al.  [35] 
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Figure  11.3:  Ionization  of  Cf ,  ROHF  basis 


Table  11.8:  Cf  Ionization,  ROHF  orbitals 


Energy  eV  (Intensity) 

C7  State 

EOMIP 

Expt. 

'^t 

3.425 

(0.907) 

3.538±  0.014'' 
3.10'' 
5.0(.9)= 

^S- 

5.130 

(0.910) 

'^t 

5.174 

(0.925) 

^s- 

5.187 

(0.919) 

5. 25 (.82  shoulder) = 

^n„ 

5.537 

(0.854) 

^n. 

5.545 

(0.855) 

'^t 

5.685 

(0.891) 

^n„ 

5.813 

(0.841) 

^n. 

5.824 

(0.838) 

6.8(.3)= 

'^t 

7.175 

(0.883) 

^s.- 

7.359 

(0.885) 

^Sp 

7.461 

(0.885) 

^^1 

7.727 

(0.868) 

"  As  reported  by  Arnold  et  al.  [36] 

''  As  reported  by  Yang  et  al.  [35] 

'^  Maxima  observed  in  plots  by  Yang  et  al.  [35] 


These  calculations  were  performed  using  ROHF  orbitals,  in  the  DZP+sp  basis 
because  of  the  size  of  the  calculation.  Using  D2/1  symmetry,  there  were  129  a  and  129 
P  functions  in  the  SCF  calculation.  Dropping  the  5  lowest-energy  orbitals  of  each  spin 
reduced  that  number  to  124  each.  Performing  the  calculation  in  the  pVTZ+sp  basis 
used  for  the  smaller  linear  carbon  clusters  would  have  meant  over  250  basis  functions 
of  each  spin. 

The  first  calculated  ionization  appears  identical  with  the  first  experimental  ion- 
ization. Although  Small ey's  group  reports  a  vertical  electron  affinity  of  3.10  eV,  the 
available  experimental  plot  has  a  maximum  ma.ximum  closer  to  3.40  eV^,  which  gives 
almost  perfect  overlap  between  the  experimental  plot  and  the  EOMIP  results.  This 
apparent  improvement  in  the  position  of  the  first  ionization  energy  compared  to  the 
C^  calculation  is  consistent  with  the  results  for  smaller  molecules,  where  reduction  of 
the  valence-type  polarization  functions  leads  to  slightly  smaller  ionization  energies, 
£ind  slightly  better  agreement  between  theory  and  experiment,  at  least  for  the  first 
ionization. 

11.4.5     Q- 

The  disagreement  of  the  UHF  EOMIP  calculation  with  experiment  (figure  11.4) 
can  be  attributed  to  the  unsuitability  of  the  orbitals.  Among  the  CCSD  amplitudes, 
there  is  one  large  TjOf  element,  a  iTg^  valence  excitation  with  amphtude  .128.  In 
the  QRHF  CCSD  calculation,  the  situation  is  similar,  and  would  presumably  lead  to 
similar  EOMIP  results 
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Figure  11.4:  Ionization  spectrum  of  Cg 


Table  11.9 

Cg  Ionization,  1 

JHF  orbitals 

Energy  eV 

(Intensity) 

Cg  State 

EOMIP 

Expt. 

3.684±0.01° 

1^.^ 

4.041 

(0.928) 

3.70* 

IS- 

4.997 

(0.909) 

b.O'^ 

3E: 

5.227 

(0.927) 

3E+ 

5.267 

(0.931) 

IS^ 

5.664 

(0.891) 

SHp 

5.693 

(0.833) 

3n„ 

5.694 

(0.833) 

m, 

6.112 

(0.840) 

in„ 

6.113 

(0.846) 

6.5^ 

3E; 

6.875 

(0.908) 

3s; 

6.985 

(0.900) 

is; 

7.109 

(0.862) 

1^1 

7.429 

(0.859) 

°  As  reported  by  Arnold  et  cd.  [36] 

''  As  reported  by  Yang  et  al.  [35] 

'^  Maxima  observed  in  plots  by  Yang  et  al.  [35] 


On  the  beta  side,  the  UHF  CCSD  also  has  a  large  tt^i  valence  excitation,  with 
amplitude  -.115,  while  QRHF  has  one  with  amphtude  -.058,  which  is  not  as  severe. 
It  would  seem  that  the  size  of  the  Ti  amplitudes  make  the  UHF  orbitaJs  the  worst 
possible  choice  from  the  point  of  view  of  the  size  of  the  correlation  correction  - 
AEqcsD  is  -1-131  for  UHF,  -1.068  for  QRHF,  and  -1.040  for  ROHF.  Only  the  ROHF 
CCSD  calculation  doesn't  show  any  T  amplitudes  over  .1,  which  should  indicate  that 
EOMIP  results  would  be  improved  if  ROHF  orbitals  were  used. 

11.4.6    Even- Length  Chains 

A  pervasive  problem  with  calculations  on  Cj,  C^,  and  C^  is  the  partially-filled 
acetylenic  valence  tt  orbital,  which  inevitably  leads  to  symmetry  breaking  because  the 
calculations  are  carried  out  in  the  D2/1  Abelian  subgroup  of  the  full  Doo/i  molecular 
symmetry.  While  an  ROHF  calculation  on  the  ion  will  ensure  that  alpha  and  beta 
orbitals  are  identical,  it  cannot  ensure  that  the  tTx  and  iTy  components  of  the  partially 
filled  valence  tt  orbital  will  be  identical.  And  while  a  QROHF  calculation  based  on  the 
ground  state  neutral  will  provide  properly  symmetric  tt^  and  tt^  orbitals,  the  mere  fact 
that  they  are  occupied  differently  in  the  reference  will  break  the  symmetry.  The  SRCC 
calculation  cannot  promote  electron  occupation  from  one  irreducible  representation 
to  another,  and  so  cannot  "fix  up"  the  broken  symmetry.  This  means  that  certain 
states  which  should  be  degenerate  in  energy  will  not  be.  For  example,  Arnold  et  al. 
[36]  cite  experimental  and  theoretical  evidence  that  C4,  Ce,  and  Cs  have  low- lying 
'A  excited  states.    These  arise  in  the  photoelectron  spectrum  when  an  electron  is 
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removed  from  a  tt  orbital  other  than  the  one  partially  filled  in  the  ground  state  ion 
or  neutral.  The  possible  results  are  a  E^,  a  E~,  and  a  A  state.  In  the  D2/1  symmetry 
group,  these  choices  reduce  to  two  E^  and  two  E"  states.  If  the  molecular  symmetry 
were  preserved,  one  of  the  S"  states  would  clearly  be  degenerate  with  one  of  the  E""" 
states,  with  characteristic  sjmimetry  between  the  tt^  and  TTy  orbitals.  In  CC  EOMIP 
calculations  on  these  even-length  chains,  the  energies  of  such  states  will  be  split, 
and  the  wavefunctions  will  not  show  exact  x-y  symmetry.  However,  in  cases  of  near 
degeneracy  where  the  wavefunctions  show  something  approaching  the  proper  form,  it 
may  be  presumed  that  the  states  are  in  fact  an  approximation  to  the  proper  A  state. 

U.4.7    Ca" 

EOMIP  calculations  in  UHF  and  QRHF  orbitals  (table  11.10  and  figure  11.4.7) 
tell  approximately  the  same  story.  The  relationship  between  the  experimental  and 
calculated  spectra  is  shght.  In  addition  to  the  afore-mentioned  multi-reference  prob- 
lem (section  10.4),  there  is  the  problem  of  spin  contamination  in  the  open-shell  singlet 
states  of  C2.  With  the  exception  of  the  ^EJ"  line,  all  the  so-called  singlet  states  are 
in  fact  severely  spin  contaminated,  with  an  S^  expectation  value  of  1.4.  This  means 
that  the  ^n„  fine  should  be  expected  to  move  up  in  energy,  narrowing  the  first  peak, 
and  moving  its  maximum  slightly  lower  in  energy,  as  well  as  decreasing  the  relative 
height  of  this  peak. 

The  two-determinant  MRCCSD  calculation,  hke  the  EOMIP  calculation,  incor- 
rectly places  the  C2  a'^II^  state  lower  in  energy  than  the  ground  X^E^"  state.  Triples 


corrections  eire  required  to  correctly  order  these  states.  The  two-determinant  calcula- 
tion gives  a  singlet-triplet  splitting  of  the  11  states  similar  to  that  given  by  the  EOMIP 
calculation,  but  gives  a  much  larger  splitting  for  the  S^  states  (4.6  eV  versus  1.17eV 
for  the  EOMIP).  With  the  ^E+  state  near  9  eV,  the  two-determinant  calculation  gives 
a  simulated  spectrum  much  more  similar  to  experiment  than  the  EOMIP  calculation 


does. 
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Figure  11.5:  C2  Ionization  spectrum,  EOMIP  calculation 
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Table  11.10:  C2" 

in  pVTZ+sp  basis. 

QRHF  orbital 

s 

Energy  eV  (Intensity) 

C2- 

C2 

Reference 

EOMIP 

Other 

Experiment 

x^e; 

X1S+ 

3< 

4.05(.849) 

3.359" 

3.273^ 
?^  4.4*^ 

a^n„ 

ItTu 

3.450(.891) 

3.059* 

3.35r 

^m 

^: 

4.600(.893) 

4.446'' 

^n„ 

l7r„ 

4.954(.863) 

4.090'' 

'K 

2< 

5.770(.640) 

9.074'' 

'^t 

2a;,  27rj2a+l7r„ 

11.952(.057) 

'^t 

2^;,  3a;t2a;r2a+ 

12.221(.100) 

'K 

2<,  27fgt27f„3a; 

12.395(.274) 

'^^ 

2a;,  27rt2<l7f„ 

14.015(.516) 

'^: 

2a;.  3ajt2a:2a,+ 

15.478(.175) 

9 

2a;,  3a;t2a;3a; 

19.69(.422) 

"  ACCSD  in  pVTZ+sp  basis 

''  TWODET  MRCCSD  calculation,  QRHF  orbitals 

'  Arnold  et  al.  [36] 

'  Yang  et  al.  [35] 


11.4.8    C 


In  the  spectrum  of  figure  11.7,  one  can  confidently  assign  the  ^Il„  line  to  the  left 
shoulder  of  the  first  experimental  peak,  but  beyond  that  it  is  difficult  to  associate 
calculated  with  experimental  features.  Arnold  et  al.  assign  a  peak  at  4.209  eV  to 
a  ^A  state  of  C4  In  the  QROHF  orbitals,  the  character  of  the  two  wavefunctions 
differ  enough  that  they  do  not  approximate  a  A  state.  The  ^E+  state,  for  example, 
dominated  by  the  •  •  •  iTTgxVffgx  configuration,  has  the  accompanying  •  •  •  iTigyVffgy  con- 
figuration with  the  wrong  sign  to  indicate  a  A  state.  And  the  ^E~  state  is  perhaps 
closer  to  being  a  low-spin  triplet  than  a  singlet  state 

In  the  ROHF  orbitals,  however,  the  spin  contamination  and  symmetry  breaking 
problems  are  somewhat  reduced  in  the  two  EOMIP  wavefunctions,  which  do,  when 
taken  together  resemble  a  ^A  state,  supporting  the  assignment  made  by  Arnold  et  al. 
While  the  correlation  correction  in  the  ROHF  orbitals  is  .03  Hartrees  smaller  than 
that  in  the  QRHF  orbital  set  (the  total  energies  agree  to  a  milli Hartrees),  some  of  the 
largest  T2  amplitudes  in  the  ROHF  set  are  up  to  50%  larger  than  the  corresponding 
QRHF  amplitudes,  indicating  a  possible  source  of  error  in  ROHF-based  calculations. 
It  is  impossible,  based  on  this  evidence  to  draw  a  firm  conclusion  on  the  assignment 
of  the  peak. 
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Figure  11.7:  Cj  ionization  spectrum,  ROHF  orbitals 


Table  11.11:  Ionization  of  C^  in  pVTZ+sp  basis,  ROHF  orbitals 


Energy  < 

sV  (Intensity) 

C4  State 

Reference 

EOMIP 

Experiment 

3.882±  .01" 

'^.- 

IWg 

3.655 

(.944) 

3.70" 

'Big 

iT^gx 

3.843 

(.935) 

^A, 

l^sy 

4.093 

(.909) 

4.209" 
4.5^^ 

^n, 

5^ 

5.057 

(.875) 

^n„ 

4^ 

5.377 

(.869) 

^n. 

5a; 

5.439 

(.855) 

^n„ 

4^: 

5.795 

(.837) 

6.4^ 

3v  + 

l7f„y 

6.598 

(.925) 

^^u 

iTfui 

7.106 

(.882) 

'K 

iTTux 

7.151 

(.880) 

'K 

ItTuj, 

7.341 

(.871) 

"  As  reported  by  Arnold  et  al.  [36]  *■  As  reported  by  Yang  et  ad.  [35] 
^  Maxima  observed  in  plots  by  Yang  et  al.  [35] 
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Figure  11.8:  C4  ionization  spectrum,  QROHF  orbitals 


Table  11.12:  Ionization  of  C^  in  p\TZ+sp  basis,  QROHF  orbitals 


Energy  ( 

3V  (Intensity) 

C4  State 

Reference 

EOMIP 

Experiment 

3.882±  .or 

^^.- 

l7fg 

3.512 

(.904) 

3.70* 

^^1. 

iT^gx 

3.698 

(.896) 

4.209" 

'A9 

iTTgy 

3.965 

(.876) 

4.5^ 

^n 
^^9 

5^ 

4.912 

(.841) 

^n„ 

4a: 

5.234 

(.841) 

^n. 

5< 

5.281 

(.814) 

^n„ 

^< 

5.637 

(.798) 

6.4^ 

'^: 

iTTuy 

6.415 

(.888) 

'K 

VKux 

6.888 

(.849) 

'^z 

ItTui 

6.92 

(.840) 

'K 

iTTuy 

7.11 

(.827) 

"  As  reported  by  Arnold  et  al.  [36]  *"  As  reported  by  Yang  et  al.  [35] 
^  Maxima  observed  in  plots  by  Yang  et  al.  [35] 
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Figure  11.9:  Ionization  spectrum  of  Cg 


Table  11.13:  Ionization  of  Cg  ,  ROHF  orbitals 


State 

EOMIP 

Expt. 

Orbitals 

4.185+0.006° 

^^. 

4.347 

(0.928) 

4.1f 

f 

27f„x 

4.459 

4.547 

(1.000) 

(0.884) 

4.36" 

5.0^ 
5.2^ 

CA) 

ItTux 

l7r„j, 

'■'■VlX 

5.344 

(0.856) 

7a, 

'n,. 

5.377 

(0.853) 

6a„ 

^n„. 

5.608 

(0.834) 

6(Tg 

'n^x 

5.644 

(0.833) 

5.7'^ 

7a„ 

^s: 

5.750 

(0.933) 

ItTjx 

^s- 

6.530 

(0.878) 

^"^gy 

'K 

6.554 

(0.868) 

2-Kgx 

'K 

6.633 

(0.841) 

2^?:/ 

°-  As  reported  by  Arnold  et  al.  [36]  *  As  reported  by  Yang  et  al.  [35] 
'^  Maxima  observed  in  plots  by  Yang  et  al.  [35] 


11.4.9  a 


The  EOMIP  value  for  the  first  ionization  of  Cg  is  within  .2  eV  of  the  experimental 
value.  Beyond  that,  the  assignment  of  peaks  becomes  more  difficult,  although  the 
theoretical  and  experimental  spectra  (figure  11.9  and  table  11.13)  are  more  similar 
here  than  they  are  for  C^.  Assigning  the  cluster  of  EOMIP  states  between  5.344  and 
5.644  eV  to  the  two  experimental  peaks  at  about  5.0  and  5.2  eV  implies  EOMIP  results 
in  error  by  .3  to  .44  eV  in  this  region.  While  it  would  be  convenient  to  assign  the  ^S;!' 
peak  at  5.75  eV  to  the  experimental  5.70  peak,  an  error  this  small  seems  unlikely  in 
light  of  the  errors  at  lower  energy.  Therefore,  it  would  seem  that  the  double-peaked 
experimental  structure  centered  at  about  5.1  eV  is  composed  of  transitions  to  the 
^IIu,  ^Ug,  ^IIu,  ^lig,  and  ^E+  states  at  energies  ranging  from  5.344  to  5.750  eV.  In 
fact,  the  ^E^  could  be  related  to  the  small  shoulder  on  the  high-energy  side  of  the 
5.2  eV  peak.  This  leaves  the  remaining  peaks  in  the  range,  ^£„ ,  ^E^ ,  and  ^E~,  to 
be  assigned  to  the  experimental  5.7eV  peak,  with  an  error  somewhere  in  the  range 
of  .83  to  .96  eV.  The  ^Eg /'E+  pair  of  peaks,  separated  by  only  .09  eV,  show  some 
A  characteristics.  Although  the  ratio  of  the  two  7r„  configurations  is  a  little  under 
10:1  for  the  Ci  versus  the  C2  amplitude,  the  signs  are  consistent  with  a  ^A  state, 
and  the  energies  are  close  enough  together  to  to  be  considered  nearly  degenerate. 
The  Eu  pair  of  states  near  7  eV,  on  the  other  hand,  despite  their  nearly-degenerate 
energies,  do  not  display  the  characteristics  of  a  ^A  state,  and  show  all  evidences  of 
being  approximations  to  the  true  ^E„  and  ^E+  states. 


Cg ionizations,  EOMIP,  QROHF  orbitals 


0.9- 
0.&- 

0.: 

0.6H 
0.1 
0.. 
0.3H 

0.; 

0. 


Smalley 

EOMIP   


3^g1^g1SS 


\        l'5        I 


2'5 


i      s's      1 


Figure  11.10:  Ionization  spectrum  of  Cg  ,  QROHF  orbitals 


Table  11.14:  Cg  Ionization,  QROHF  orbitals 


Energy  eV 

(Intensity) 

State 

EOMIP 

Expt. 

^^.- 

4.498 

(.889) 

4.379±0.006° 
4.42*- 

4.566 
4.623 

(.882) 
(.871) 

4.449°  (^A) 
5.2'= 

'K 

5.370 

(.866) 

5.4'= 

5.654 
5.657 
5.849 
5.853 
6.184 

(.839) 
(.843) 
(.814) 
(.806) 
(.852) 

6.200 
6.266 

7.757 

(.843) 
(.826) 
(.871) 

*  As  reported  by  Arnold  et  al.  [36]  ^  As  reported  by  Yang  et  al.  [35] 
■^  Maxima  observed  in  plots  by  Yang  et  al.  [35] 
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11.4.10  g- 


The  resemblance  between  the  theoretical  and  experimental  pealcs  for  the  ionization 
of  Cg  is  even  more  striking  than  it  is  for  the  smaller  even- length  chains.  The  difference 
between  the  calculated  first  ionization  and  that  of  the  next  two  calculated  energies, 
at  4.566  and  4.623  eV,  belong  under  the  first  experimental  peak.  For  C^,  as  for  the 
C^  QROHF  calculation,  these  two  states  do  not  exhibit  ^A  characteristics,  and  in 
act  appear  to  be  much  more  spin  contaminated  and  symmetry  broken  than  the  ^A 
Ce  state  wavefunctions  were.  It  appears  that  the  choice  of  QROHF  orbitals  causes 
additional  difficulty  in  overcoming  the  spin  and  symmetry  contamination  problems 
of  the  reference. 


11.4.11     Summarv 


The  ROHF  orbitals  give  good  agreement  with  experiment  for  the  first  ionization 
of  C4  ,  Ce  ,  and  Cg  in  the  DZP+sp  basis  set.  For  higher-energy  ionizations,  there 
seems  to  be  a  problem  with  the  distribution  of  the  peaks.  These  even-length  chains 
show  considerable  differences  in  bond  lengths  and  the  amount  of  bond-length  alter- 
nation between  the  anion  and  the  ground-state  neutral.  There  is,  as  yet,  no  in-depth 
theoretical  analysis  of  the  bonding  character  or  the  electron  density  distribution  of 
any  of  the  excited  states  of  the  neutrals.  Watts  and  Bartlett  [65]  observed  that  the 
anions  exhibit  considerable  bond-length  alternation,  midway  between  the  cumulenic 
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neutral  ground  state  structure  and  the  acetylenic  dianion  structure.  This  suggests 
that  there  are  considerable  differences  in  tt  orbital  densities  among  these  states,  and 
likely  among  the  various  neutral  excited  states.  Bernholdt  noted  that  the  size  and 
direction  of  the  FSCC  third-order  triples  approximation  differed  between  n  and  a 
ionizations  of  C5.  It  would  require  a  more  detailed  examination  of  electron  density 
and  bonding  in  these  states  to  answer  these  questions. 


CHAPTER  12 
SUMMARY 

In  this  work,  the  author  has 

•  Derived  the  CC  EOMIP  method  for  the  study  of  electron-detached  states.  The 
implementation,  in  the  CC  EOMIP  program,  uses  Nakatsuji's  generalization  of 
the  Davidson  diagonalization  method  in  a  thrifty  and  fast-converging  approach 
to  studjang  several  states  of  a  given  molecule  in  a  consistent  framework. 

•  Developed  procedures  to  follow  desired  roots  from  one  iteration  to  the  next  and 
to  search  for  the  "shake-up"  states  which  are  important  to  the  interpretation 
of  photodetachment  spectra,  but  difficult  or  impossible  to  obtain  from  other 
coupled-cluster  based  methods  such  as  FSCC. 

•  Developed  formulae  for  transition  intensities  based  on  CCLRT  and  the  sudden 
approximation,  and  implemented  two  levels  of  approximation. 

•  Evaluated  the  predictive  abihty  of  EOMIP  by  applying  it  to  several  well-studied 
molecules,  and  compared  its  accuracy  to  that  of  ACC,  EE-EOM,  and  two- 
determinant  coupled-cluster  approaches. 
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•  Investigated  photodetachment  spectrum  of  the  linear  carbon  cluster  anions  C^ 
through  C^,  with  special  attention  to  its  ability  to  interpret  available  experi- 
mental spectra.  Various  basis  set  and  molecular  orbital  effects  were  observed. 

•  Demonstrated  the  formal  correspondence  between  the  EOMIP  and  FSIP  meth- 
ods, giving  simple  proofs  of  the  invariance  of  the  Fock-spax;e  method  to  the 
choice  of  active  space,  and  of  the  size  intensivity  of  the  EOMIP  method  for 
principal  ionizations. 

As  timings  in  table  12.1  indicate,  the  computational  the  required  for  converging 
these  12  roots  is  a  very  small  compared  to  the  time  required  for  a  CCSD  calculation. 
Only  3.2%  of  the  total  central  processing  unit  (CPU)  time  for  the  C^  job  is  spent 
on  the  EOMIP  calculations.  These  timings  are  comparable  to  those  for  the  FSCC 
calculation,  which  is  to  be  expected  since  the  formal  equivalence  of  the  two  methods 
means  that  the  terms  to  be  encoded  in  either  case  are  very  similar. 

From  essentially  single-reference  starting  points,  the  EOMIP  method  provides 
valence  ionization  energies  useful  for  comparison  to  experiment,  when  used  with  high 
quality  basis  sets.  Sufficient  polarization  functions  are  important  to  the  representation 
of  the  wavefunction,  but  addition  of  diffuse,  Rydberg-like  functions  to  these  basis 
sets  maJce  only  very  minor  contributions.  Using  pVTZ+sp  and  DZP-Hsp  basis  sets 
in  conjunction  with  the  2hlp  CCSD-EOMIP  method  for  small  closed  shell  molecules 
gives  quite  reasonable  results.  There  is  a  tendency  for  higher-energy  roots  to  show 
increasingly  positive  errors.  This  tendency  can  be  counteracted  to  some  extent  if 
core  hole  orbitals  are  used,  bringing  the  core  hole  peak  into  better  agreement  with 
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Table  12.1:  Computational  times  for  EOMIP  and  FSCC  calculations 

for  electron  detachment  of  C^  and  N2  using  the  PXTZ-f-sp  basis  set.  An  RHF  and 

a  UHF  reference  are  used  for  N2  and  C^,  respectively.    Timings  are  in  seconds  of 

Cray-YMP8 CPU time. 

Timings 

Computational  step  N2  C^ 

Integral  evaluation  24  409 

SCF  calculation  6  39 

Integral  transformation     19  236 

Integral  sorting  12  203 

CCSD  calculation  27  845 

Subtotal  89  1732 

FSIP  4     (13  iters,  4  roots)         21     (17  iters,  12  roots) 

EOMIP  singlet                     5     (12  iters,  4  roots)         33     (15  iters,  6  roots) 
EOMIP  triplet 24     (15  iters,  6  roots) 

experiment  as  might  be  expected,  but  also  improving  agreement  in  the  inner  and 
outer  valence  regions.  Using  smaller  basis  with  fewer  polarizing  functions  can  also 
counteract  this  tendency  for  the  outer  valence  and  core  peaks,  but  accuracy  in  the 
inner-valence  region  suffers.  The  shake-up  spectrum  is  poorly  represented,  with  too 
much  intensity  concentrated  in  regions  that  are  not  important  in  the  experimental 
spectrum. 

While  it  might  seem  reasonable  that  some  level  of  triples  approximation  would 
improve  the  position  of  the  higher  EOMIP  peaks,  perhaps  lowering  the  energies  too 
much,  as  is  the  case  with  the  SRCC  CCSD-fT(CCSD)  approximation  (and  less  so 
with  CCSD(T)),  an  examination  of  FSCC  triples  approximations  does  not  bear  this 
out.  In  Bernholdt's  work,  triples  corrections  increased  the  energies  more  often  than 
decreased  them.  Either  the  triples  approximations  are  poor  approximations  to  the 
true  triples  (as  is  indicated  by  the  size  of  fourth-order  triples  corrections  in  recent 
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work  by  Pal  and  Mukherjee[43]),  or  some  higher  order  of  excitation  is  as  important 
as  or  more  important  than  triples.  In  SRCC,  for  example,  the  quadruples  which  arise 
from  T2  are  often  more  important  to  the  wavefunction  than  the  triples  corrections, 
and  such  contributions  are  missing  in  the  EOM  or  FSCC  hnear  ansatz. 

For  cases  (such  as  that  of  a  multi-reference  starting  point)  where  double  excitations 
not  involving  the  annihilated  electron  have  very  different  magnitudes  in  the  N-1-  than 
in  the  N-electron  wavefunction,  the  method  might  be  improved  by  the  inclusion  of 
some  approximation  to  spectator  triples  in  the  EOMIP  manifold  space.  This  kind  of 
triples  correction  might  be  most  useful  in  cases  such  as  the  C2  ajid  C^  ionizations. 

One  source  of  difficulty  with  EOMIP  calculations  on  open-shell  systems  is  the  spin 
contamination  problem,  especially  for  small  molecules  where  the  sphtting  between 
states  with  different  S^  eigenvalues  is  largest.  There  are  two  approaches  which  can  be 
taken  to  this  problem.  One  is  simply  to  include  one  more  level  of  excitation  operator. 
For  example,  the  UHF  reference  for  C2  shows  (0|  S"^  |0)  =  .770,  with  an  average 
multiplicity  of  2.020.  The  CCSD  wavefunction,  because  of  the  exponential  form  of 
the  wave  reaction  operator,  contains  determinants  of  all  excitation  levels  possible  for 
the  system,  and  effectively  has  had  projected  out  of  it  the  worst  spin  contaminants. 
The  expectation  value  of  5^  is  .751,  with  the  approximate  spin  multiplicity  in  error  by 
only  .05%.  Adding  triple  excitations  to  the  EOMIP  wavefunction  would  improve  the 
spin-contamination  situation  by  pushing  the  problem  up  one  more  level  of  excitation, 
where  the  djoiamic  correlation  contributions  should  be  smaller.  This  is  a  brute  force 
approach,  one  which  significantly  increases  the  size  of  the  EOMIP  vector  while  adding 


many  new  terms  to  the  EOMIP  equations,  terms  which  require  more  operations  to 
evaluate  than  coupled-cluster  T2  terms,  although  fewer  than  coupled-cluster  T3  terms. 

The  spin  contamination  problem  can  also  be  solved,  somewhat  more  thriftily  in 
terms  of  computational  time  and  disk  space,  by  the  exphcit  spin  adaptation  of  the 
EOMIP  equations.  This  task  must  be  undertaken  separately  for  every  spin  case  of 
interest,  resulting  in  a  proliferation  of  code.  The  easiest  and  most  commonly  required 
case  is  the  open-shell  singlet.  Conceptually,  one  can  think  of  creating  two  types  of 
determinants  at  the  EOMIP  singles  level  -  one  based  on  removing  an  alpha  electron 
from  a  reference  with  one  unpaired  alpha  electron,  the  other  based  on  removing  a  beta 
electron  from  a  reference  with  one  unpaired  beta  electron.  Since  the  terms  involving 
only  those  operators  associated  with  one  or  the  other  of  the  references  are  identical 
to  what  is  already  coded,  the  only  portion  remaining  is  a  block  of  terms  involving 
an  alpha-beta  projection  of  H .  The  amount  of  work  involved  in  fully  implementing 
these  terms  is  similar  to  that  involved  in  the  implementation  of  the  original  terms, 
and  is  another  possible  course  for  future  development. 

Although  the  relevant  condition  for  the  EOMIP-SD  to  work  well  is  that  the  amph- 
tudes  obtained  from  Ti  *  S\  aren't  too  different  from  the  final  (Ti  *  5H-  52)  amplitude 
for  the  same  determinants,  it  seems  that  simply  having  smallish  TI  and  T2  ampli- 
tudes is  a  big  help  in  general,  something  that  can  be  most  fully  appreciated  after 
a  study  of  difficult  molecules  such  as  the  small  Hnear  carbon  cluster  anions  studied 
here.  No  doubt  important  reasons  for  this  are 


•  The  method  is  most  applicable  when  the  S3  and  higher  amplitudes  are  pretty 
nearly  small  enough  to  neglect  and 

•  large  Ti  and  T2  amplitudes  can  create  large  contributions  at  the  S3  and  higher 
level  (as  Ti  *  S2,  T2*  Si,  T2  *  S2,  etc.).  in  contradiction  to  the  first  point 

There  might  be  exceptional  cases  where  a  large  Tj  *  52  or  T2  *  Si  contribution  to 
the  ^s  part  of  the  wavefunction  is  just  exactly  the  one  needed,  or  an  erroneous  53  con- 
tribution to  the  wavefunction  is  fortuitously  in  just  the  right  direction  energetically, 
but  these  cases  are  most  likely  exceptions  to  the  rule. 

The  coupled- cluster  EOMIP-SD  has  been  shown  useful  for  producing  ionization 
spectra  which  compare  well  with  experiment  for  inner  and  outer  valence  regions,  in- 
cluding inner  valence  shake-up  effects,  as  well  as  core  hole  primairy  ionizations,  for 
both  saturated  and  unsaturated  closed  shell  molecules  when  high-quahty  basis  sets 
are  used.   For  the  hnear  carbon  cluster  anions,  all  of  which  are  open  shell,  most  of 

Table  12.2:  EOMIP  errors,  pVTZ-h  basis,  small  closed-shell  molecules 


first  outer        inner 

ionization     valence     valence         core 


closed  shell 
core  hole 

0.2 
0.1 

0.5 
0.3 

1.0  to  1.5 
0.5  to  1.3 

0.8 
0.2 

which  are  characterized  by  partially- filled  it  orbitals,  the  situation  is  different.  The 
EOMIP  approach  has  several  advantages.  First,  it  can  be  used  to  directly  calculate 
many  ionizations  quickly  and  cheaply  in  terms  of  computer  time,  as  compared  to 
single-state  methods  which  require  an  expensive  correlated  calculation  for  each  state. 
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Second,  EOMIP  can  be  used  to  search  the  shake-up  space.  It  is  difficult  or  impossible 
to  converge  on  desired  shake-up  peaks  using  conventional  N-electron  methods,  and 
the  use  of  the  well- correlated  CCSD  reference  state  makes  EOMIP  more  accurate 
than  the  2hlp  CI  approach.  The  Davidson  diagonalization  method  for  the  extraction 
of  roots,  along  with  the  root  finding  and  root  following  algorithms,  makes  the  EOMIP 
approach  less  susceptible  to  the  intruder  state  problems  which  can  make  it  difficult  to 
impossible  to  converge  inner  valence  FSIP  states.  While  highly  correlated  calculations 
such  as  CCSD(T)  will  provide  better  information  about  specific  roots,  EOMIP  gener- 
ally provides  more  information  overall  about  more  portions  of  the  ionization  spectrum. 
EOMIP  can  be  reliably  applied  to  systems  with  conventional  bonding  characteristics. 
More  experience  with  unconventional  bonding  situations  such  as  the  even-length  lin- 
ear carbon  cluster  chains  should  shed  some  light  on  which  types  of  roots  are  reliable 
in  these  situations,  and  what  can  be  done  to  improve  EOMIP  performance.  Co- 
development  of  triples  approximations  from  FSCC  and  EOMIP  standpoints  should 
lead  to  better  performance  of  both  methods  in  most  situations. 


APPENDIX  A 
INVARIANCE  OF  FOCK-SPACE  ENERGIES  TO  CHOICE  OF  ACTIVE  SPACE 

Invariance  of  the  Ionization  Energies 

Prom  the  definitions  given  in  section  4.1,  the  FSIP  equations  have  been  shown  to 
satisfy 

Q[HN^-m^ff)P    =    Q(^H^  +  [Hr,S)^-(SH^ff)^y  =  0  (A.l) 

where   Pi/g_^P    =    P[Hj,  +  {Hs)c- {SH!,),)P  ^  P^Hr,  +  {Hj^S),)^.2) 

The  operator  equation  PH^txPC  =  PCW,  as  well  as  the  identity  QH^rcPC  —  0  can 
be  written  in  matrix  form: 

(HNn-QH^ff)^^    =    {HQp  +  HQQSQp-SQpH^ff)Cpp  =  0         (A.3) 

=    CppW  (A.4) 

Consider  a  partitioning  of  the  active  space  into  A  and  its  complement  in  P,  B 
{=  P  -  A).  The  complement  to  A  in  the  full  (P  +  Q)  space  is  Q'  {=  Q  +  B).  The 
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matrices  appearing  in  the  H^a  expression  can  be  partitioned  into  A  and  B  parts, 


Hpp 


Haa     Hab  \  fi       _  (  ^AQ 


Hba    Hbb 


Sqp  —  (^  Sqa    Sqb  j      Cpp  —  \   n 


Caa    Cab 

A      CbB 


w  = 


Wa      0 

0    y\!B 


etc. 


which  leads  to 


H 


eff 


Haa    Hab    I   I   (    ^aqSqa    HaqSqb 
V  Hba    Hbb  j      \  HbqSqa    HbqSqb 

CA.4W4    CabWb  \ 


Caa    Cab 
Cba    Cbb 


(A.5) 


The  AA  part  of  this  equation  can  be  written  as 


{(^eff)^^  +  {tieff)^B^BACA\)CAA  = 

[Haa  +  HabCbaCaa  +  HaqSqa  +  HaqSqbCbaCaajCaa  —  Caa^^a,  (A. 6) 
and  the  BA  part  as 

[Hba  +  HbqSqa  +  HbbCbaCaa  +  ^bqSqbCbaC^ajCaa  =  C'b^VV.4.(A.7) 

The  Bloch  equation  can  also  be  expanded: 
Q[Hn-QH^ff)PCpp  =  0 

[{{h^^)qa  m^)QB)  +  {i^^H^effhA  mtff)QB)]( '^^■^  ^'^^ 


Cba    Cbb 


0; 


which  gives,  in  the  AA  block  of  the  product. 


[Hqa  +  BqbCbaCaa  +  HqqSqa  +  HqqSqbCbaCaa 

-SqaH^^^^  -  SqaH^^^JCbaCa'a 
-SQBH^ff^^  -  SQBH^ff^^CBAC-A\j  Caa 
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=      [fiQA  +  ^QqSqa  +  tiQBSsA  ~  SqACAAy^ACAA  ~  SqbCbaV^aCaa)  ^ AA 

=    0  (A.8) 

Consider  now  the  solution  to  the  FS  equations  in  the  A  model  space 

[Hq'a+ Hq'q'Sq,a- Sq,aH-^jjjCaa    =    0  (A.9) 

H^effCiA={RAA  +  HAQ'S^p)ci^    =    Ci,W^.  (A.IO) 

It  will  be  shown  that  the  solution  to  the  FSIP  equations  in  the  smaller  space  is 

C^A    =    Caa  (A.ll) 

S^SA    ^    CbaCa\  (A.12) 

^QA      —      \^QA  +  SqbCbaC'aa)  —  ^QA  +  SqbSqa  (A. 13) 

W^    =    W4.  (A.14) 

In  order  to  do  this,  two  conditions  must  be  satisfied: 

1)  H'^ff  =  Haa  +  HabS^a  +  HaqS^q  =  CaaWaCaa;  (A.15) 
and 

2)  (m^-Q^H^ff)       =0.  (A.16) 

The  satisfaction  of  the  first  condition  relies  on  nothing  more  than  the  fact  that  H 
does  not  depend  on  the  choice  of  active  space.  Expanding  A(Q^~  HujA  gives 

-^gff    =    Haa  +  Haq'S'q,a 

=    Haa  +  HaqS'qa  +  HabSba  (A- 17) 

=     i"eff)^^+{Heff)^sS^A 
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Substitution  of  the  definitions  given  in  equations  A.  11  through  A.  14  into  equation 
A. 18  gives  exactly  the  parenthetical  expression  in  equation  A. 6. 

For  the  proof  that  the  second  condition  is  satisfied,  the  QA  portion  of  equation 
A. 15  is  compared  to  equation  A. 8,  and  the  BA  portion  is  shown  to  be  zero. 

=     {Sqa  +  HqqSq^  +  HqBSgji  -  Sq^H^^jCaa 

—  [Hqa  +  HqqSq^  +  HqbSqj^  -  SqACAAy^ACAA  "~  'S'qbC'b^C^^(^CaaW^C'X4 j  JC'aA 

—  \Hqa-\-  HqqSqji^  +  HqbSqj^  —  SQACAAy^ACAA~  ^QBCsAy^ACAAJ^AA  (A. 18) 

Since  this  is  exactly  the  same  expression  as  that  given  in  equation  A. 8,  the  QA 
projection  of  the  Bloch  equation  is  satisfied.  Finally  using  equation  A. 7  as  well  as 
the  the  definitions  in  equations  A.  11  through  A.  14, 

=    [Hba  +  SbqSq^  +  HbbSba  ~  ^BA^ejf)^^'^ 

=      \CBAy^ACAA~  ^BACAAipAAy^ACAAjJ^AA 

=    0,  (A.19) 

the  BA  projection  is  shown  to  be  zero  as  well.  Therefore,  the  definitions  in  equations 
A.  11  through  A.  14  give  the  solution  to  the  Fock-space  equations  in  the  smaller  model 
space  A  and  the  energies  are  invariant  to  a  reduction  in  the  size  of  the  model  space. 
Any  two  model  spaces  for  the  IP  problem,  whether  or  not  one  is  a  subset  of  the  other, 
are  subsets  of  the  "complete  active  space" .  The  FSIP  energies  obtained  in  either  case 
are  identical  to  those  obtained  for  the  same  model  functions  in  the  complete  active 


space,  and  therefore,  if  there  are  any  model  functions  in  common,  to  those  obtained 
in  each  of  the  smaller  model  spaces. 

A  Modification  of  the  Bloch  Equation  to  Preserve  Invariance  in  Any  Sector 

In  the  higher  sectors,  fi^"'-"'  (=  {e^'""'"'}  )  and  C^'"-")  give  the  solution  to  the 
pim.n)  jnodel  space  problem,  while  Q^"?"?  is  the  product  of  all  J^^'-*^^  where  l+k  <  m+n. 
In  the  discussion  which  follows,  all  the  {m,n)  superscripts  will  be  dropped.  The 
operators  fi,  f^^/d'  ■^'  Q'  ^^^  ^  ^^^^  ^^^  ^^  assumed  to  be  those  which  correspond  to 
the  sector  being  discussed.  Hf^  is  the  familiar  Q'^°'°'>~  HnQ^°-^'>  —  e~'^Hse^.  Prom 
the  definitions  given  above,  it  can  be  shown  that 

PH^ffP  =  P{R  +  [RS],  -  {SR)c)P  =  P{R  +  iRS)c)P,  (A.20) 

where 

R  =  H^^Qm. 

The  equation  PH^aPC  =  PCW,  as  well  as  the  identity  QH^^PC  =  0  can  now  be 
written  in  matrix  form: 

PH^ffPC    =>    f  Rpp  +  RpqSqp  j  Cpp 

=    CppW  (A.21) 


where 


Rpp    =    Hpp  +  HpQin^i^ -  Ij^^ 

RpQ         =        HpQ 
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Since  the  S  operators,  by  definition,  can  not  cause  scattering  into  the  P  space,  the 
P  projection  on  the  left  prevents  any  contribution  from  f2~^  The  right  projection 
by  P  and  the  normal  ordering  of  ^I^qi  prevents  many  products  of  S  amplitudes  from 
contributing.  The  surviving  contributions  from  ^tot  ^^  e^(l  +  S  +  i^old  ~  ^"^r 

The  model  space  spanned  by  P  is  defined  by  the  choice  of  active  orbitals.  The 
usual  development  of  the  FS  equations  assumes  that  this  same  choice  is  used  for  all 
sectors.  Suppose,  however,  that  it  were  desirable  to  solve  the  lower-rank  sectors  in  one 
active  space,  and  the  current  one  in  a  smaller  model  space  than  that  defined  by  the 
original  choice  of  active  orbitals.  The  following  discussion  will  show  that  it  is  possible 
in  this  smaller  model  space,  by  using  all  the  S  amplitudes  from  the  lower-ranked 
sectors,  to  arrive  at  the  same  energies. 

The  new  model  space  is  a  subset  of  the  P  model  space,  called  A;  its  complement 
in  F  is  B  {=  P  —  A).  The  complement  to  A  contains  both  Q  and  B,  and  will  be 
referred  to  as  Q'.  The  matrices  appearing  in  the  H^a  expression  can  be  partitioned 
into  A  and  B  parts, 

RaA      RaB   \  D    _  /  ^AQ    \ 


Rpp  — 


W  = 


Rba    Rbb 

Sqa    Si 

yvA    0 

0   H' 


R 


^^-   i? 


Sqp  —  [^  Sqa    Sqb  j      Cpp  —  (  ^ 


BQ   / 

Caa    Cab 
Cba    Cbb 


b 


etc. 


which  leads  to 


H 


eff 


(  Raa    Rab   ]   I   [   RaqSqa    RaqSqb 
{  Rba    Rbb  J      \  RbqSqa    RbqSqb 

CaaV^a    Cab^^b  \ 


Caa    Cab 
Cba    Cbb 


(A.22) 


The  AA  part  of  this  equation  can  be  written  as 


{{H^ff)AA  +  {H^ff)ABCBACAA-')CAA     =  (A.23) 

yRAA  + RabCbaCaa  + RaqSqa  + RaqSqbCbaCaaJ^aa    =    C'aaWa-(A.24) 
If  we  define 

S'ba    =    CbaCaa  (A.25) 

and   5-;^^    =    {Sqa^SqbCbaCa\),  (A.26) 

it  is  possible  to  simplify  the  expression  to 

[Raa  +  RabS'ba  +  RaqS'qajCaa  =  Caa^a-  (A. 27) 

Similarly,  the  BA  portion  of  this  equation  becomes 

{H^ff)BA  +  {H^jj)bbCba  = 

[Rba  +  RbbS'ba  +  RbqS'qajCaa  =  Cba^^a-  (A.28) 

If  Q'  =  {exp(5')}  is  defined,  the  AA  expression  has  the  form  of  an  effective 
Hamiltonian,  and  can  be  written 

A{Rn')^C'A  =  C'AWa  (A.29) 

If  we  define 

Kff    =    {(^r'R^%^  (A.30) 

=      {H^ff)AA  +  {H^ff)ABCBAC'AA  =  CAAWACX'A  (A.31) 

=    Haa  +  Haq{QoI(1  -  i)q^  +  [^AB  +  Haq{^oI(1-  ^)qb)S'ba  +  HaqS'qa 


{H'^ff)BA  =  {{nr'Rn^^  (A.32) 

=    {H^jj)ba  +  {H^ff)BBCBAC-A  =  Cba\^aC-A  (A.33) 

=    Hba  +  Heqipoid  "  l)g^  +  [^bb  +  Hsqi^old  ~  ^)qb)S'ba  +  ^bqS'qa 

we  can  develop  a  method  which  is  invariant  to  this  reduction  of  the  model  space. 

A  new  type  of  term  has  been  introduced  in  the  H'  rr  expression  (see  the  underhned 
term  above).  The  {H'  ff)BA  term  is  created  to  appear  in  the  modified  FS  equation, 
which  is 

Q'(^Hr,n  -  [n^i^  -  l){H'^ff)BA  -  ^Kff)^  =  0,  (A.34) 

and  breaks  down  into 

Rqa  +  RqbS'ba  +  RqqS'qa  -  i^old  ~  ^)qb{H'^jj)ba 

-{^old-^)QAKff-SQAKff    =    0  (A-35) 

and 

Rba  +  RbbS'ba  +  RbqSqa  —  ^BA^eff    ~    ^'  (A. 36) 

More  expUcitly  with  all  the  "new"  terms  (excepting  those  in  H'  rr)  underhned  this 
becomes 

Hqa  +  Hqq(Q.q^  -1)qa-V  [Hqb  +  Hqq{Q.qI^  -  1)qb)S'ba 

+HqqS'qa  -  i^old  ~  ^)QB{H'^ff)BA 

-i^old  -  ^)QAH',ff  -  S'QAH',ff  =  0  (A.37) 

and 

Hba  +  ^BQi^oid  -  1)qa  +  {Hbb  +  HBgi^pid  -  1)qb)'5'ba 
+HbqS'qa  -  S'qaH'^^  =  0.  (A.38) 


It  is  relatively  easy  to  show  that  the  5^^  and  5q^  defined  above  satisfy  these 
modified  Bloch  equations.  Returning  to  the  Bloch  equation  in  the  larger  space, 

Q{H^tot-^totHeff)P    =    0 
When  partitioned,  the  Q-A  block  is 

[RqA  +  RqqSqa  -  i^old  -  ^)QB{HQff)BA 

-{n^l^  -  l)QA{H,ff)AA  -  Sqa{H^^)aA  -  SQB{H,ff)BA}  =  0       (A.39) 

and  the  Q-B  block  is 

[RqB  +  RqqSqb  -  i^old  ~  ^)QB{Hgff)BB 

-{n^l^  -  l)QA{H,ff)AB  -  Sqa{H^jj)aB  -  SQB{H,ff)BB}  =  0.      (A.40) 

Multiplying  the  Q-B  block  by  CbaCaa  ^^'^  adding  the  two  expressions  yields 
\Rqa  +  RqbCbaCji^js^  4-  Rqq  \Sqa  +  SqbCbaCaa\ 
~^^old~  ^)qb  [{H(^ff)BA  -  {H^jj)bbCbaCaa[ 

~(^0ld'  ^)«-4  l^^eff')AA  -  {H^ff)ABCBACAA] 


-s, 


QA 


{H^ff)AA  +  {H^g)ABCBACj^\ 


-Sqb    {H^g)BA  +  {Hf,ff)BBCBACAA    }  =  0 


(A.41) 


This  expression  simplifies  to 


\Rqa  +  RqbS'ba  +  RqqSqa 
~i^old~  ^)QB{ff'eff)BA 
-(^old-^hAH'^ff 

-SqACAAy^ACAA  ~  SqBCBAy^ACAAj  —  0- 


The  last  two  terms  reduce  to 


which  allows  the  original  Bloch  equation  to  be  expressed: 

|i?Q.4  +  RqbS'ba  +  RqqS'qj^  -  {^old  ~  '^)qb{H'^jj)ba 

This  is  identical  to  the  modified  Q-A  projection  of  the  Bloch  equation  expression  in 
the  reduced  space. 

The  B-A  part  of  the  Bloch  equation  can  be  simplified  to 

{H'^ff)BA  -  S'baH'qjJ 
=    CBAyVAC2\-iSBACA)CAA'^AC2\  =  0  (A.43) 

Therefore,  the  S'b^  and  S'q^  satisfy  the  modified  Bloch  equations,  and  define  a 
modified  effective  Hamiltonian,  H'  rr.  The  eigenvalues  are  identical  to  a  subset  of 
those  in  the  larger  model  space  and  H^fc  eigenvectors  are  the  projection  into  the 
smaller  model  space  of  those  from  the  original,  larger  model  space. 

The  Bloch  Equation  and  Modified  Bloch  Equation 

To  make  clear  the  differences  between  the  equations  used  in  an  implementation 
of  the  original  and  modified  Fock-space  method  that  has  been  presented  here,  the 
equations  for  both  will  be  presented  in  as  much  detail  as  possible. 


Equations  in  the  Original  Model  Space 
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The  P-P  projection  of  the  effective  Hamiltonian  is: 


{^eff)pP  ~  ^PP  +  ^PQi^old  ~  ^)qp  '^  ^PQ^QPr 


and  the  equations  which  the  Sqp  amplitudes  must  satisfy  are 


Hqp  +  Hqq{9.q^  -  l)gp  +  HqqSqp 

-i^old  ~  ^)qp  [^pp  +  ^PQi^old  -  ^)qp  +  ^^pq^Qp] 
-SQp{H^ff), 


(A.44) 


'  pp 


(A.45) 


Equations  in  the  Reduced  Model  Space 


The  A-A  projection  of  the  modified  effective  Hamiltonian  is 


H'^fr  =  Haa  +  Haq  ip-old  ~^)oa'^  i^^^  "^  ^AQi^old~  ^)qbjS'ba  +  HaqS'qa- 


'QA 


(A.46) 


The  equations  which  the  S'q^  and  Sg^  amplitudes  must  satisfy  axe 


Hqa  +  Hqq[QqI^  -  1)qa  +  HqqSqp  +    Hqb  +  HqqJQ^i^  -  1)qB 


9' 


BA 


~(^old  ~  ^hp  Hba  +  Hsqi^old  ~  ^)qa  +  ^bqS'qa 


+  [Hbb  +  ^Bqi^oid  -  ^)qb}S'ba\ 

'^^old  ~  ^hA  [Haa  +  HAQi^old  ~  ^)qa  +  ^aqSqa 
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+  [Hab  +  HAQJ^oid  -1)qb  \S'ba 


-^QA^'eff 


(A.47) 


and 


Hba  +  Hsqiytold  ~  ^)qa  +  ^bqS'qa  +   Hbb  +  ^Bqi^pid-  ^)qb\  S'ba 


-^BA^ff 


(A.48) 


where  the  terms  that  would  not  arise  in  the  usual  treatment  have  all  been  underhned 


for  clarity. 


APPENDIX  B 
THE  IMPORTANCE  OF  TRIPLES  -  CC  EOMIP-SDT  PERTURBATION  THEORY 

As  Mukherjee  et  al.  have  pointed  out [10],  the  CC  EOMIP  equations  for  the  prin- 
cipal ionizations  are  equivalent  to  the  Fock  Space  IP  equations.  With  a  convenient, 
connected  perturbation  theory  thus  defined,  it  is  a  simple  matter  to  use  this  pertur- 
bation theory  to  examine  the  importance  of  a  triples  correction  to  the  CC  EOMIP-SD 
theory  presented  and  used  in  this  thesis. 

The  CC  EOMIP-SDT  Equations 

The  full  CC  EOMIP-SDT  equations  can  be  written  as 

M  MA  MMA 


MVST 


MB  M 


^  N'M  M 

-IT.(i:  {^M I \B£)  t'A  c^^  +  E  (^1  /  \^)  ^tYv 
^  mmb\  s  I  vr 

-  ^(1  -  Pxj)Y.  {m\J^)4irv  +  \  E  {Anej^)^%       (B.2) 


2 
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IC^ICAB      _ 


UJ    c 


■XJK. 


- 1:  {MAB\]\IJ)C)cl,  +  ^  E  (E  {MM\\SC)fc^)c'^f, 

M  ^  MNS\  C  / 

+  AxjK  E  {AB\\eK)4'j  -  (1  -  Pab)  Axjk  E  {iB\\JK)c^^ 

e  M 

-  Axjk  E  (-^I  / 1^)  ^MfK  +  (1  -  Pab)  E  (^I  /  \^)  ^xYk 

M  T 

+  \  E  {AB\\E7)cYf^  +  \a^jk  E  {mM\\1J)cI,^k 

+  (1  -  Pab) Axjk  E  (P^ri^^)^^:/^  +  \Mjk  E  (^7^B||  |£: J-X:)c^§; 

PF  ^  PEF 

-  \{\  -  P^b)  AxjkY:  {MVB\\\JKT)c^ifT,.  (B.3) 


NPF 


For  the  purpose  of  developing  an  EOMIP  perturbation  theory,  it  is  convenient  to 
write  this  in  a  compact  matrix  notation, 


M.MVZT 


(B.4) 


mnvst  ^-^^ 


c'm+  E 


yl 


>1 


^ 


>fA^)'^^ 


MATT')'^^^^ 


(B.5) 


,  .KJCAB 


I    AB\ 


M 


I    AB 


/  AB  .     ey=\KeT 


MM 


^mm 


MNVST   ^ 


XJK 


J^J^T.j^'MNV 


(B.6) 


A  review  of  the  Ravleigh-Schrodinger  Perturbation  Theory 


In  a  nutshell,  Rayleigh-Schrodinger  perturbation  theory  (RSPT)  is: 


{H-E)\I)    =    0 


{h°  -  E') 


7°      =    0 


H    =    H^  +  H 


0    ,     t/l 


{/°|/°)  =  1,    (/"|/)  =  1,    (/°|/^)  =  0(iV^O) 


E    =    E^'  +  E' +  E^  +  E^  +  ■■■ 
{{E°  -  H")  +  iE'  -H')  +  E'  +  E'  +  ---)  (|/°)  +  |/')  +  |/')  +  |/')  +  •••)=  0 

In  order  for  the  final  equahty  above  to  hold  true  the  equahty, 

{E'  -  H')  |/^)  +  {E'  -  H')  j/'^-i)  +  Y^E""-'  \l')  =  0, 

fc=0 

must  hold  true  for  all  A'',  A'^  =  1,  oo.  This  means  that 
(/°|i/^|/°)-(/°|£^|/°)    =    0 

{i'\e'\i''-')  =  {i'\h'\i''-')-{i'\J2e''-'\i') 


jN^ 


=    (E'  -  H')-'    {H'  -  E')  /^-^)  -  Y.  E"-'  I/') 

\  Jk=0  / 

=    (£0  -  H'r'  ({H'  -  E')  l/'^-i)  -  |/°)  £^  -  Ye"-'  \i')^ 

=    {E'  -  H")-'  ({H'  -  E')  I""-')  -  /°)  (/°  {H'  -  E')  /^-i)  -  Y  E""'"  |/') 

=    [E'  -  H')-'  (r'{H'  -  E')  |/^-^)  -  YE""-'  |/')] 


N-2 


k=\ 


EQMIP  Perturbation 


It  is  conceptually  simplest  to  choose  a  Koopmans'  or  extended-Koopmans'  zeroth- 
order  wavefunction,  and  a  zeroth-order  array  which  is  diagonal  in  these  single-ionization 
wavefunctions.  This  means  choosing  ^°  to  be  either  {i\  f  \m)  Sim  or  (i|  /  1^)  ^im-  If 
the  latter  is  chosen,  the  first  order  energy  will  disappear.  The  author  has  chosen  the 
former.  The  quantities  to  be  used  in  the  perturbative  development  are: 

A'    =    {i\f\m)5,m  {i\A'\m)    =    {i\  f  \m)  -  {i\  f  \i)  S,m 


{Ky\  =  x:|o) 

u'^'    =    {K'\f\K') 
In  the  first  order, 


R^    =    1  -  P'^ 


KQ  A<i\  uK. 


Z^    =    {uj^''-A')R^. 


(B.7) 


u'''    = 


C-'    = 


{i\A'\)C)/Af 


'I  J 


"i  jk 


ab 
ijk 


A' 
A' 


^)l^! 


x:  /Agf , 


where  Af  =  u'^^  +  fa,  Af/  =  u^^  +  U  +  f^^  -  f^a,  etc. 
In  the  second  order, 


cu 


.a:2 


}C' 


1  Die 


A'R 


K') 


I  ija    \  J  I 

.E{'C|.4l.;:)cf,f.  (B.S) 


,Jk 


The  final  term  of  equation  B.S  represents  the  first  triples  correction  to  the  EOMIP-SD 
energy.  In  a  perturbation  analysis  of  the  CCSDT  or  CISDT  equations,  there  would 
be  no  triples  contribution  to  the  energy  until  4th  order.  Because  H  is  not  hmited  to 
1-  and  2- body  terms,  the  EOMIP  triples  amplitude  appears  in  order  1.  And  because 
only  one  {ijWab)  term  is  required  to  close  the  energy  term  (as  opposed  to  one  (ii||a6) 
and  one  {ij\\kb)  or  {cj\\ah)  term  for  CC  or  CI),  this  triples  term  appears  in  the  energy 
in  order  2.  In  more  detail,  this  term  is: 

=    {Mef)  (-  T.  (jk\\ai){Kmm)tZ'  +  E  W\ab}{Ka\\ICc)t^J; 

\        m  e 

-  Y:  {jk\\ab){K:b\\ck)t^^  +  •  •  -^  /A^f  *  S)Cz  (B.9) 

The  lowest  contributions  to  the  T"^  M^l  ^)  term  are  of  order  2  in  W,  which  makes 
this  energy  term,  of  order  2  in  yl,  minimally  of  order  3  in  W.  This  is  still  an  order  lower 
than  the  first  triples  contribution  appearing  in  CCSD  or  CI,  implying  that  triples  are 
more  important  in  the  EOMIP  problem  than  they  are  in  the  CC  problem. 
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